UNCLAr:j;iH'iKD 


?JIT-Plli7  (5 


DEPARTMENT 
OF  PHYSICS 

DDC 


p  NOVI  1965 
ilibv^UatJ  1/ 1 


'ICGATHON  ACCFLP’.HATOR  FINAL  REPORT 
-ipril  15,  1961  -  June  15,  19bli 


SIGNAL  CORPS  CONTRACT  NR.  DA  36-n39-sc-872i42(E) 
/vRPA  Order  Nr.  112-6?,  Project  Code  Nr.  3720 


Hnited  iltates  Anay 
Electronics  Materiel  Agency 
Fort  fKnmouth,  New  Jersey 


BEST 

AVAILABLE  COPY 


SIT-Fll+T  (5/6'^) 


Information  for  Abstract  Card 
Department  of  Physics 
'.evens  Institute  of  Technology 
Hoboken,  New  Jersey 


Megatron  Accelerator  Final  Report 

K.  C.  Rogers,  S.  J.  Lukasik,  L.  Ferrari 

Final  Report,  April  15,  1961  -  June  .5,  196** 

223  pp  illustrated 

Signal  Corps  Contract  No.  DA  36-039-3c-872**2(e) 
Unclassified 

ARPA  Order  Nr.  112-62  Project  Nr.  3720 


This  is  the  final  report  on  the  Megatron  (plasma  betatron)  Project, 
Construction  details  of  the  accelerator  are  given,  and  methods  of  shaping 
and  measuring  the  magnetic  guide  field  are  described.  Methods  of  generating 
K  j  plasma  and  its  initial  properties  are  listed.  Application  of  the  time- 
varying  magnetic  guide  field  results  in  the  production  of  a  runaway  electron 
beam  of  peak  energy  1  Mev  and  beam  intensity  approximately  one  ampere. 
Methods  of  detecting  the  runaway  electrons  are  descrioed,  and  the  electron 
life  histories  are  studied.  It  is  found  that  ultimate  destruction  of  the 
runaway  electron  beam  results  from  a  collective  interaction  of  the  beam 
electrons  aind  the  background  plasma.  The  observations  are  compared  with 
theoretical  results  of  Field  and  Fried.  Also  included  are  the  results 
of  a  computer  study  of  early  time  collective  effects  in  a  cold  plasma 
subjected  :o  a  rapialy  rising  betatron  guide  field. 


TABLE  OF  CONTENTS 


Introduction  ,  "A  Survey  of  the  Megatron  Project" 
K.  C.  Rogers 


"Behavior  of  Runaway  Electrons  in  a  Plasniii  Betatron" 
L.  A.  Ferrari,  K.  C.  Rogers 


"Early-Time  Collective  Effects  in  a  Plasma  Bete.tron" 
S.  J.  Lukasik,  K.  C.  Rogers 


Papers  and  Publications  Produced  on  the  Megatron 


Personnel 


Distribution  List 


A  SURVEY  OF  THE  HBQATRON  PROJECT 


j. 

In  1956, ac.  thn  CERN  Accelerator  Conference  ,  Q,  I,  Budker  suggested 
that  the  acceleration  of  electrons  within  a  plasma  by  means  of  magnetic 
inducticn  may  be  a  way  to  create  an  intense  (kiloampere)  beam  of  rela¬ 
tivistic  electrons.  His  suggestion  was  accompanied  by  a  preliminary  re¬ 
port  on  a  betat;,*on  device  in  which  he  claimed  to  have  generated  relativistic 
electron  beams  of  approximately  20  amperes.  The  following  year  D.  Flnkelsteln 
ax  SteveiiJ  in  u  brief  i,;ote  outlined  a  mega-gauss  betatron  (Megatron)  scheme 

for  production  of  relativistic  electrons  in  a  plasma.  The  Megatron  design 

3 

was  guided  by  the  earlier  work  of  Forth,  Levine  and  Wannlek  on  the  pro¬ 
duction  of  ve:ry  intense  pulsed  magnetic  fields.  At  about  the  same  time 
(195?)  Flnkelsteln  at  Stevens  initiated  the  construction  of  a  low  inductance 
rapid  discharge  capacitor  bank  (100  K  Joule)  suitable  for  driving  an  air-core 
betatron  flela  ceil.  He  was  Joined  by  K,  G.  Rogers  in  the  work  In  early  1958 
and  a  proposal  for  support  cf  the  work  was  submitted  to  the  U,  P,  Atomic 
Energ}’’  Commission,  Construction  of  the  bank  proceeded  in  1958  under  the  di¬ 
rection  of  K.  C,  Rogers  while  D,  Flnkelsteln  spent  a  year  at  CERN,  In  the 
Spring  of  1959  the  Atomic  Energy  Commission  awarded  a  one  year  contract  to 
Stevens  for  construction  of  the  plasma  betatron  accelerator.  Work  carried 

out  under  this  contract  was  reported  at  the  1959  CERN  Accelerator  Conference 

'\ 

and  aroused  considerable  interest.  Further  support  of  the  Stevens  work  was 
requested  of  the  Advanced  Research  Projects  Agency  and  was  awarded  through  a 
plasma  physics  contract  at  Stevens  with  the  United  States  Army  SJ.gnal.  Corps, 
This  agency  took  over  support  of  the  project  until  its  termination  in  June 
196i|,  Dr,  Flnkelsteln  joined  the  faculty  of  Yeshiva  University  In  I960 


and  terminated  his  connection  vriLth  the  project  in  i96l.  Since  that  time 
the  work  has  been  carried  out  under  the  direction  of  K.  C.  Rogers. 

Many  technical  problems  arose  during  the  course  of  the  project  and 
there  were  serious  delays  until  they  could  be  solved.  The  original  scheme 
for  generating  the  betatron  field  proved  unworkable  because  of  difficulties 
in  producing  the  desired  field  shape.  Another  method  was  found,  but  this 
placed  limitations  on  the  maximum  volts/tum  that  could  be  .applied  to  the 
plasma.  Construction  of  the  low-inductance  high  voltage  capacitor  bank 
was  difficult  and  s.low  because  of  dielectric  breakdown  prob'.  ams,  The  most 
difficult  and  time  consuming  task  was  that  of  shaping  the  air-core  betatron 
field  so  as  to  have  all  of  the  desired  properties  l.e»,  suitable  vector  po¬ 
tential  well  depth,  freedom  from  azimuthal  inhomogeneities,  etc.  The  shape 
of  pulsed  magnetic  fields  varying  on  a  microsecond  time  scale  had  to  be 
measured  with  a  precision  of  ~  0.1^.  The  production  of  large  runaway  electron 

beams  in  the  plasma  required  that  the  betatron  field  be  applied  to  a  low 
10  -3 

density  10  cm  '^)  plasma  \d.th  no  appreciable  neutral  particle  background 
density.  Production  of  a  suitable  plasma  in  the  acceleration  chamber  re¬ 
quired  special  techniques  which  in  turn  introduced  problems  of  their  own 
regarding  betatron  guide  field  perturbations.  All  of  these  probler.is  were 
eventually  solved  and  runaway  electron  beams  were  consistently  produced  with 
energies  approaching  the  peak  energy  of  the  machine,  approximately  1  Kev. 

(The  machine  peak  energy  was  set  by  the  dielectric  strength  of  the  flux 
guides  in  the  betatron  field  coil. )  However,  they  were  disappointingly 
smiall  (~  1  ampere).  In  spite  of  the  considerable  care  taken  to  optimize 
the  conditions  for  runaway  electron  production  the  beam  Intensities  produced 
remained  siaall.  The  original  design  goal  of  the  project  was  the  production 

-ii- 
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of  klloawpere  electro'"  beans  with  applied  eleotrio  fields  of  2  x  10'  voits/cn. 
Dielectrio  strength  United  the  applied  electric  fields  in  our  machine  to  ~ 

20  "Tolts/om  and  resulted  in  beaias  of  approximately  one  ampiere.  Even  if  our 
results  could  be  oxtrapodated  to  the  original  design  figure  of  2  x  lO"^  volts/ cm, 
we  would  expect  only  100  amperes  of  beam.  The  technical  difficulties  la  ac¬ 
hieving  kilovolt/ cm  applied  electric  fields  in  a  betatron  guide  field  are  so 
great  that  we  do  not  consider  it  reasonable  to  attempt  to  meet  them. 

Our  results  further  show  that  the  runaway  electron  beams  created  iri  the 
plasma  eventually  are  destroyed  by  collective  effects.  The  beam  survives  to 
nearly  the  end  of  the  acceleration  cycle  before  its  destruction  however,  so 
that  the  electrons  are  accelerated  to  approximately  the  peak  energy  of  the 
machine.  We  conclude  that  our  results  offer  little  encouragement  for  further 
development  of  plasma  betatrons  as  sources  of  intense  relativistic  electron 
beams. 

Many  colleagues  and  students  have  contributed  from  time  to  time  to  this 
project.  This  final  report  was  prepared  by  L.  Ferrari,  S,  J,  Liuvcsik,  and 
K,  C,  Rogers,  but  it  rests  heavily  on  the  many  contributions  of  Q*  Brucker, 

J,  R,  M,  Coulter,  D,  Finkelstein,  H,  Huber,  A,  Jeraiakian,  C,  T,  Lunghard, 

I,  Mansfield,  and  0,  Zepko, 

K,  C,  Rogers 
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2,  D.  Finkelstein,  The  Megatron,  NYO-7735>  (19^7)  >  D.  Finkelstein,  The  Megatron, 
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"Behavior  of  Runaway  Electrons 
in  a  Plasma  Betatron" 

L.  A.  Ferrari  and  K.  C.  Rogers 


ABSTRACT 


BEHAVIOR  OF  RDIJAWAY  ELECTRONS  IN  A  PLASMA  BETATRON 


Runaway  electrons  generated  by  induction  acceleration  in  a 
plasma  betatron  with  axially  symmetric  guide  field  (r,z^t)“B^ 

(  >$/m)  <  -f  ij)  Vl'  (n=0,6,'t  =i6/4  sec  ,  «!4.8  cm)  together 

with  auxiliary  azimuthal  guide  field  B^(  ■^)“(k/r)  ,(2  kG  ^  B^(  r)  ^  5kG) 
have  been  studied  for  a  wide  range  of  induction  electric  fields 
E^(2,5  v/cm  i (E)  426v/cm)^  Experiments  have  been  carried  cut  in  argon 
and  krypton  at  ambient  pressures  of  0.5-2  Torr.  The  ratio  of  the 
applied  electric  field  to  the  critical  electrical  field  for  runaway 
was  typically  /■./  500  so  that  most  of  the  electrons  should  run  away. 
Runaway  currents  of  lA  were  observed  immersed  in  a  background 
conduction  current  of  25-100  A.  The  runaway  electrons  st.ike  the 
Vacuum  chamber  walls  at  a  time  t^  after  betatron  acceleration  has 
begun.  The  resulting  x-rays  are  detected  with  a  scintillation  counter. 
The  x-ray  time  t^  has  been  measured  as  a  function  of  and  B^  for 
several  half-cycles  of  the  guide  field.  It  is  found  that  t^  decreases 
with  but  increases  with  B^.  It  is  also  found  that  t^  depends  on 
the  relative  orientation  of  E^  with  respect  to  It^.  This  is  shown  to 
be  due  to  a  small  transverse  magnetic  field  associated  with  the  B^ 
guide  field.  The  small  runaway  current  is  consistent  with  that 
determined  from  the  self  electric  field  generated  by  the  adiabatic 
constriction  of  the  runaway  stream  or  the  Negative-Mass-Instability . 


Tiie  reason  for  the  beam  disruption  mid  subsequent  x-ray  emission  are 
unexplained  but  the  mechanism  of  Field  and  Fried  appears  to  describe 
the  process. 

The  betatron  guide  field  was  generated  in  a  low  inductance 
(50  X  10  single-turn  coil  using  flux-concentrators.  Techniques 
used  to  construct  the  flux-concentrators  euid  measure  the  resulting 
pulsed  magnetic  I'leld  distribution  are  also  described. 


Lawrence  A,  Ferrari,  Author 

Professor  Kenneth  C.  Rogers 
Thesis  Advisor 

Depai-tment  of  Physics 
Stevens  Institute  of  Technology 
Castle  Point  Station 
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INTRODUCTION 

(1) 

In  1956  Budker  suggested  that  the  principle  of 

magnetic  self »f ocualng  of  a  relativistic  stream,  first 

(2) 

discovered  by  Bennett,  might  be  used  to  build  a  new  type 
of  high  energy  particle  accelerator.  This  principle  is 
based  on  the  fact  that  it  takes  a  relatively  small  amount 
of  positive  charge  added  to  an  intense  relativistic 
electron  beam  to  replace  ..he  Coulomb  repulsion  by  strong 
attraction.  The  magnetic  self  focuslnf^  produces  electro¬ 
magnetic  radiation  which  then  damps  the  transverse 
oscillations  of  the  electrons.  This  would  then  cause  the 
beam  to  shrink  to  a  thin  thread  with  very  large  electric 
and  magnetic  fields  on  the  surface. 

Budker  proposed  that  the  self  magnetic  field  of  the 
beam  so  formed  could  then  be  used  as  the  guld  field  for 
the  acceleration  of  ions.  The  electrons  would  travel 
large  distancos  while  the  beam  diameter  was  shrinking, 
therefore  it  was  obvious  that  a  circular  type  of  particle 
accelerator  such  as  a  betatron  would  be  necessary  to 
hold  the  electrons  a,  while  the  stream  was  shrinking  and 
b,  while  the  ions  were  being  accelerated, 

Budker  envisioned  a  circulating  beam  current  of 
17,000  amp  (obtained  when  the  number  of  electrons  per 
length  of  stream  equal  to  the  classical  radius  of  the 
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electron  was  equal  to  one  )  however,  beam  current 

(3) 

caloulatlona  by  Flnkelstein  and  Kaleonnler  and  finally  by 

(A) 

Schmidt  showed  that  the  beam  diamagnetism  would  limit  the 
maximum  beam  current  in  plasma  betatrons  to  several 
hundred  amperes. 

In  spite  of  this  drastic  reduction  In  maximum  beam 

current  from  that  envisioned  by  Budker  the  plasma  betatron 

(5) 

Is  a  device  which  still  warrants  Interest.  Falnborg  In 
hie  review  article  lists  approximately  twenty  different 
Instabilities  that  can  be  generated  by  the  passage  of  a 
charged  particle  stream  through  a  plasma.  The  energy  lost 
by  the  beam  particles  In  exciting  oscillations  can  be 
considerable;  when  the  number  of  particles  per  bunch 
/«  - //  the  energy  loss  can  be  as  high  as  -  /o 
per  bunch.  Since  this  energy  lose  can  be  so  large  It  Is 
obvious  that  these  Interactions  may  play  an  Important  role 
In  plasma  physics  and  hone o  be  worthy  of  study.  In  a  plasma 
betatron  the  electrodes  which  are  present  tn  most  beam- 
plasma  studies  are  absent.  This  coupled  with  the  toroidal 
geometry  enables  a  very  long  path  length  for  the 
circulating  electrons.  In  principle  It  Is  also  possible 
to  accelerate  the  electrons  to  any  velocity  between  their 
Initial  thermal  velocity  and  essentially  the  velocity  of 
light .  The  magnetic  guide  fields  are  sufficiently  simple 
so  that  the  single  particle  motion  can  be  studied  analytically 


and  thuB  may  help  to  understand  the  experimental  beha\ior 


of  the  device, 

(6) 

Plasma  betatrons  with  low  Inductance  b.b  well  as  high 

inductance  can  be  constructed.  High  Inductance  plasma 

(7-10) 

betatrons  consist  of  a  symmetric  arrangement  of  a  group  of 

single  turn  colls  arrap'^od  geometrically  so  that  the 

resulting  magnetic  field  In  the  median  plane  between  them 

satisfies  the  betatron  condition.  A  low  Inductance  plasma 

betatron  Is  c>onatructed  within  a  single  turn  solenoid. 

For  the  latter  case  the  vacuum  magnetic  field,  which  Is 

essentially  uniform  must  be  distorted  (shaped)  In  a 

cyllndricaily  symmetric  fashion  until  a  betatron  field  Is 

obtained.  The  possibility  of  utilizing  the  high  magnetic 

field  that  can  be  generated  In  such  a  coll  for  a  small 

orbit  accelerator  has  been  pointed  out  by  Furth,  Levine 
(11)  (12) 
and  Wanlek  and  also  by  Flnkelstoln, 

(7) 

The  betatron  studied  by  Linhart  and  Reynolds  and 

(8) 

Skarsgard  at  CERN,  has  an  azimuthal  magnetic  field  In 

addition  to  the  betatron  field.  Results  were  discouraging 

since  runaway  currents  of  only  /  amp  were  generated, 

(9) 

Drees  at  Bonn  has  constructed  a  hl^  inductance  plasma 
betatron  that  uses  a  high-frequency  quudrupole  electric 
field  for  additional  stabilization.  The  maximum  runaway 
current  observed  in  this  devise  Is  only  ^  ^  amp>, 
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\10) 

Shepherd  and  Skarsgard  have  studied  a  plasma  betatron 
similar  to  the  CERN  machine.  Again  the  runaway  current  is 
of  the  order  of  a  few  amperes. 

The  low  inductance  plasma  betatron,  the  subject  of 
this  study,  is  considerably  more  difficult  to  construct 
than  the  high  inductance  type.  The  advantage  of  this  type 
of  betatron  is  that  the  low  inductance  (  <*.  h  ) 

will  in  principle  allow  the  field  at  the  orbit  to  reach 
several  kilogauss  in  sec  or  less  and  it  is  for  this 
reason  that  this  type  was  chosen  for  this  study  rather 
than  the  high  Inductance  device.  Although  the  low 
inductance  is  highly  desirable  in  a  plasma  betatron,  the 
construction  of  the  field  coll  and  allied  a.pparatus  poses 
more  experimental  problems  than  are  encountered  in  the 
high  inductance  type.  Among  these  are:  a,  obtaining  a 
high  precision  magnetic  field  shape  in  a  small  device, 
b,  voltage  insulation  because  of  the  combination  of  small 
size  and  high  electric  fields,  c,  production  of  a  hignly 
ionized  low  density  plasma  in  a  small  vacuum  chiraber, 
d,  a  high  degree  of  azimuthal  uniformity  in  the  plasma 
density  is  also  required,  therefore,  the  plasma  must  be 
generated  within  the  torus,  e,  all  measurements  are  more 
difficult  because  of  the  small  size,  f,  the  time  required 
to  build  up  a  steady  curre  layer  in  the  field-coil  takes 
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about  three  half  cycles  of  the  betatron  field,  therefore, 
magnetic  field  measurements  with  A.C,  excitation  cannot 
be  used  (appendix  #1).  Other  techniques  for  making 
precise  measurements  of  the  spatial  distribution  of  a 
pulsed  magnetic  field  had  to  be  developed. 

Construction  of  a  low  inductance  plasma  betatron  at 
Stevens  Institute  of  Technology  was  initiated  in  1957  by 
Prof*  D,  Plnkalstein.  The  present  betatron,  which  is  quite 
different  from  the  device  originally  designed  by  Finkelstoln, 
gradually  evolved  to  its  present  state  because  of  the 
technical  problems  encountered  in  the  earlier  ideas. 

In  this  report  we  describe  the  construction  and 
operating  characteristics  of  a  low  Inductance  plasma 
betatron.  A  cyllndrically  symmetric  magnetic  field 
distribution  that  satisfies  the  betatron  condition  is 
obtained  by  placing  flux  concentrators  at  appropriate 
positions  within  the  interior  of  a  single  turn  solenoid. 

In  addition  to  the  betatron  field  there  is  an  azimuthal 
magnetic  field  that  is  required  for  runaway  production. 

It  is  observed  in  this  experiment  as  well  as  In  the 
other  plasma  betatrons  that  beam  disruption  and  subsequent 
x-ray  emission  begins  before  the  end  of  the  accelf ration 
cycle  has  been  reached.  The  runaway  electrons  hit  the 
walls  of  the  vacuum  chamber  without  benefit  of  the  beam 
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disturbing  techniques  used  in  ordinary  betatrons.  The 
time  at  which  x-ray  emission  starts,  and  the  number  of 
runaway  electrons  have  been  studied  as  functions  of  the 
various  experimental  parameters. 

We  find  that  the  runaway  current  is  limited  to  about 
one  ampere  by  either  one  of  two  mechanisms;  the  self 
electric  field  that  is  generated  by  the  adiabatic  collapsf 
of  the  runaway  stream  or  the  Negative  Mass-Instability, 

The  streaming  velocity  appears  to  be  limited  by  a 
beam  plasma  interaction  that  has  been  treated  theor-  tically 
by  Field  and  Fried  (see  section  ifZ)  > 

The  maximum  electron  energy  at  the  onset  of  x-ray 
emission  i^as  approximately  3^0  Ae  /  ,  Although  some  x-ray 
emission  occurs  at  the  peak  of  the  acceleration  cycle 
indicating  that  the  full  energy  of  the  machine  was  reached 
(  ^  /.5  ma\/.  )  the  Intensity  was  greatest  at  the  beam 
disruption  time.  If  we  choose  the  experimental  parameters 
that  give  the  greatest  x-ray  energy  for  a  given  accelerating 
field  strength,  wo  find  that  E,  ^  lo  v/om  is  required  if 
the  electrons  are  to  have  a  kinetic  energy  of  when 

they  strike  the  walls  of  the  vacuum  chamber.  This  electric 
field  is  about  three  times  greater  than  the  experimentally 
applied  fields.  The  extrapolated  runaway  current  with 
El  =70  v/cra  is  only  5  amperes. 
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In  section  #2  we  present  a  brief  review  of  the  theory 
of  runaway  electrons.  Section  #3  deals  with  the  single 
particle  motion  In  a  betatron  magnetic  field  and  In  a 
combined  betatron  and  azimuthal  magnetic  field.  The 
particle  motion  Just  after  the  betatron  field  Is  turned 
on  is  also  discussed  here. 

In  section  #4  the  generation  of  a  betatron  magnetic 
field  using  flux  concentrators  is  described.  Techniques 
for  making  precise  measurements  of  pulsed  magnetic  fields 
are  also  described  in  this  section. 

Section  //5  deals  with  the  generation  of  the  Bq  field 
and  the  method  used  for  prelonlzation.  The  electron  density, 
electron  temperature,  runaway  current,  x-ray  measurements, 
and  electric  field  measurements  are  described  In  section 
/5«6. 

Finally  the  results  of  the  measurements  are  analyzed 
and  discussed  in  section  //7. 
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II  THEORY  OP  RUNAWAY  ELKCTHONS 


Tha  subject  of  electron  runaway  In  Ionized  cases  has 
u  .  (13-19) 

been  studied  by  many  authors.  In  all  of  the  calculations 
the  electric  field  ,  which  accelerates  an  electron  by 

its  thermal  velocity  in  one  collision  period,  enters  as  a 
parameter.  If  the  applied  electric  field  F/4  then 
the  resulting  current  is  linearly  related  to  the  perturbing 
field.  Tha  stationary  methods  of  solution  of  the  collision 
dominated  Boltzman  equation  ip., 


^  ^  J.e  •  ^  =  /2£\  2.1 

are  assumed  to  yield  correct  answers  provided  the  average 
random  electron  speed  is  much  larger  than  the  electron 
drift  velocity.  In  this  limit  the  electrical  conductivity 
follows  the  well  known  7*  law  (Spltzer  conductivity). 

The  calculations  are  based  on  the  assumption  that  a 
steady-state  velocity  distribution  is  attained  several 


mean-free  collision  times  after  the  electric  field  is 

4.  .  (13-18) 

turned  on.  The  present  calculations  relax  this  last 

restriction  on  the  problem,  l.e.,  a  time  dependent  solution 

of  aq.  (2,1)  is  sought.  In  the  next  paragraphs  we  give 

tha  essential  results  of  all  the  runaway  electron 

calculations  rather  than  discuss  the  work  contained  in 

each  paper  separately. 


Whon  £* ,  the  application  of  the  external 
electric  field  causes  the  electron  distribution  function 
to  be  only  slightly  disturbed  from  its  initial  equilibrium 
situation.  If  the  distribution  f\mction  resembles  a 
Maxwellian,  ie,  if  it  has  a  high  energy  tall,  there  wll] 
always  be  some  electrons  with  velocities  large  enough  to 
runaway  no  matter  how  small  the  applied  field  is.  As  far 
as  the  main  bulk  of  electrons  are  concerned,  the  application 
of  the  electric  field  causes  heating,  thereby  enabling 
more  electrons  to  move  into  the  runaway  region,  'Tie 
velocity  of  any  given  electron  gradually  increases  with 
time.  Thus  all  the  electrons  eventually  runaway. 

When  I  ^Me  electric  field  term  in  the 

Boltzman  equation  dominates,  all  the  electrons  runaway 
immediately,  and  the  runaway  current  increases  linearly 
with  time  according  to 


/w 


(2,2) 


(19) 

Recently  Field  and  Pried  have  attacked  the  problem 
of  electron  runaway  in  a  way  that  is  very  different  from 
the  previous  calculations.  Instead  of  the  conventional 
Fokker-Planck  collision  term  for  the  right-hand  side  of 
2,1  they  assume  that  arises  from  the  generation 


(-] 

of  ion  acoustic  waves.  The  calculation  is  done  in  the 
limit.  Thus  initially 


is  of  the  same 


order  of  magnitude  as  the  conventional  collision  term  so 


that  the  current  density  -j  increases  linearly  with 
tim;’».  As  soon  as  it  passes  a  small  threshold 

>lx 


y  ^  /rtc.  AT 


/  \ 
[  j 


however,  the  (9^]  term  begins  to  rise  exponentially 
and  eventually  becomes  comparable  with  the  driving  term. 
The  current  density  reachss  a  maximum  and  then  decreases 
rapidly. In  time. 

To  suinmarize,  for  efficient  runaway  production  the 
induction  electric  field  £r^  should  be  much  larger  than 
f  where 

•  /4  \ 

F  it  t.H  /O  W 

‘  Tc 


In  this  experiment  the  ratio  of  is  very  large, 

typically  500  so  that  all  the  electrons  should  runaway 
Immediately.  Furthermore  the  ?oulomb  collision  time  is 
(0  7^  =  10  cl/,  ),  This  is 

more  than  twice  as  long  as  the  acceleration  cycle,  hence 
the  runaways  will  move  through  the  backround  plasma  in  a 
oolllslonless  way. 

Further  discussion  of  the  work  by  Field  and  Fried 
will  be  left  to  section  #7« 


Ill  EQUATIONS  OF  MOTION  FOR  SINGLE  PARTICLES 


In  this  seotlon  we  derive  the  equations  of  motion 
for  single  particles  In  an  axially  symmetric  betatron 
magnetic  field  and  In  a  betatron  field  with  an  additional 
azimuthal  guide  field. 

3.1  The  non-relatlvlstlc  L^grangian  for  a  charged 
particle  with  mass  m,  and  charge  q  In  a  magnetic  field 
with  vector  potential  is  given  by 

/  -  -f 

In  cylindrical  co-ordinates  )  the  Lagranglan  takos 

the  form 

/  =  lym  -t  3,2 

The  equations  of  motion  for  the  particle  are  given  by 
Lagrange’s  equations : 

^/iL)  -  iA,  =0 
(  ifj  3/;  3-3 

VUiere  the  ^’5  are  the  generalized  co-ordinates.  Thus 
the  Jt  ,  S  e-i^d  ^  equations  of  motion  are 


f » 


yiyw  III 


12 


0  V  I  20  t  29  I  ^ 


If  we  perform  r.he  differentiation  with  respect  to  time 
and  make  use  of  B  -  V ,  equatlon53.4  become 

/m  'ji  *  ,tft  Ji  9^  -  j  ^  ^  ^Jt  &  B^ 

^  A  A  e  Bji 

,.  '  3.5 

,m  A  9  -i-  ^  /rn  A  6  — 

iwe  have  neglected  the  time  dlpencfence^  of  A,} 


,m  A  9  -i-  ^  /rnA  6  —  "/•^Sa 
iwe  have  neglected  the  time  dlpencfence^  of  %) 

The  first  field  configuration  we  will  study  will  be 
Just  the  betatron  field.  For  this  situation  6^=0  so 
that  9q.'n  (3.5)  reduce  to 


Ar>  =  yViA0^+aJt0B^ 


?  "  -r 


Bx 


*  *  •  • 

yw  A  9  -t  ^  /n\  k  Q  =  -  ^  ^  ^ 

Before  we  analyze  the  equations  of  motion  any  further  It 
will  be  helpful  If  we  flret  derive  the  neoeaeary  conditions 
for  an  equlllbrluia  orbit.  Only  one  oompon«nt  of  the  vector 
potential,  ,  la  necessary  to  generate  an  axially 
symmetric  betatron  field.  Thus  the  Lagrangian  (3,2) 


reduces  to 
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We  will  assume  that  is  cylindrical ly  symmetric 

and  we  will  also  assume  that  is  time  independent. 

Since  fig  is  symmetric,  the  $  -  component  of  the 
canonical  momentum  Pm  z  will  be  a  constant  of  the 

®  de 

motion.  We  choose  =0.  Thus 


p.  s  =  ym  ^  -f-  9 

0  f 


Afia  ~  0 


Solving  for  (Ad  )  in  eq,  3*6  and  substituting  into 
©q*(37),  the  Lagranglan  becomes 


Thus  vrlth  this  Vector  Potential,  the  ^  ^  ^  equations  of 


miOtlon  are: 


/«  3  ±-  f)e  --0 

/«  3^ 


3,10 


Therefore  at  the  equllibrluni  orbit  both  and 

•aiiS  must  vanish.  The  familiar  A  ^  /  betatron 
condition  is  obtained  by  calculating  the  magnetic  flux 
^  that  links  the  equilibrium  orbit  with  radius  X^Ao  • 


lA 


J=  /I  6-i’l  =  ^6>^V 

where  ^^>l8  the  average  value  of  the  magnetic  field 
within  the  orbit.  We  also  have  that 

if  J 

fjf*)  <=  jjTl  3.12 


where  A^(i^)  la  the  value  of  the  vector  potential  at  the 
equilibrium  orbit.  From  B^^XA 


j. 

?/^6'  f  fle 
'nr  ^  -jr 

3.13 

since  1  f 

;; — 1  c  0  ,  we  have 

*  ^0  M 

J>o 

3.IA 

thus 

=  jTTJ, 

finally 

^3  =  i 

3.15 
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Now  we  return  to  the  equations  of  motion  'to 

calculate  the  particle  orbits.  The  •'z”  component  of  the 

magnetic  flei§  neighborhood  of  the  equilibrium  orbit  is 

(20) 

assumed  to  be  of  the  form 


/n. 


3.16 


The  radial  field  component  associated  with  this  is 

obtained  from  B  -  ^  ie,, 

Now  we  linearize  the  equations  of  motion  ie.,  we  let 

whe  re  I  ^  I  ^  ^  jind  ^  ^  ^  ^ 

To  first  order,  we  find  that 


UJo 


K 

2  ^ 

« *  * 

^  q/  +  14)^ yn  /f 


3.18 


finally,  the  radial  and  vertical  equations  of  motion 
become 


'i 


4  cJ,  /n  /a  =  d 


1 


3.19 


>1 


U; 
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The  solutions  of  the  radial  equation  will  be  harmonic 
providea  /y}  /  I  •  Similarly  the  solutions  of  the  *’  ” 

equation  will  be  harmonic  if  /?»  >  /  •  Hence  we  have 
arrived  at  the  familiar  "weak  focusing  condition",  le.,  in 
a  weak  focusing  field  (a  betatron  field),  the  field  index 
"n"  at  the  orbit  must  lie  between  o  and  /  • 

3*2  At  this  point  we  will  digress  from  the  main  purpose 
of  this  section  to  calculate  the  effect  that  azimuthal 
inhomogenaities  in  the  guide  field  have  on  the  betatron 
orbits.  This  is  a  very  important  consideration  in  the 
design  of  the  betatron  field  because  as  we  shall  show, 
azimuthal  inhomogendtles  increase  the  amplitude  of  both 
the  radial  and  vertical  betatron  oscillations  thei^oby 
placing  restrictions  on  the  vacuum  chamber  cross-section. 

We  assume  that  the  z-component  of  the  magnetic  field 

(21) 

Is  given  by 


3.20 

Ard  that 


,.<1^  P  r 


3.21 


The  equations  of  motion  for  the  radial  and  vertical 


directions  are  now 
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/jt  +  a;/  2  ‘Cn 

/ffft 


3.22 

where  ^  =  cJ^  i: 

The  complete  solution  to  the  radial  equatior.  Is 


4  =  ^  +  B  CAhoY^iOtii-^) 


^2 
/7Tt  =  l 


J>e  /  ytn  uf a -i  -h  o!^  ) 


3.23 

The  c one t ant 8 'a' and *b"  are  determined  by  initial  conditlona. 
The  solution  to  the 'z"*  equation  is  obtained  from  eq.  (23) 
by  simply  replacing  (  )  by  (>»f  )  and 

1^)  <(V  terms,  1,0., 


w 


■r-  iiB  H  I  ^ 
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<:  r  b  c^(io,ynt) 


•o 


Thus  we  see  that  the  harmonloa  that  will  give  the  largest 
displacement  from  the  equilibrium  orbit  are  those  with 
low  'V  values,  ,  The  higher  m-values  will 

not  play  an  important  role  because  of  the  '/a 

/  /itt 

dependence , 

Equations  3.23  will  be  used  later  on  to  calculate 
the  particle  displacements  in  the  field  coll. 

3.3  We  will  now  study  the  solutions  of  the  equations  of 
motion  (3.5)  when  in  addition  to  t:.e  betatron  magnetic 
field,  there  is  an  azimuthal  guide  field  ,  por  this 

case  eq's  3.5  retain  their  full  form  with  0^  ^  the  same 

as  before.  The  linearized  equations  of  motion  about  the 
equilibrium  orbit  \  are 


19 


^  -t  4  /ft  ^ 


-  OtJ^  ^ 


(u  = 


3.2)4 


where  ^^4  =  ^  >  and  Bg  la  the  magnitude  of  the 

azimuthal  field  &t  A  =  Jt,  ,  The  solutions  of  eq  (3.2i|) 


are  assumed  to  be  of  the  form* 


3.25 


substituting  (3.25)  into  (3.21+)  we  get  the  pair  of  equations 
-  (l-ztC)  'f.  -  u)  uJc  -  o 


I  A 

-  oJ  i-  oJt  /n 


jl  a)  ^  ^ 


3.26 


The  condition  for  a  non-trlvial  solution  Is  that  the 
determinant  of  the  cc-t3ff iclants  vanish,  le., 


vOf  ( I  -yfi )  ~  ~  A  u)u)c 

X  (j)  u)^.  /ft  ' 


~  o 


3.27 
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Thus  ’V'  Is  detei-mlned  by  the  solution  of 


3.23 


Finally 


*/k 


3.29 


Thus  in  general  the  electron  motion  is  harmonic  with 
frequencies  that  are  complicated  functions  of 
and  UJ^  ,  However,  when  *Jc  >>  the  roots  of  3.29 


become 


3.30 


3.31 


Prom  (3.26)  we  find  that 


/w  I 

/ 


when  U)  ^  a),  ^  a. 


~  X  !  with  u)  zi  m) 2 


3.32 


3.33 
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The  motion  lo,  therefore,  compoeed  cf  two 
elementary  motions  (Flg.#I),  one  being  nearly  a  cyclotron 
motion  with  anecuiar  frequency  a>j^  'it  t^c.  the  other 

a  slow  elliptic  motion  whoa©  angular  frequency  is 


3*^  In  the  proceeding  section  we  have  studied  the  motion 
of  charged  particles  under  the  combined  influence  of 
axially  symmetric  betatron  and  azimuthal  magnetic  fields 
(  and  ).  The  motion  was  studied  for  the  statio 
case,  l.e.,  both  and  remained  unchanged  during  the 
times  of  interest.  For  the  case  when 

O 

found  that  the  motion  in  the  (x,z)  plane  was  composed  of 
two  elementary  notions  (Fig.#!),  a  nearly  cyclotron  motion 
with  <4^  04.  and,  a  alow  elliptic  motion  with 

»  This  solution  is  correct  as  long 
as  both  )'»' are  sufficiently  large.  However, 

this  is  not  always  the  case  for  a  plasma  betatron.  The 
time  scale  of  the  present  experiment  la  as  follows? 
la  approximately  trapezoidal  In  shape  and  lasts  for  about 


100  yti  i6c. ,  on  the  other  hand  rings  with  a  period  of 
16/1  <ec  .  Therefore  during  an  acceleration  cycle  {a 
qimrter  period),  B^  ( )  rises  from  zero  to  its  maxlnnnn 
value,  while 


22 


Be  is  essentially  unchanged.  The  equations  of  motion 
discussed  in  the  previous  section  are  not  valid  when 
=0,  When  =0  eq.' s  (3*2)4. )reduce  to 

yjt  j 

‘A  « 

^  3.35 

Eq,*  a (3. 35) are  the  equations  of  motion  for  a  charged 
particle  in  a  uniform  magnetic  field  which  Is  not  the  case 
for  8e  since  U  jji  • 

Since  eq.*  8(3.24)are  valid  when  B^yo  ’aot  valid 
when  B«’  “0  we  now  seek  to  determine  a  time  tar  such 

that  the  equations  of  motion  are  valid  for  all 
and  further  more  we  must  also  investigate  the  electron 
motion  for  -6  < 

The  solution  to  (3.24)  with  a  time  varying  is 
very  difficult  to  do  analytically,  however,  we  will  make 
an  estimate  the  range  of  validity  In  the  following  way. 

The  terras  and  Ir.  the  equations 

of  motion  are  measures  of  the  focusing  force  produced  by 
the  betatron  field.  If  the  ^  field  was  Independent 
of  "r”  the  equations  of  motion  would  be  valid  for  all 
time  including  t=c.  Therefore  tne  eqxiations  of  motion 
should  be  valid  for  all  times  that  the  magnetic  focusing 


ELECTRON  MOTION  IN  THE  (X.Z)  PLANE 


I 


FIG.#  I 
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terms  are  larger  than  the  higher  order  terms  In  the 

expansion  of  Bg  about  the  equilibrium  orbit.  If  the 

first  two  terms  of  the  expansion  of  3^  are  retained  the 

equations  of  motion  become, for  B^  yO  >  (consider  only  the 

d 

radial  equation). 


+  uft  (i~yn)  yf-  -  ^  (f  ~  '  /^t) 


3.36 


then  the  solution  (3.31)  and  (3.3if)  should  be  valid  when 


(jol  d-^)  >?  0O<.  ijji^ 


Since  i^o  = 


3.37 

’  Xkt  oq,  (3.37)  reduces  to 

/m  * 


t  » 


-I 


Be  '% 
3^ 


3.38 


Typical  experimental  parameters  are:  =5  k  gauss, 

/y)  =  flj  =1  k  gauss  and  il  39x10^  rad/sec.  If  we 

assume  ^  is  the  velocity  of  a  1C  volt  electron,  then 

tv-  is 


££C. 


W«  must  now  Investigate  the  motion  for  O  ^  ^ 

to  make  sure  that  all  the  electrons  will  not  drift  out  of 

the  acceleration  region  during  this  time.  For  times  less 

than  the  electrons  will  drift  in  the  field  with 

(23) 

their  guiding  center  velocity  which  is  given  by 


whore  =  JL  t  is  the  particle  velocity  parallel 

/f*i 

to  Bfi  ,  Since  increases  more  rapidly  than  either 

yjL  or  ^  we  will  neglect  the  first  terra  in  eq,  (39), 
then  in  a  pure  3^  field  the  drift  is  in  the  "z”  direction 
and  is  equal  to 


/ir»  = 


/m  u/,t 


t"' 


-  - 


3.1^0 

The  distance  that  an  electron  travels  during  the  time 


Is 


3  =  f  J  t  =  X  (  i'  H 

'o  ^ 


3.1+1 


3 


-tJ 


3.1*2 


25 


-8 

With  fs-  30  v/cm,  =  5  ora,  =  lx.10  sec.  and 
=5  k  gauss,  wo  find  that 

Since  the  inside  diameter  of  the  vacuum  chamber  is 
X  /  r  CM.  we  can  expect  that  the  magnetic  focusing 
forces  will  be  large  enough  to  stop  the  "z"  drift  of  most 
of  the  electrons  before  they  hit  the  vacuum  chamber  walls. 
This  transverse  drift  during  the  time  o  /.  t  /L  will 
cause  the  loss  of  some  particles  to  the  walls  but  certainly 
not  all  of  them.  Therefore,  the  fact  that  the  equations 
of  motion  are  not  valid  for  the  entire  part  of  tiie 

acceleration  cycle  does  not  seem  to  lead  to  any  dlaasterous 


effects 


IV  THE  M.\GNh:TIC  FIELD 


General  considerations: 

Probably  the  major  difficulty  encountered  in  the 
design  of  a  betatron  particle  accelerator  is  that  of  the 
magnetic  guide  field.  If  the  magnetic  guide  field  is 
"right”,  then  properly  injected  particles  will  usually  be 
accelerated  to  the  full  energy  of  the  machine,  A  plasma 
betatron,  in  addition  to  being  subject  to  all  the  difficultie 
that  are  encountered  in  electron  betatrons,  will  be  at  the 
mercy  of  any  collective  effects  that  may  develop  during 
the  acceleration  cycle.  We  have  theref':-''e  paid  a  great 
deal  of  attention  to  the  betatron  guide  field  in  this 
experiment,  with  the  hopes  that  we  might  be  able  to 
separata  single  particle  effects  from  collective  effects. 

The  degree  of  precision  with  which  the  guide  field 
must  be  measured  can  be  estimated  from  the  vector  potential 
well  depth.  The  vector  potential  well  depth  X  is 
defined  as: 

f  =  fig  (■»,+(*)  - 
fig  o. ) 


4.1 

whare  la  the  equilibrium  orbit  radlua  and  la 

a  small  displacement  from  , 


if  we  assume  that  the  ''n" -value  is  unifom  in  the  or'bit 
region  then  it  is  easy  to  show  that  the  well  depth  is  given 
by  (to  second  order) 


h.2 


With  n=0.5  and  =2",  the  well  depth  at 

=0.  2'*  (determined  by  the  vacuum  chamber  size)  is 

X= 

Since  the  well  depth  in  a  betatron  turns  out  to  be  very 
small  an  accuracy  of  t  o.t  %  in  all  magnetic  field 
measurements  was  decided  upon,  A  further  complication 
arises  because  of  the  time  that  it  takes  to  build  up  a 
steady  current  layer  in  the  coil.  A  simple  calculation 
of  the  magnetic  field  penetration  time  into  a  conductor 
for  the  frequencies  we  are  interested  in  (  ^  65  k  c ) 

shows  that  at  least  one  complete  cycle  of  the  betatron 
field  is  required  to  build  up  a  steady  state  ''skin  depth*’ 
in  the  field  coll,  (see  appendix  #  1).  This  means  that 
the  field  distribution  in  the  coll  on  the  second  half 
cycle  will  probably  be  different  from  that  on  the  first. 

At  tne  start  of  the  third  half-cycle  (fifth  quarter  period) 
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the  magnetic  field  distribution  will  be  close  to  Its  steady 
state  A.C,  value, 

Tlierefore  If  the  betatron  field  Is  determined  with 
continuous  A,  C.  excitation.  It  Is  quite  likely  that  the 
field  distribution  on  the  first  half-cycle  may  not  satisfy 
the  betatron  condition.  We  have  therefore  m''de  all 
magnetic  field  measurements  under  pulsed  operations.  In 
fact  all  of  the  magnetic  measurements  were  made  under  some 
of  the  conditions  that  were  used  for  experiments  with 
plasma. 


4.2  Field  Coil  Design 

The  design  of  the  plasma  betatron  field-coil  has 
been  determined  largely  by  considerations  of  magnetic 
field  shape  and  accessabllity  rather  than  mechanical 
strength.  Fortunately,  the  final  design  is  sufficiently 
simple  so  that  both  mechanical  and  magnetic  requirements 
are  easily  satisfied. 

The  coil  shown  In  piga j^2&3  1b  machined  from  two 

blanks  of  zirconium-copper  {99 Cu ,0,1$%  )supplied 

by  the  Philadelphia  Bronze  and  Brass  Company,  The  blanks, 

originally  work  hardened,  were  recovered  by  heat  treating 

o 

for  four  hours  at  600c  F  and  then  oven  cooled  for  24  hours 
to  prevent  distortion  during  machining  and  aging. 

Aligr-  u  the  four  "  c's"  in  the  vertical  direction 
is  by  precision  ground  alumina  spacers  resting  in  Jig  bored 
wells.  Proper  separation  between  the  two  halves  of  the 
coil  is  maintained  by  inserting  insulating  spacers  into  the 
gaps  at  the  current  feed  points. 

The  physical  dimensions  of  '  coll,  Inside  diameter 
4»875"»  separation  between  halves  0,750",  overall  height 
0.750",  were  arrived  at  after  magnetic  measurements  in  a 
mock-up  coil  indicated  that,  with  flux-concentrators,  the 
betatron  condition  might  be  satisfied  in  this  type  of 
arrangement , 
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4»3  The  Concentrators 

Flux  concentrators  r  ke  use  of  the  fact  that 
good  conductors  tend  to  exclude  time  varying  electromagnetic 
fields  from  their  interior,  .The  physical  mechanism  for 
the  exclusion  of  the  flux  is  quite  simple.  A  time  changing 
magnetic  field  has  associated  with  it  an  induction  electric 
field.  The  electric  field  in  the  conductor  generates  eddy 
curremcs  that  are  Just  large  enougii  to  cancel  the  applied 
magnetic  field. 

Each  flux  concentrator  is  composed  of  16  copper 
discs  1/16"  thick  with  an  outside  diameter  of  4,8l3" 
t  <?.001"  and  an  Inside  diameter  of  1,5000"  i- ,  000^" , 

A  .020"  slot  in  each  discs  permits  magnetic  flux  to  enter 
the  central  hole.  (Fig. #5)  'The  discs  are  Insulated  from 
each  other  by  a  ,010"  thick  mylar  sheet.  Each  of  the  16 
copper  pieces  and  mylar  insulators  is  coated  with  Emerson 
&  Cummings  1264  epoxy  resin  and  then  mounted  on  a  mandrel 
which  maintains  alligmnent  between  the  discs.  The  mandrel 
is  designed  so  that  no  disc  is  out  of  alignment  from  any 
other  by  more  than  0.001".  The  entire  assembly  is  then 
clamped  and  allowed  to  harden.  Staggering  the  gaps  in 
alternate  layers  reduces  the  azimuthal  field  inhomogenS/ties 
that  are  caused  by  the  slots  in  the  discs. 


FIG.  9^^  2 


FIELD  COIL 
(side) 
FIG. ^3 
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ij.»4  Difference  Amplifier  Technique  For  Measuring  The 
Magnetic  Field 

All  of  the  magnetic  field  measurements  were  made 
with  two  magnetic  pick-up  loops.  One  of  the  probes,  fixed 
in  space,  was  used  to  monitor  the  current  through  the 
Aield-coll.  The  second  probe  {  in  some  cases  a  moving 
probe)  was  used  to  map  the  magnetic  field.  Each  signal 
was  integrated  In  an  R-G  network  and  fed  into  a  Tektronix 
Type  ^  differential  preamplifier  mounted  in  a  Tektronix 
Type  555  dual  beam  catbode-ray  oscilloscope.  The  difference 
between  each  pick-up  loop  signal  and  the  variable  D.C, 
bias  voltage  built  into  the  preamp,  was  recorded  on 
Polaroid  film, (In  principle  it  is  possible  to  see  a  change 
of  S  /my.  in  a,  100  volt  signal). 

The  current  through  the  field-coil  varied  from 
5-10^  between  successive  capacitor  bank  discharges  because 
of  th‘/  uncertainty  in  the  charging  voltage.  Under  these 
conditions  one  must  keep  discharging  the  capacitor  bank  at 
each  probe  position  until  the  signals  from  both  pick-up 
coils  were  "on  screen".  Installation  of  the  capacitor 
bank  voltage  measuring  circuit  shown  in  Fig, 4  removed 
this  difficulty.  With  this  circuit  it  was  possible  to 
discharge  the  capacitor  bank  so  that  signals  from  both 
the  fixed  probe  and  moving  probe  were  "on  screen"  every 
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discharge.  The  circuit  is  balanced  so  that  there  is  zero 
current  in  the  middle  leg  wtien  a  calibrated  high  voltage 
source  is  connected  across  the  input  (in  our  case  a  John 
PluJte  .Precision  H.V,  Supply,  with  a  0,1^  accuracy).  The 
bridge  will  tnen  indicate  zero  current  when  the  capacitor 
bank  voltage  reaches  the  voltage  used  for  balance.  With 
an  RCA  ultra  sensitive  microammeter  in  the  middle  leg  it 
la  possible  to  see  a  difference  cf  one  part  in  i|.000  in 
th  charging  voltage.  Normally  the  capacitor  bank  would 
be  allowed  to  overcharge  approximately  10  to  7^  volts. 

The  charge  was  then  allowed  to  leak  off.  Switch  Swap 
opened,  disconnecting  tne  microammeter,  and  simultaneously 
generating  a  pulse  which  initiated  the  bank  discharge. 

This  entire  process  was  repeated  at  least  three 
times  at  each  probe  position.  All  the  data  ( li.  a  given 
set  of  circumstances)  were  then  normalized  to  the  same 
value  of  field-coil  current  (provided  the  location  of  the 
fixed  coll  was  not  changed).  The  spread  in  the  data 
obtained  this  way  was  typically  +  ,0$%, 


TO  CAPACITOR  BANK 
100  m 


\ 
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4.5  Shaping  the  Field 


The  radial  distribution  of  the  2-component  of 
tiiO  magnetic  field  in  the  median  plane  of  the  field  coll 
la  determined  by  the  following  field-coil  and  flux  concen¬ 
trators  parameters 

a=  Tlie  diameter  of  the  central  hole  in  the  flux 
concentrators 

b-  The  separation  between  the  flux  concentrators 
and  the  interior  wall  of  the  field  coll 
c-  The  height  of  the  flux  concentrators  themaelves 
d=  The  vertical  gap  between  the  two  halves  of  the 
f ^eld-coil 

e=  The  axial  spacing  between  the  two  flux  concentrators. 
Suitable  choices  for  dimensions  "a"  through  "d"  were 
determined  by  magnetic  field  measurements  in  a  model 
field-coil. 

In  order  to  discuss  the  roles  the  different 


parameters  play  in  determining  the  field  dl strlbutlonp  it 
is  convenient  to  introduce  the  parameter  "C',  0  is  defined 

as  the  ratio  of  the  magnetic  field  strength  a^o)  on 

the  axis  of  the  field  coll  in  the  median  plane,  to  the 
field  at  the  equilibrium  orbit  radius.  With  the  present 
geometry  the  vector  potential  has  a  minimum,  with  the 
"n"  value  betw'-  n  /.jro  and  one  at  Mq  =1,9  inches  for 


values  of  G  betvieen  4  and  5. 

P’leld  distributions  with  G>5  generate  a  vector 
potential  that  decreases  nonotonicaly  in  the  orbit  region* 
For  a  given  set  of  flux  concentrators  "G"  decreases  with 
increasing  ’’a".  In  general  the  larger  "a"  is  with  respect 
to  the  inside  diameter  of  the  field  coil  the  srnallejG  will 
be.  The  dimensions  "b",  "d"  and  "e"  control  the  fine 
details  of  the  magnetic  field  in  the  orbit  region,  ie.,  the 
"n"-value.  Once  a  set  of  flux  concentrators  has  been 
constructed  the  only  parameter  that  can  be  changed  easily 
is  "b"*  This  is  accomplished  by  adding  an  additional  ring 
to  the  innermost  flux-concentrator  disc  (see  flg.£). 
Obviously  this  procedure  only  works  when  "b"  is  to  be 
made  smaller.  An  increase  in  the  value  of  "b"  or  a 
decrease  in  '*d"  will  cause  the  "n''-value  to  become  negative 
for  values  of  (r)  greater  than  the  equilibrium  orbit 
radius. 

In  the  few  proceeding  paragraphs  we  have  described 
in  general  terms  the  roles  the  various  parameters  play 
in  determining  the  radial  profile  of  the  magnetic  field  in 
the  coll.  A  more  quantitative  description  would  require 
an  extensive  set  of  magnetic  measurements  covering  many 
different  values  of  all  the  parameters,  most  c?  which  did 
not  lead  to  a  suitable  betatron  field. 
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Since  we  are  only  interested  in  those  experimental 
situations  that  produce  an  acceptable  betatron  field  shape, 
a  thorough  set  of  magnetic  measurements  was  not  made  until 
it  was  reasonably  clear  that  a  given  set  of  parameters 
would  produce  the  desired  result.  Table  ^  lists  the 
values  of  the  different  parameters  which  produced  a 
betatron  field  in  the  median  plane  of  the  field  coll. 


Table  #  1 


a-  1.500" 

b=  0.031" 

c=  1,250" (approximate  ) 

d=  0.750" 
e=:  1,750" 


I4..6  Radial  Distribution 


The  radial  distribution  of  the  "z"  component  of 
the  magnetic  field  A)  in  the  field-coil  was  measured 

with  a  small  pick-up  coil  (1/8“  O.D,xl/l6''  high)  that 
was  constrained  to  move  along  a  diameter.  The  pick-up 
coll  was  mounted  on  one  end  of  a  length  of  bakelite  rod 
that  had  a  square  cross-section  (0,250"  x  0,250"),  A 
probe  holder  that  could  be  mounted  in  the  gap  between  the 
two  halves  of  the  field  coil  had  ll}.  radial  grooves 
machined  into  it.  The  grooves  were  machined  so  that  the 
probe  rod  could  move  no  more  than  t  .002"  in  either  of 
the  vertical  or  lateral  directions,  Aligni.ient  of  the  probe 
holder  with  the  geometric  axis  of  the  field  coil  was 
accomplished  by  drilling  holes  for  the  alumina  ceramic 
spacers  that  are  used  for  vertical  alignment  of  the  field 
coil.  The  vertical  positioning  of  the  probe  holder  was 
done  by  adding  shims  until  the  desired  "z"  position  was 
reached  (most  of  the  measurements  wore  done  at  z=o",  the 
geometric  mld-pleane  of  the  coll).  Radial  motion  and 
position  was  obtained  by  connecting  the  probe  rod  to  a 
precision  screw  thread  (10  threads  /cm).  The  screw  thread 
was  graduated  in  units  of  (1/200  mm)  so  that  precise 
radial  positioning  could  be  maintained. 


4,7  The  Vector  Potential 


The  vector  potential  in  the  fleld-coll 

was  obtained  two  different  ways.  First,  by  numerically 
Integrating  the  expression 


X  /  -A' 


4.3 


where  {  A' )  is  the  radial  distribution  of  the  magnetic 

field,  and  secondly,  by  measuring  the  voltage  induced  in 
a  series  of  single  turn  pick-up  coils  concentric  with  the 
geometric  axis  of  the  field  coil.  The  first  method, 
because  it  is  more  time  consuming  was  only  used  to  check 
the  second,  and  then  only  after  the  single  turn  coils  had 
indicated  that  a  vector  potential  minimum  existed. 

The  desired  accuracy  of  our  magnetic  measuremente, 
(1  ,1%)  requires  that  tbcradlus  of  each  of  the  pick-up  colls 
be  known  to  at  least  t  ,001",  'fhe  single  turn  colls  were 
made  of  0,007”  copper  wire  which  was  glued  Into  0,007" 
grooves  machined  in  a  piece  of  plexiglass.  The  diameter 
ef  each  groove  was  measured  with  an  optical  comparator 
and  found  to  be  wl th into, 000^”  of  the  nominal  value,  A 
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length  of  micro-dot  cable  was  connected  to  each  loop  via. 
a  1/16"  X  1/16"  radial  slot  machined  In  the  plexiglass 
under  the  wires.  This  type  of  connection  allows  the 
perturbation  at  the  connection  to  be  minimized, 

A  second  source  of  error  arises  If  there  la  a  lateral 
displacement  of  the  set  of  pick-up  colls  with  respect  to 
the  axis  of  the  machine.  One  can  show  (appendix  #2)  that 
the  fractional  error  In  the  signal  from  a  pick-up  coll 
with  radius  "r"  due  to  a  lateral  displacement  "a"  is 
given  by 

-Z  (l-^) 

4.4 

If  we  assume  that  n=^,  and  take  r=2"  then  with  a=,010" 


The  lateral  displacement  of  0,010"  assumed  above  would 
only  result  from  a  severe  mis-alignment  In  the  fleld-coll. 
The  measured  displacement  Is  ,006",  Since  the  error  due 
to  small  lateral  displacements  Is  extremely  small,  It  has 
been  neglected  In  all  of  our  treatment  of  the  experimental 
data. 


4*8  Azimuthal  Field  Measurements 


The  azimuthal  variation  of  the  magnetic  guide 
field  0^  i  Jl  )  was  measured  along  the  equlllbrlura  orbit 
radius.  The  probe  O.D,  was  1/8'*  and  its  height  was  l/l6'** 
The  probe  was  mounted  in  a  plexiglass  disc  that  was  free 
to  rotate  in  a  shallow  well  machined  into  another  piece 
of  plexiglass.  Alignment  of  the  probe  holder  in  the 
radial  and  vertical  directions  was  done  exactly  the  same 
as  it  was  for  the  radial  probe. 

The  shallow  wall  holder  was  also  used  to  locate 
the  single  turn  pick-up  coils  in  the  field  coll. 
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l4.«9  Results  of  the  Magnetic  Field  Measurements 

The  radial  distribution  of  the  betatron  guide 
field  3^  on  the  first  half-cycle  is  shown  in 

Flg,#^,  The  '’n”-value,  measured  at  twelve  different 
azimuths  is  displayed  in  Flg*jS^2.»  vector  potential 

{Ji  )  that  corresponds  to  the  guide  field  of  Fig. #6 
Is  shown  in  Fig. #8.  The  minimum  in  ,  hence  the 

equilibrium  orbit,  occurs  at  r=1.9''.  The  vector  potential 
at  the  peak  of  the  second  half -cycle  is  alb‘>  shown  in 
Flg.j^S,  As  we  can  see,  fig  on  the  second  ha3.f -cycle 
is  not  the  same  as  it  was  on  the  first.  All  of  the  magnetic 
measurements  that  we  have  described  so  far  have  been  made 
at  the  peaks  of  the  field  on  the  respective  half-cycles. 

The  fact  that  the  vector  potential  at  the  peak  of  the 
second  half -cycle  is  different  from  that  on  the  first 
raises  the  question  of  whether  or  not  the  field  shape  is 
correct  during  the  entire  first  quarter  period.  With  the 
differential  amplifier  techniques  previously  discussed  it 
is  possible  to  measure  only  peaks  (in  time)  with  high 
precision  hence  a  different  procedure  had  to  be  used. 

As  a  first  step  in  this  direction  an  x-y  densitometer  was 
constructed  eo  that  the  amplitude  of  the  signal  from  a 
magnetic  oick-up  loop  (recorded  on  Polaroid  film)  could 
be  measured  as  a  function  of  time. 


FIELD  STRENGTH  (ARBITRARY  UNITS) 
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RADIAL  DISTRIBUTION  OF  BETATRON  FIELD 
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The  densitometer  consists  of  a  Derbyshire  lathe  cross-head 
(for  the  j-j  motion),  modified  to  hold  the  Polaroid  film, 
and  a  microscope.  The  cross-head  and  microscope  were 
rigidly  attached  to  a  one  inch  thick  steel  plate. 

Operation  as  a  densitometer  is  carried  out  by  connecting 
a  93IA  photomultiplier  tube  to  the  microscope  eyepiece. 

The  field  of  view  could  be  changed  by  placing  appropriate 
apertures  in  the  plane  of  the  cross-hairs.  Aperture  sizes 
ranged  from  ,001"  to, 005"  dia.  With  a  sufficiently  small 
aperture  the  accuracy  of  the  device  is  determined  by  the 
cross-head  and  the  trace  width  on  the  film.  With  the 
Derbyshire  cross-head  co-ordinates  could  be  measured  to 
within  f  .P005",  Since  this  is  a  fixed  number  it  is 
obvious  that  the  error  associated  with  an  amplitude 
measurement  depends  on  the  amplitude,  ip.,  the  larger  the 
amplitude  the  smaller  the  percentage  error,  A  second 
complication  arises  because  of  electrical  noise  on  the 
signal  which  las' 3  for  approximately  of  the  first 
quarter  cycle.  This  noise,  and  width  of  the  trace,  though 
very  small  for  most  purposes,  was  still  large  enoughto 
make  the  error  at  the  start  of  the  signal  larger  than  the 
0,1^  we  discussed  previously.  Nevertheless,  we  still 
measured  Rg  {A  )  at  different  Instants  of  time  during 
the  first  quarter  period. 
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The  accui>acy  of  these  measurementa  was  only  0,75-1,0^ 
but  served  to  rhow  that  nothing  drastic  was  happening  to 
^0  during  this  time. 

Because  of  the  extreme  difficulties  Involved  in 
measuring  the  magnetic  field  at  times  other  than  at  the 
peaks,  no  further  attempts  were  made  to  improve  the  accuracy. 
The  azimuthal  variation  of  the  betatron  guide 
field  at  the  equilibrium  orbit  radius  In  the  '*untrimmed'’ 
condition  is  shown  in  Pig,#^.  Ttie  corrected  field 
variation  is  also  shown  on  the  same  plot.  The  corrections 
were  made  by  adding  copper  shims  to  the  flux  concentrator 
disc  that  is  nearest  to  the  coil  mid-plane.  The  shims 
are  placed  at  those  azimuths  where  the  field  is  too 
large.  Piling  then  permits  finer  field  adjustments  to 
bo  made,  (Plg.#5),  The  amplitude  of  the  first  four 

harmonica  of  the  corrected  field  were  computed  and  are 
listed  below. 

Harmonic  # 

1 
2 

3 

4 


Amplitude 

0,^3% 

0,22% 

0.48^- 

0,21^ 
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if  we  insert  these  values  into  the  equation  for  the 

particle  displacement  from  the  equilibrium  orbit  in  a 

betatron  field  (•q#3*23 )  we  find  that  the  maximum  amplitude 

of  oscillation  will  be  X  0,02".  Allowing  a  factor  of 

(21) 

five  to  include  the  effects  of  the  windings  the  maximum 
amplitude  becomes  X  0,10"  which  is  well  within  the  limits 
of  the  vector  potential  well  (see  Fig.^)» 

All  the  magnetic  field  measurements  Just 
described  were  performed  under  the  same  experimental 
conditions  that  were  used  for  experiments  with  plasma  save 
one.  All  of  the  measurements  were  made  without  the  Sb 
windings  in  place.  The  presence  of  the  extra  copper,  for 
the  Sq  field,  that  is  required  for  runaway  production 
leads  one  to  suspect  that  the  field  conf igiu?ation  with 
the  Bq  coils  in  place  may  be  quite  different  from  the  one 
without  them, 

A  second  field-coil  together  with  flux  concentrators 
was  constructed  in  order  to  answer  questions  such  as  this 
one  while  experiments  wore  being  conducted  in  the  first 
coil.  The  only  difference  between  the  two  systems  was  in 
the  design  of  the  flux-concentrators.  The  separation 
between  successive  discs  in  each  flux  concentrator  was 
Increased  from  0,010"  in  the  first  set,  to  0,030"  in 
the  new  set.  The  gap  in  each  disc  was  increased  to  ,030", 
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the  mylar  insulation  between  the  discs  removed,  and  the 
entire  space  between  the  discs  filled  with  epoxy  to 
improve  the  bond* 

A  betatron  field  in  the  mid-plane  of  the  field 
coil  was  obtained  with  a  flux  concentrator  separation 
Just  slightly  greater  than  that  in  the  experimental  coil 
(approximately  The  vector  potential  in  the  mid-plane 

and  1/8“  above  and  below  the  mid -plane  of  the  field  coil  Is 
shown  in  Fig«#10«  The  equilibrium  orbit  is  again  at  a 
radius  of  1,9  Inches,  The  vector  pot>entlal  well  in  this 
case  is  slightly  deeper  than  in  Pig. #8,  because  the  “n" 
value  is  smaller  (close  to  zero,  but  not  equal  to  zero). 

The  effect  of  the  Bq  coils  on  the  betatron 
field  shape  was  determined  by  placing  on  the  single  turn 
coil  holder,  a  system  of  colls  that  duplicated  the  coils 
used  to  generate  the  field.  In  addition  a  type  ”c" 
network  Indentical  to  the  network  used  to  drive  the 
coils  was  placed  across  the  dummy  system.  The  resulting 
vector-potential  agreed  with  the  vector  potential  without 
the  colls  within  the  acc\irecy  of  our  measurements. 

Thus  wo  can  conclude  that  the  presence  of  the 
coil  does  not  seriously  affect  the  betatron  field 
shape,  at  least  not  the  azlmuthally  averaged  value. 
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In  the  preceeding  paragraphs  we  have  described 
In  detail  the  entire  process  of  constructing  the  magnetic 
guide  field  for  a  low  inductance  plasma  betatron,  using 
flux  concentrators.  The  techniques  used  to  roake  precise 
pulsed  magnetic  field  measurements  in  the  betatron  coil 
wer‘0  also  described.  In  the  fi  al  few  sections  we  have 
presented  experimental  data  that  show  that  the  magnetic 
field  distribution  does  indeed  satisfy  the  betatron 
condition  a  for  both  equilibrium  and  stability. 


V  ADDITIONAL  APPARATUS  AND  PLASMA  GENERATION 
5.1  The  3^  Field 

The  magnetic  field  was  generated  by  the 

discharge  of  a  low  impedance  transmission  line  through  the 

toroidal  windings  wound  on  the  vacuum  chamber.  The  line, 

(24-26) 

shown  in  Fig .#11,  which  is  designated  as  a  Type  *^0*  network, 
is  formed  by  the  parallel  connection  of  "n”  individual 
-C^  elements  each  of  which  is  resonant  at  an  odd 
harmonic  of  the  fundamental  frequency  of  the  first  section. 
This  type  of  network  has  the  advantage  that  it  incorporates 
constructively  the  internal  Inductance  present  in  all  hlgb 
voltage  capacitors.  This  feature  makes  the  line  very  easy 
to  build. 

The  design  parameters  for  the  Type  "c” 
drive  were  determined  primarily  from  the  requirement  that 
the  Bq  field  remain  essentially  constant  during  several 
half  cycles  of  the  betatron  field.  Preliminary  experiments 
on  runaway  generation  with  a  faster  6q  drive 
25  .tt  sec)  system  indicated  the  magnitude  of  the  field 
that  was  necessary  for  x-ray  production.  These  conslderatlona 
led  to  the  following  specifications  for  the  type  C-network, 


SWITCH  (5C22  HYDROGEN  THYRATRON) 


TYPE  *'C”  NETWORK 


FIG.  9^1  I 
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peak  current  = 
current  wave  form  = 
rise  t  .ae  = 
decay  time  = 
pulse  duration  = 
characteristic  Impedance  = 


5000  amp, 
trapezoidal 
10  M  sec 
10  M  sec 
100  M  sec 

.286  iv. 


The  design  synthe;  3  is  based  upon  Gullleman's 
technique  of  the  Fourier  analysis  of  the  desired  continuous 
wave  train  from  an  ideal  short-circuited  lossless 
transmlaalon  line,  le.,  the  current  thr^^ugh  the  load  la 
given  by 


00 


0 


5.2 


/. 


=  peak  current 
=  Fourier  coefficient  of 
the  n^^  section 


t 


=  pulse  length  in  seconds 
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The  Inductive  and  capacitive  values  of  each 
section  are  determined  by  the  desired  characteristic 
impedance  and  ,  namely 


The  calculated  values  of  the  b  ,  and  L- 
parameters  for  a  four  section  line  with  the  current  and 
waveform  that  was  specified  earlier  are  presented  in 
Table  #2. 


Table  #  2 
Calculated  Values 


Measured  values 


Section 

L  (uh) 

C  (uf ) 

fr(k( 

1 

7.3 

lil.0 

4.95 

3 

Q.2 

II4-.5 

14.7 

5 

11.2 

3.51 

25.5 

7 

20 

1.06 

34.8 

Because  of  the  low  frequencies  involved 
the  series  inductances  were  generally  quite  larger 
than  the  Internal  inductance  of  the  capacitors  ,  The 
additional  inductance  was  achieved  by  winding  #18  copper 
wire  on  1"-0.D,  poly  vinyl  chloride  tubing  until  the 
measured  resonant  frequency  -tf  each  section  agreed  with 
the  calculated  frequency.  The  sections  were  then  Joined 
together  with  RG-8u  coaxial  cable,  and  switched  into  the 
Bq  solenoid  with  a  $022  hydrogen  thyratlon.  The 
magnitude  of  vs  the  line  charging  voltage  is  shown 


5*2  The  Vacuum  System 

The  acceleration  chamber  was  a  pyrex  glass  torus 
with  major  and  minor  diameter  of  3.0”  ajid  5/0”  respectively, 
A  single  radial  port  was  used  both  to  pump  out  the  chamber 
and  to  provide  easy  access  to  the  chamber  Interior,  The 
chamber  was  connected  to  twS'^^nch  oil  diffusion  pumps 
placed  In  series.  With  the  second  diffusion  pump  blanked 
off  the  base  pressure  of  the  system  was  10  ^mm  Hg  measured 
at  the  pump.  Tiie  base  pressure  at  the  acceleration 
chamber  was  5^10  ^  mm  Hg.  No  provision  to  bake  the  system 
was  made  because  the  small  clearances  between  the  vacuum 
chamber  and  the  flux  concentrators  would  allow  too  much 
heat  transfer  to  the  flux  concentrators,  A  variable  leak 
connected  to  the  vacuum  system  between  the  diffusion  pumps 
and  acceleration  chamber  was  used  to  admit  the  experimental 
gas  and  control  the  backreund  pressure.  The  Initial  neutral 
pressure  was  vrrled  between  ,5~2  Torn.  (Argon  and  Krypton). 


AZIMUTHAL  GUIDE  FIELD  e_{kgaus*) 
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5.3  Plasma  Generation 

Considerations  of  maximum  beam  current  leads  to 

an  upper  bound  on  the  electron  density  that  can  be  present 

In  the  vacuum  chamber  Just  prior  to  the  application  of  the 

betatron  field.  Tlie  maximum  beam  current  In  a  plasma 

(4) 

betatron  has  been  calculated  by  Schmidt,  for  the  case  of 

zero  azimuthal  magnetic  field.  Since  the  beam  Is 

diamagnetic  the  position  of  the  equilibrium  orbit  will  be 

shifted  radially  inward  from  the  single  particle  position. 

If  N  la  the  number  of  electrons  per  unit  length  of  stream 

and  A.  the  vector  potential  due  to  the  betatron  field 
^  (4) 

the  maximum  beam  cuiTont  x  »  Is  determined  from 


5.4 


where  i)  *  position  of  the  new  equilibrium 

orbit,  and  /\  la  related  to  the  self  Inductance  of  the 
beam  via.  A'  •  Using  values  of  determined 

from  the  shape  of  the  vector  potential  on  the  Ist  and  2nd 
half  cycles  fFlg,#8)  we  find  that  the  maximum  possible 
beam  current  on  the  first  half  cycle  Is  about  70  amp.  and 
about  300  amp.  on  the  second.  Thus,  the 
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Initial  electron  density  should  be  on  the  order  of  lO^^-lO^^ 
per  cm^. 

In  order  to  keep  neutral  gas  scattering  losses 
to  a  minimum  the  backround  neutral  pressure  should  be  as 
lou  as  possible.  This  means  that  we  should  have  jH  lo”*^ 
mm  Hg.  At  these  backrouiid  pressures  the  mean  free  path 
i  for  an  ionizirg  collision  is  much  greater  than  the 
vacuum  chamber  dimensions  (d).  When  1  d  breakdown  is 
most  easily  accomplished  by  means  of  electron  mult Ipact lin V 
This  process  makes  use  of  the  secondary  emission  properties 
of  the  walls  of  the  vacuum  vessel.  Briefly  the  breakdown 
process  proceeds  as  follows.  An  R-F  electric  field  is 
applied  between  the  walls  of  the  vacuum  chamber  as 
sketched  below  (care  must  be  taken  so  that  walls  with  good 
secondary  emission  properties  are  used). 


a  b 
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(One  can  also  place  the  electrodes  Inside  the  vacuum 
chamber.  For  this  case  the  electrodes  must  be  good  secondary 
emitters).  We  start  by  assuming  that  tlere  is  a  free  electron 
near  "a".  The  initial  polarity  is  chosen  so  that  tho 
electron  is  accelerated  toward  "b".  Tlie  frequency  of  tne 
R-F  field  is  determined  by  the  requirement  that  the 
electric  field  change  sign  when  the  electron  reaches  ’’b'*. 

Then  if  the  electron  has  gained  sufficient  energy  to 
generate  more  than  one  secondary  at  "b”,  the  electrons  at 
*'b”  will  bo  accelerated  toward  "a".  In  this  manner  wo  can 
build  up  a  cloud  of  electrons  oscillating  between  the 
vacuum  chamber  walls#  Soma  of  the  electrons  make  ionising 
collisions  and  thus  a  discharge  is  gradually  built  up.  The 
process  of  electron  multlpactlon  takes  place  whenever  tho 
mean  free  path  is  greater  than  the  vacuum  chamber  dimensions 
and  the  frequency  and  amplitude  of  the  R-F  field  aatlsfy 
the  conditions  Just  described.  In  fact  in  some  cyclotrons 
electron  multlpactlon  can  prevent  the  build  up  of  the  R-F 
Dee  voltage  because  of  tho  additional  oscillator  loading 
at  resonance.  This  is  overcome  by  making  the  R-P  voltage 
build  up  sufficiently  rapidly  so  that  the  resonance  region 
is  swept  through  quickly,  or  by  applying  a  d.c,  bias  voltage 
to  the  electrodes  which  destroys  the  resonance. 

An  elementary  calculation  of  the  resonance 
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conditions  proceeds  as  follows;  we  assume  that  the  R-P 
electric  field  is  given  by  sin  wt,  Hie  equation  of 
motion  for  an  electron  in  this  field  is 


thus 


4”  “  ujt 

it  /m 

/rtf,  (V  5.5 


If  we  take  v-o  at  ts.o,  then  A=  f  .  Integrating  a 

second  time  we  get 


/f  a 


X  xiin  a/t  )  4  8 
(V  ) 


similarly  we  take  /)t  =o  at  t=o,  then  B=0.  Hence  the 
velocity  "v"  and  dleplacement  "x”  are  given  by 


/\r  - 


0^- 


/nttc^ 


h- 


U/'^) 
Uf  / 
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The  resonance  conditions  are  obtained  by  requiring  that 
the  electric  field  change  sign  at  x=d,  this  means  that 
when  x=d,  we  must  have  i.?t= 'jr  .  The  elc  '.tron  energy 
Wo  at  x=d  must  be  large  enough  so  that  the  secondary 


electron  emission  coefl'iclent  t  is  greater  than  unity 
(  is  the  number  of  secondaries  released  per  primary)* 
Substituting  into  (5«6)  we  find  that  resonance  occurs 
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when  y  satisfies  the  pair  of  equations 


J[=  24- 


/yri  (jJx 


5.7 


Solving  for  U)x  we 


-yr  \A/, 

'fa'  /m''> 


5.8 


For  pyrex  glass  Wt  ?00  volts  for  good  secondary 
emission.  Than  with  d=l  cm  we  find  that 


=  u)j  ^  150  m  c 

A  high  frequency  quadrupole  electric  field,  was 
used  to  prelonlze  the  gas  within  the  vacuum  chamber.  To 
minimize  the  effects  of  the  quadrupole  field  during 
betatron  acceleration,  the  quadrupole  electrodes  were 
designed  so  that  the  position  of  the  R-P  potential  minimum 
la  close  to  the  position  of  the  equilibrium  orbit  in  the 
betatron  f  leld .  (This  ''alculation  la  described  in  appendix 
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The  quadrup'" ■’.e  electrodes  and  acceleration  chamber  are 
shown  In  The  quadrupole  wires  were  brought  out 

through  the  gap  in  the  median  plane  of  the  field  coil  and 
connected  to  a  two  wire  transmission  line  of  variable 
length.  The  transmission  line  was  shorted  at  the  far  end 
so  that  It  could  also  be  used  to  couple  power  fion  the 
orclllator.  The  quadrupole  electric  field  has  no  effect 

on  the  runaway  behavior  since  the  x-ray  emission  time  with 

-6 

the  quadrupole  electric  field  turned  off/xlO  sec  before 
tlie  application  of  the  betatron  field  is  the  same  as  it  is 
when  the  electric  field  is  run  continuously. 


FIELD  COIL 


BLANK  PAGE 
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VI  DIAGNOSTICS 


6,1  Electron  Density 

Electron  density  measurements  in  the  R-P 
generated  plasma  are  difficult  because  conventional 
diagnostic  techniques  such  as  Langmiur  probes,  double 
Langmuir  probes  or  microwaves  (cavity  techniques  excepted) 
cannot  be  used  effectively  at  the  electron  densities  and 
temperatures  in  hhe  discharge.  Absence  of  a  reference 
electrode  prevents  the  use  of  the  single  Langmuir  probe. 
Attempts  to  measure  the  density  v/ith  double  probe 
technique  of  Johnson  and  Walter  were  frustrated  because 
the  radius  of  the  space  charge  sheaths  turned  out  to  be 
larger  than  any  convenient  probe  separation.  At  the  low 
electron  densl'  ies  present  in  the  discharge  (  'k.  10  /cm  ) 
the  required  microwave  wave  length  is  muci.  greater  than  the 
dimensions  of  "he  plasma  so  that  phase  shift  measurements 
are  Impossible. 

In  recent  years  Dattner  and  othel’^^ liave  generated 
new  int'^^rest  in  the  dipole  resonances  first  studied  by 


Tonks.  Briefly  if  one  places  a  plasma  in  an  oscillating 


electric  field,  a  resonance  takes  place  when  the  frequency 


of  the  field  w  is  related  to  the  plasma  frequency 


via.  the  equation 
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where  k  is  a  constant  that  depends  on  the  plasma  geometry. 
(k=l  for  plane,  k=2  for  cylindrical  and  k=3  for  spherical). 
There  are  other  resonances  present  but  the  one  just 
mentioned  is  the  strongest.  Thus  by  observing  the 
frequency  at  which  resonance  occurs  one  can  determine  the 
electron  density,  at  least  in  principle.  The  only 
requirement  on  the  frequency  is  that  the  wave  length 
should  be  greater  than  the  dim.ensions  of  the  plasma. 

An  experiment  was  set  up  to  see  if  the  Tonks- 
Dattner  resonance  could  be  observed  in  th.s  R-P  exdlted 
plasma.  Plasma  was  generated  by  applying  R-P  voltage  at 
120  me,  to  a  vacuum  chamber  and  quadrupole  system  identical 
to  the  system  used  in  the  field  coll.  A  strip  line  was 
placed  across  one  half  of  the  vacuum  chamber  as  sketched 
below. 
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The  3trip  line  was  then  connected  to  a  signal 
generator  whose  frequency  could  be  varied  between  ^00  and 
1000  me.  The  experiment  proceeded  as  follows.  At  a  given 
frequency  the  line  was  tuned  to  give  a  low  standing  wave 
ratio,  (V.S.W.R.)  typically  l.O^  with  no  plasma  in  the 
discharge  vessel.  A  plasma  was  then  generated  and  the 
V.S.W.R.  remeasured.  We  would  expect  to  see  a  dramatic 
change  in  the  V.S.W.R.  as  the  frequency  of  the  R-F  field 
went  through  the  plasma  resonant  frequency,  Tlie  frequency 
was  varied  from  ^00-1000  me.  in  steps  of  20  me.  No  change 
in  the  V.S.W.R,  indicative  of  a  plasma  resonance  was 
observed.  Since  the  discharge  could  not  be  run  continuously 
(because  of  vacuum  chamber  heating)  it  is  possible  that 
there  was  enough  variation  in  the  electron  density  between 
discharges  to  mask  out  the  effects  of  the  plasma  resonance. 
It  is  also  possible  that  the  frequency  steps  were  too 
wide.  If  the  "Q"  of  the  resonance  was  greater  than  50 
then  the  band  width  A-f  ='^/(3  would  be  less  than  the 
frequency  steps  used. 

In  another  experiment  the  frequency  was  varied 
continuously  from  500-1000  me,  and  the  reflected  power 
from  the  strip  line  was  recorded  on  a  chart  recorder.  No 
significant  change  in  the  reflected  power  was  observed 
between  the  plasma  and  no  plasma  conditions. 
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I.i  the  former  experiment  the  transmission 
system  Is  not  "flat”  over  the  frequency  range  ^00-1000  me. 
There  are  many  resonances  present  without  a  plasma  in  the 
vacuum  chamber.  Under  tbosts  clrciiitistances  it  is  quite 
possible  that  a  system  resonance  could  completely  overcome 
any  effect  of  a  plasma  resonance.  Thus  the  next  step  was 
to  shut  off  the  quadrupole  electric  field  and  look  for 
resonances  during  the  plasma  decay.  Efficient  detection 
was  accomplished  by  terminating  the  strip  line  with  the 
characteristic  impedence  of  the  cable,  and  measuring  the 
power  reflected  from  the  plasma  column  back  toward  the  signal 
generator,  A  directional  coupler  was  inserted  between  the 
plasma  column  and  generator  with  a  20  db  plckoff  to  a 
crystal  detector.  Approximately  10-l5db  attenuation  was 
used  between  the  directlordcoupler  and  the  U.H.F,  generator 
for  isolation.  The  120  me,  breakdown  oscillator  generates 
large  D.C,  level  signal  in  the  crystal.  The  oscilloscope, 
used  to  monitor  the  reflected  power,  was  adjusted  so  that 
the  D.C.  signal  from  the  reflected  U.H.P.  signal  is  on 
screen.  Then  the  signal  induced  by  the  quadrupole  electric 
field  (q,e,f, )  will  cause  the  D.C,  level  of  the  trace  to 
move  off  screen.  Off  resonance,  or  with  no  plasma,  a 
step  response  is  seen  as  the  q.e.f.  is  turned  off.  The 
output  changes  from  the  negative  D.C.  level  of  both  the 
q.e.f,  and  the  U.H.P.  to  the  less  negative  value  of  Just 
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the  U.H.F.  At  the  resonance  a  pulse  appears  r'^lative  to 
the  U.K.P,  level  with  no  plasma  (in  the  positive  direction). 
Initial  experiments  were  performed  in  argon  at  high  presoures, 
typically  100  microns,  in  order  to  make  detection  of  the 
resonance  easier.  Figure  #14  Illustrates  a  typical  curve 
of  electron  density  versus  time  during  the  afterglow.  The 
argon  pressure  was  held  constant  at  100  microns  and  the 
time  of  resonance  was  measured  at  frequency  increments. 

The  electron  density  was  calculated  from  the  frequency 
relation  for  a  dipole  resonance 

UJx  «  ~^-L- 

,?/m  6c  6,2 

with  corrections  made  for  the  glass  vacuum  chamber  walls. 
Results  show  a  decay  to  a  fixed  value  (not^O)  of  electron 
density  which  is  believed  to  be  that  maintained  by  the 
U.H.P,'  The  initial  decay  was  evaluated  in  terms  of  the 
volume  recombination  coefficient  (  "D  ),  If  the  particle 
loss  is  due  to  volume  recombination  only  then  the  temporal 
behavior  of  the  density  is  given  by 

-L _ L-  =3* 

0.3 

where  is  the  initial  electron  density.  If  we 
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extrapolate  the  density  data  in  Flr*//l]-i  to  zero  time  we 

9  1 

find  that  ^>c.  7,5  x  10  electrons  /cm"^.  The 

-7 

first  four  points  result  in  >  values  of  9.19  x  10  , 

••  Y  « Y 

9*43  X  10  and  9.3i|  x  10  /sec  respectively.  The 

difference  betwsen  the  ^  values  obtained  here  and  the 

-73  (28) 

value  of  6,7  i  0,9  x  10  cm  /sec  reported  recently 

for  argon  is  probably  due  to  diffusional  losses. 

At  the  neutral  pressure  used  in  this  exD(  rlment 

the  space  charge  sheaths  surrounding  Langmiur  probes  are 

small  enough  so  that  I-V  characteristics  can  be 

interpreted  without  a  great  deal  of  er,  ^r.  Thus  the 

density  during  the  quadrupole  electric  field  discharge 

was  determined  with  double  Langmuir  probes  and  compared 

with  that  obtained  via. the  resonance  technique.  The  ion 

(hence  electron)  density  was  obtained  from  the  well  known 

relation  ^  xr  /y  where  is  the  saturated 

ion  current,  A  is  the  effective  probe  area,  {  AT  )  is  the 

ion  thermal  velocity  (the  ion  temperature  is  assumed  to  be 

the  same  as  that  of  the  neutral  gas).  The  electron  density 

obtaiiied  this  way  is  about  a  factor  of  ten  higher  than 

that  obtained  with  the  resonance  technique.  As  the 

neutral  gas  pressure  is  lowered  the  electron  temperature 

increases  and  the  density  probably  decreases  thus  making 

double  probe  data  impossible  to  Interpret  and  furthermore, 


I  ELECTRON  DENSITY  DECAY  AFTER  SHUT  OFF 
j  OF  QUADRUPOLE  ELECTRIC  FIELD 
I  ARGON,  INITIAL  PRESSURE  =  100  MICRONS 
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the  resonance  becomes  more  difficult  to  find  for  reasons 
that  are  so  far  unexplained.  Because  of  the  large 
aifference  between  the  two  density  results  at  high  neutral 
pressure,  no  extensive  Investigations  were  carried  out  at 
lower  pressures. 

An  attempt  was  made  tj  measure  the  electro  i 

density  spectroscopicali f  using  the  method  of  Inglla  and 
(29,30y 

Teller,  If  one  examines  the  broadening  of  individual 
lines  of  a  series  such  as  the  Balmer  series,  one  finds 
that  the  higher  members  of  the  series  are  broadened  more 

than  the  lower  members.  Since  the  higher  members  of  the 

series  are  also  closer  together  than  the  lower  members, 
one  eventually  reaches  a  point  where  the  broadening  of 
the  lines  is  comparable  to  the  spacing  between  them,  and 

at  this  point  the  lines  in  the  series  merge  together  and 

are  no  longer  visible  as  distinct  lines.  This  process 
manifests  Itself  In  the  spectrum  as  a  depression  of  the 
series  limit  and  is  equivalent  to  saying  that  the 
perturbation  of  the  energy  levels  due  to  the  Stark  effect 
becomes  comparable  to  the  energy  difference  between  two 
adjacent  levels.  By  equating  those  two  energy  terms  one 
cani>8olve  for  the  principle  quantum  number  n^of  the  upper 
level  of  the  last  distinct  line  of  the  series.  The  ion 
density  is  then  given  by 
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/Tf 


/Hi 


6.4 


This  method  has  been  used  sucessfuliy  to  determine  the  ion 

(31) 

density  in  a  decaying  hydrogen  plasma  at  relatively 

,  3 

high  ion  densities  (  10  /cm  j.  At  low  electron  densities 

(10  /cm^)  the  upper  quantum  number  /n^  of  the  last  line 

of  the  series  is  verj  large  {  -iz,  90)  therefore  a 

spectroscope  with  large  dispersion  is  required.  After  the 

betatron  field  la  turned  on  the  electron  (  and  ion) 

density  increases,  the  discharge  "blorms”.  It  is  quite 

possible  that  the  electron  density  may  reach  a  value  of 

12  3 

10  /cm  during  this  discharge.  At  thin  density 

-Jt  20  and  hence  may  be  observable  without  resorting 
to  sophisticated  experimental  apparatus.  Although  all  of 
our  experiments  were  done  with  argon  and  krypton,  helium 
was  used  to  test  this  method  because  of  Its  simpler  spectral 
structure,  Tlie  spectral  lines  of  interest  vieve  from  the 


aeries  /w  0, 


2'P. 


(series  limit-3679,. 7  A°),  This 


series  was  chosen  because  of  its  close  resemblence  to  that 
,  (32) 

of  Hydrogen, e specially  at  high  values  of  the  principle 


quantum  number.  The  discharge  was  viewed  with  a  Hllger 
Medium  Quartz  Spectrograph  (for  this  experiment  a  glass 
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prism  .:as  used).  The  spectrum  was  recorded  on  Kodak  Type 
7  Spectroscopic  plates  and  developed  In  Acufine  to  Improve 
the  spee  The  experiment  proceeded  as  follows,  A  plasma 

was  generated  using  the  high  frequency  quadrupole  electric 
field.  The  spectroscope  shutter  was  opened  manually  just 
prior  to  the  application  of  the  betatron  and  8^  magnetic 
fields,  and  closed  again  as  soon  as  possible  after  the 
discharge.  Although  the  shutter  is  open  for  a  period  of 
time  considerably  longer  than  the  duration  of  the  betatron 
field,  the  maxlmiom  light  intensity  occurs  during  this 
time,  V/hen  the  spectrum  obtained  this  way  is  compared 
with  the  spectrum  from  the  R  -F  excited  plasma  alone,  we 
find  additional  linos  which  we  believe  are  due  to 
additional  excitation  during  the  pulsed  discharge.  If 
the  electron  density  is  in  the  region  li*^  -10^^/cm  ^  the 
last  dlatlnct  line  of  the  series  will  be  in  the  wave  length 
region  between  3757.1  and  3679.1  A°,  We  observe  what 
appears  to  be  Nitrogen  band  s tract ui*e  in  this  wave  length 
region  which  makes  interpretation  of  the  series  limit 
impossible.  Spectra  v/ere  again  taken  after  the  R-P 
discharge  had  been  run  for  some  time,  but  the  results 
were  the  same.  In  order  to  make  the  discharge  spectroscop¬ 
ically  clean  it  is  necessary  to  run  the  discharge  with  the 
and  8^  fields,  for  extended  periods  of  time. 
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rince  this  nulae  cleaning  would  have  required  extensive 
revision  of  the  electronic  controls  this  was  not  pursued. 
In  the  preceeding  paragraphs  we  have  discussed 
several  attempts  to  measure  the  electron  density  in  the 
R-F  excited  plasma.  The  reasons  for  not  obtaining  a 
satisfactory  result  for  the  electron  density  with  the 
various  techniques  were  also  discussed.  So  far  no  mention 
has  been  made  of  still  another  experimental  parameter 
from  which  at  least  a  lower  bound  on  the  electron  density 
can  be  obtained;  namely  the  induced  current  in  the 
acceleration  chamber.  The  induced  current,  1  measured 
with  the  Rogowsky  coil  wrapped  around  the  vacuum  chamber, 
is  related  to  the  electron  density  via^the  relation 


were  A  is  the  vacuum  chamber  cross-sectional  area  and 
^  and  are  suitably  averaged  velocities. 

Since  the  electrons  are  much  more  mobile  than  the  ions 
we  can  assume  that  they  carry  all  the  current;  then  we 
have 

f  = 


6,6 
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If  we  let  =C;  the  velocity  of  light,  we  obtain  a 

lower  bound  on  the  electron  density  (assuming  /♦;«  is 

uniform  over  the  entire  cross-section).  Since  the 

maxlmyrn  x-ray  yield  occurs  for  peak  values  of  conduction 

currents  in  the  vicinity  of  100  amperes  we  find  that  with 

10  3 

1=100  amp,  the  electron  density  is  greater  than  lO/cm  • 

As  we  will  see  in  section  i^6,3  the  runaway  current  only 
makes  up  a  small  fraction  of  the  total  induced  current, 
therefore  we  should  put  C  .  Since  the  azimuthal 

velocity  of  the  bulk  of  the  electrons  is  very  uncertain 
we  will  assume  that  K^Aftyxfj  c  .  with  this  assumption 
the  electron  density  during  the  acceleration  cycle  is  in 
the  neighborhood  of  lO^Vci^*'^*  This  is  the  value  of  electron 
density  that  we  will  use  in  all  future  discussions.  This 
assumption  may  not  be  as  bad  as  it  seems  at  first  because 
collective  electronic  processes  usually  Involve  the  plasma 
frequency  U)^  which  depends  on  • 

The  possibility  of  measuring  the  electron  density 
during  the  acceleration  cycle  using  double  Langmuir  probes 
will  be  discussed  in  section 


6,2  Electron  Temperature 

The  electron  temperature  In  the  R-F  excited 
plasma  was  determined  from  the  1-v  characteristic  of  the 
double  probe  date.  The  measured  T<,  was  approximately 
20  G,/0:  at  an  initial  argon  pressure  of  5  ^  10  mm  Hg. 
This  temperature  seems  to  be  much  too  large  and  could 
result  from  either  the  poor  sampling  of  the  electron 
population  or  from  the  interpenetration  of  the  probe 
sheaths.  Therefore  we  have  assumed  that  the  initial 
electron  temperature  is  ^  10  ev.  and  that  it  is  also 

10  c.v  at  the  start  of  each  of  the  successive  half  cycles 
of  the  betatron  field. 
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6.3  The  Runaway  Current 

The  runaway  current  ia  a  very  important  piece  of 
data  in  this  experiment.  With  a  knowledge  of  the  runaway 
current  we  can  determine  the  total  number  of  runaways, 
therefore  the  beam  plasma  frequency,  and  the  relative 
efficiency  of  runaway  production.  The  Rogowsky  coil  that 
is  wrapped  around  the  vacuum  chamber  measures  the  total 
current  that  la  flowing  through  the  chamber.  This  current 
is  composed  of  two  parts;  a  conduction  current  (resistive) 
and  a  runaway  current  (inductive).  One  hopes  in  an 
experiment  such  as  this,  that  the  runaway  current  will 
Indicate  its  presence  by  a  "udden  dip  in  the  total  current 
when  x-ray  emission  stfirts  (ip.,  the  energetic  runaway 
electrons  generate  x-rays  at  the  vacuum  chamber  wall), 

Ttiia  will  certainly  be  true  if  the  runaway  current  is  a 
substantial  fraction  of  the  total  current.  Unfortunately 
the  runaway  current  in  this  experiment  has  been  too  small 
to  be  detected  easily  by  means  of  the  Rogowsky  coil,  (the 
minimum  current  we  can  detect,  is  %  10  amp.).  Therefore 
we  have  had  to  obtain  the  numerical  value  of  the  runaway 
current  from  x-ray  data.  One  important  step  in  this 
procedure  is  the  measurement  of  the  azimuthal  variation 
of  the  x-ray  emission,  Fig.;^l5-l8  show  the  scintlllirtor 
pulse  height  versus  peak  conduction  current  for  four 
different  azimuthal  positions  of  the  detector. 


. 
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The  position  of  the  collimating  lens  was  chosen  so  that 
the  scintillator  would  view  only  one  quadrant  of  the 
acceleration  chamber  at  a  time.  Because  of  the  considerable 
variation  in  the  scintillator  pulse  height  from  discharge 
to  discharge,  no  attempt  was  made  to  plot  the  x-ray  signal 
verses  the  measured  peak  current,  but  instead  all  the 
scintillator  signals  that  corresponded  to  currents  within 
a  10  ampere  interval  were  plotted  on  the  same  ordinate 
(at  the  center  of  the  interval). 

The  variation  in  scintillator  signal  from 
discharge  to  discharge  does  net  necessarily  mean  that 
tho  number  of  runaway  electrons  is  fluctuating  by  the  same 
amount,  A  limited  number  of  experiments  were  performed 
with  two  counters,  each  one  viewing  a  different  position  of 
the  acceleration  chamber.  Again  there  was  variation  in 
the  signals  from  each  counter,  but  there  was  no  correlation 
between  the  pulre  heights,  le,,  if  one  counter  had  a  large 
signal  there  might  have  been  zero  signal  on  the  other  and 
vice-versa.  Because  of  this  type  of  variation  in  the 
x-ray  output  from  different  positions  of  the  vacuum  chamber, 
we  will  assume  that  f  -  average  scintillator  signal  at  any 
position  will  be  repr osertatlve  of  the  x-ray  flux  from 
that  position.  As  we  can  see  from  Figs ,;^15-18  the  x-ray 
emlaslon  varies  about  a  factor  of  two  around  the  vacuum 
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chamber.  Since  there  ie  no  azimuth  at  which  x-ray 
emission  is  extremely  pronounced  we  will  assume  that  it 
is  uniroriu.  Then  the  number  oT  photons  N  that  correspond 
to  a  scintillator  pulse  height  (Ph)  is  given  by 

/V  =  (n)  t' &' 

where  is  the  counter  response  (pulse  height  for  a 

single  photon  of  a  given  energy),  C  is  the  counter 
efficiency  and  G  is  a  geometry  factor  (essentially  the 
solid  angle  subtended  by  the  scintillating  crystal). 

The  counter  efficiency  E  was  determined  with  a 
calibrated  rftdlo-actlve  soui'ce  ,  100-101  uc,2/26/63). 

The  counting  rate  versus  pulse  height  for  two  different 
source  to  counter  distances  is  shown  in  Fig. 19.  Since 
the  1,28  mev  x-ray  peak  occurs  at  a  pulse  height  of  1,4 
volts,  the  counter  response  ^  =0.41  v/375kev  photon(the 
scintillating  crystal  was  thallium  activated  sodium  iodide 
whose  response  to  photons  is  linear  with  energy). 

We  observe  that  the  counting  rate  for  the  1,28 
mev  line  fo.  the  two  different  source  to  counter  distances 
does  not  behave  as  l/r^.  The  ratio  of  the  counting  rates 
is  approximately  3  whereas  it  should  be  4  if  the  1/r^  law 
were  followed.  We  believe  that  this  discrepancy  is  due  to 
a  too  low  discriminator  level  in  the  pulse  height  analyzer. 
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Since  this  error  will  be  small  compared  to  other  errors 
that  will  be  made  in  determining  the  runaway  current  we 
will  neglect  it. 

Tlie  counter  efficiency  E  is  defined  as  the  ratio 

of  the  counting  rate  to  the  rate  at  which  x-rays  are 

incident  on  the  scintillating  crystal.  ITie  number  of  x-rays 
22 

from  the  Na^  source  that  are  incident  on  the  crystal  is 


A/l 


r,  c/yp[- 


6.8 


where  D  Is  the  diameter  of  the  crystal  x  Is  the 

crystal  to  source  distance  (6"),  is  the  original  source 

strong  1  (100-101  u  c),  '2^,j  Is  the  half  life  of  the  source 

22 

(2,6  years  for  Na^  )  and  t  is  the  time  at  wnlch  the 

calibration  is  performed.  If  we  insert  these  numbers  into 

equation  we  find  that 

I  f.  rF 

I  (0  C 

/t  (  O 

Hi  =  LSI^  I  o'* 


Since  only  2p,8%  of  the  Incident  y*  -rays  are  stopped  in 
the  crystal 
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3.  3  ^  X  ^  t 


rays/sec  at  1.28  oev 


from  Fig.#_^  the  counting  rate  at  1.28  mev(pulse  height 
of  l.i;  volts)  la  5o00  counts/min  or  93  counts/aec, 
therefore  the  counter  efficiency  is 

E  -  -J-L-  -  ;i.y/)c/a 

3. 3  <»>/<» 


The  geometrical  factor  G  is  Just  the  fractional 
solid  angl^  subtended  by  the  crystal. 
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viT 


6.9 


where  R  is  now  the  distance  from  the  crystal  to  the 
vacuum  chamber  (a  more  elaborate  calculation,  which  takes 
the  finite  size  of  the  acceleratlot  chamber  into  account 
changes  the  G  value  from  eq,6,9  by  ouly  S%) *  Thus 


0.75)^ 


G.H  t  iO 


Finally  for  a  pulse  height  of  50  volts  (375  kev 
photo:. j)  the  number  of  x-raya  generated  at  the  vacuum 


chamber  wall  is 


,,v«  (^0)  Chi)' (i  ^ 

//  =  a.l'fJd^  photons 


Now  all  we  need  do  to  obtain  the  niomber  of 
runaway  electrons  that  correspond  to  N  photons,  is  to 
determine  the  efficiency  of  x-ray  production  by  the 
a  lectron-s , 

The  long  succession  of  gradually  Improving 
theories  of  bremsstrahlung  has  dealt  almost  exclusively 
with  the  thin  target  case.  Here  one  calculates  the 
radiation  expected  from  a  collimated  beam  of  mono-energetic 
electrons  while  passing  through  a  target  which  is  so  thin 
that  the  electrons  lose  no  appreciable  energy  by  ionization, 
suffer  no  significant  elastic  deflections,  and  suffer  only 
one  radiative  collision.  The  thick  target  case,  although 
far  less  adequately  covered  by  theory  because  of  its 
complexity,  in  nevertheless  the  usual  laboratory  situation. 
In  this  experiment  we  will  assume  that  the  x-radiations 
detected  wloh  the  scintillation  count^’^  are  produced  via. 
thick  target  brem.sstrahlung. 

The  average  bremsstrahlung  energy  emitted  in  an 

(34) 

element  of  path  length  Js  is 


where 


la  the  energy  lost  per  cm  via,  bremaatrahlung 


and  (^\  la  the  energy  lost  per  cm  via  Ionization. 

The  average  energy  J  radiated  by  an  electron  of  Initial 
energy  7  In  being  atopped  la 


I  • 


6.11 


Integration  ualng  the  theoretical  values  of  radiative 
and  ionization  losses  Is  valid  only  over  the  portion  of 
the  path  for  which  the  electrons  velocity  la  largo 
compared  with  the  velocity  of  the  atomic  electrons.  In 
most  cases,  the  integration  must  be  ’formed  graphically 
because  of  tbo  complicated  form  of  ^  * 

Tie  spectral  distribution  of  thick  target 

(34) 

brempstrahlung  from  non-relatlvlatlc  electrons  can  be 
thought  of  as  the  sum  of  the  contributions  from  a  number 
of  thin  target  cases  of  various  electron  energies  aa 
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shown  in  Flp, /20 


The  total  energy  per  frequency  interval 
proportional  to  2:  (  ^^4,^  -  ^  ),  Iq, 


13 


=  const,  h 

6.12 

If  eq,  6,12  is  integrated  over  all  frequencies  from 

'^=■0  to  =  T/h  {h=Planck '  s  constant),  we  find  for  the 

total  bremsstrahlung  energy  I,  in  mev  per  incident  electron 

j  =  Jii 


6.13 
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No  simple  dependence  of  the  x-ray  Intensity  in  any  given 

(34) 

direction  on  ^  and  E  la  found  experimentally.  The 
total  bremsstrahlung  Integrated  over  all  angles  la  found 
to  be  proportional  to  -5  ,  within  the  accuracy  of 

measurement,  and  also  to  be  proportional  to  E  within 
1  10  percent,  in  agreement  with  eq.  6,13. 

The  numerical  value  of  the  constant  k  is  knovm 
only  approximately,  the  values  ranging  between  0,3  and 

-3  -1  (34) 

1,1  X  10  (mev)  .  Following  Evans  we  will  assume  that 

I  =  0.7, 2 

(l  and  E  in  mev).  Thus  the  fraction  of  the  Incident  electrons 
energy  which  is  converted  to  bremsstrahlung  In  a  thick 
target  is  approximately 

-3 

j  0.1  IL  E 
£ 


For  glass  (S ^  effective  2.  is  Therefore 

375  hev  electrons  lose  0,26^  of  their  energy  via. 

bremsstrahlung.  Finally,  the  number  of  runaway  electrons 

7 

that  correspond  to  the  flux  of  ^xlO  photons  is 


,/  .1  %  fo 

=:  .21  t/o 

1^10 


The  runaway  current  is  (typical) 
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As  we  have  Just  seen  the  total  runaway  current 
la  quite  amall.  This  result  is  ccnsistant  with  our 
Rogowsky  coll  measurements.  One  obtain  greater 
sensitivity  with  the  Rogowsky  coil  if  one  measures 
Instead  of  the  current  i.  Tbiis  was  not  done  in  the  first 
place  because  the  signals  become  very  difficult  to 

interpret  as  soon  as  the  current  waveform  has  some  high 
frequency  structure  on  it.  A  further  complication  arises 
because  cf  the  noise  signal  that  is  generated  by  the 
and  fields.  In  our  situation  the  signals  generated 
by  this  pickup  was  approxlm.ate ly  30  times  greater  than 
Che  signal  from  the  plasma  current.  The  pickup  from  the 
and  fields  is  easier  to  eliminate  if  all  the 

signals  are  integrated. 

Although  there  were  only  two  stray  magnetic  fields 
present,  three  distinct  signals  had  to  be  removed 

from  the  Rogowsky  coil  circuit,  l.e.,  the  signals  from 

and  Bp  ,  and  the  signal  from  the  current  that  was 
induced  in  the  coll  driving  circuit  by  the  time  varying 
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8^  field.  The  latter  signal,  which  should  have  been 
proportional  to  8^  ,  was  distorted  by  the  large 

capacity  (ll+O  uf )  present  In  the  Type  "G"  network.  Tne 
effect  of  this  signal  was  ellrninatad  by  placing  in  series 
with  the  Bq  colls,  a  back-winding  which  linked  flux  from 
•  Signals  from  pick-up  loops  that  linked  flux  from 
and  separately,  were  Integrated  and  fed  In  parallel 
Into  one  channel  of  a  Tektronix.  Type  2  differential 
preamplifier.  Each  integrating  network  had  in  addition 
to  the  usual  Integrating  resistor  R,  a  small  variable 
resistor  in  series  with  R,  so  that  the  attenuation  could 
be  varied  by  a  small  amount.  The  integrated  signal  from 
the  Rogowsky  coll  was  then  applied  to  the  other  channel  of 
the  preamp,  and  the  difference  taken.  The  amount  of  flux 
linking  the  two  auxiliary  colls  wa.s  changed  until  the 
effects  of  and  Bs  on  the  difference  signal  was 

almost  completely  eliminated.  Changing  the  integration 
time  constant  with  the  variable  res  is.  helped  to  make 
the  cancellation  better.  The  cancellation  cannot  be  made 
complete  because  the  pick-up  loops  are  not  the  same 
electrically  and  therefore  introduce  some  distortion.  In 
this  manner  the  signals  from  the  stray  fields  wore  reduced 
so  that  plasma  currents  as  sma.i.l  as  10  amp.  could  bo 


detected. 


As  we  just  seen  typical  runaway  currents 

were  cf  the  order  of  one  amp.  therefore  too  small  to  be 
observed  on  the  Rogowsky  coll  current  trace.  Since  the 
plasma  current  turned  out  to  be  a  fairly  smooth  function 
of  time  any  small  kink  In  the  current  would  be  magnified 
on  the  signal.  Rather  than  going  through  the 

complicated  process  of  eliminating  pick-up,  the  Rogowsky 
coil  slgna]  was  differentiated  electronically  and  then 
displayed  on  the  C,R,0,  screen.  The  Type  jS  amplli'ler 
was  plugged  into  a  portable  preamp,  power  supply.  The 
oucput  ^rom  the  Type  ^  ,  (the  pla  .-a  current),  was  then 
differentiated  with  a  Tektronix  Type  '’0"  operationel 
amplifier.  Displaying  instead  of  1  also  failed 

to  reveal  the  rijnaway  current.  Since  ^he  peak  x-ray 
emission  occurs  at  plasma  conauctlon  currents  of  80-100 
amperes,  the  failure  to  observe  a  runaway  current  of  only 
one  amp,  la  not  surprising. 
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6,i|-  X-Ray  Measuremants 

The  most  extensively  studied  parameter  in  this 
experiment  has  been  "t^^  "  the  x-ray  emission  time,  "tx  " 
is  defined  as  the  time  delay  between  the  start  of  the 
betatron  acceleratior.  cycle  (le,  when  goes  through  zero) 
and  the  on-set  of  x-ray  emission.  It  is  important  that 
this  quantity  be  measured  accurately  so  that  the  effects 
of  varying  the  other  experimental  parameters  can  be 
determined.  The  x-ray  detector  (W<*-  -I  crystal,  RCA  6810-A 
photo-multiplier)  was  designed  with  enough  gain  so  that 
the  start  of  the  scintillator  pulse  could  be  made 
unambiguous.  Thus  a  good  measurement  of  "t^  "  now  only 
required  that  the  zeroes  of  the  betatron  field  be  measured 
accurately.  All  the  difficulties  associated  with 
measuring  the  zeroes  of  an  integrated  signal  or  with 

measuring  the  positions  of  the  maxima  on  a  /Ji 

signal  are  eliminated  if  one  uses  a  peaking  strip  pi 
placed  in  the  fringing  field  to  measure  the  zeroes  of 
.  The  peaking  probe  consists  of  a  few  turns  of 
copper  wound  on  a  short  length  of  Mo-permalloy  wire. 

The  Mo-permalloy  wire  saturates  at  only  few  orsteds  of 
total  field.  Thus  if  the  peak  time  varying  fields  are 
large  a  very  precise  measurement  of  the  zero  of  the  total 
field  can  be  made, 

'The  output  pulse  from  the  peaking  strip  probe 
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used  for  these  measuromence  has  a  full  width  at  half 
mRxii,iiOT  r  0,0^  XL  see  so  that  the  position  of  the 

peak  could  be  determined  to  within  0^02  ax  sec.  The 

scintillator  signal  was  extrapolated  back  to  zero  and  its 
displacement  from  the  peaking  strip  maximum  was  measured 
on  a  photograph  of  the  CRO  trace  with  a  calibrated  steel 
scale  and  magnifying  glass,  A  typical  scintillator  and 
peaking  strip  signal  is  shown  in  Fig./S^^,  The  total 
error  associated  with  measuring  "t^^  "  is  believed  to  be 
approximately  t  0,03  M  sec. 

Since  the  peaking  probe  measures  the  zero  of  the 
total  field  care  must  be  taken  to  make  sure  the  stray 
field  from  the  3^  coils  does  not  have  any  effect  on  the 
probe.  The  position  of  the  peaking  strip  pulse  relative 
to  the  start  of  the  oscilloscope  trace  was  measured  with 
5^  5000  gauss  ftnd  3^  =0  using  the  measuring 

techniques  just  described.  The  peak  betatron  field  was 
varied  over  all  of  the  values  used  for  x-ray  measurements. 
No  systematic  dlflarences  between  the  5^  =5k  gause  and 
=0  vconditions  was  observed,  Ttius  we  conclude  that  the 
Bq  field  has  no  effect  on  the  peaking  strip  signal, 
within  the  accuracy  of  our  measurements. 


UPPER  TRACE  =  TOTAL  INDUCED  CURRENT  !IO  amp/cm 
SWEEP  ==  2X  id'sec./cm.  BRIGHT  PORTION  INDICATES 

DURATION  OF  LOWER  TRACE. 

LOWER  TRACE  =  PEAKING  STRIP  SIGNAL  AND 
SCINTILLATOR  SIGNAL. 

SWEEP*  5X  lO^sec./cm 


TIME - »- 

PEAKING  STRIP  SCINTILLATOR 


TYPICAL  SIGNALS 


FIG. -#21 


83 


6.5  X-Ray  Absorption  Mrasurements 

The  x-ray  energies  Indicated  In  Pigs  ,lf21  - 
were  obtained  from  the  value  of  the  betatron  magnetic 
field  at  the  orbit  at  "t,y:  ”,  the  time  x-ray  emission 
begins,  via.  the  relation 

6.15 

7-  £('•3  p  ,4^  ^  yH1,C^ 

6.16 

where  e,  is  the  peak  field  strength  at  the  orbit  in 
kilo-gauss,  P  is  the  orbit  radius  in  cm,  and  oJ  is  the 
frequency  of  the  betatron  field,  (the  energies  T  and 
nig  c,^  in  eq.(6,/6)  are  in  units  of  mov).  In  order  to  make 
sure  that  the  x-ray  energy  does  indeed  correspond  to  the 
electron  energy  in  the  betatron  field  we  have  m^de  x-ray 
absorption  measurements  on  both  the  first  and  second 
half  cycles.  Fig,/;^  shows  the  variation  of  the 
scintillator  pulse  height  as  a  function  of  absorber 
thickness  on  the  second  half  cycle  for  the  conditions 
indicated.  The  error  flags  represent  the  iriaxlrn.um  and 
minimum  of  the  data  points  whose  average  is  given  by  the 
circles.  Vhe  slope  of  the  straight  lirie  drawn  through 
tVte  lolntfi  corre.sponds  tu  a  mars  absorotlon  coefllcient 
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("“/fj  "0.095  cm/gm.  The  x-ray  energy  corresponding 

to  this  absorption  coefficient  la  375  With  the 

values  of  the  experimental  parameters  chosen  for  these 
measurements  the  electron  energy  In  the  betatron  field  at 
the  time  the  x-ray  emission  starts  Is  350  kev. 

Since  there  is  considerable  variation  in  the 
data  from  the  absorption  measurements  the  two  energies 
must  be  considered  to  be  in  agreement. 

On  the  first  half  cycle  of  the  betatron  field 
interpretation  of  the  absorption  data  is  not  as  easy  as 
it  was  for  the  second  half  cycle.  The  scintillator  pulse 
height  vs.  absorber  thickness  on  the  first  hali*  cycle  is 
also  shown  in  The  slope  of  the  straight  line 

through  the  points  corresponds  to  a  mass  absorption 
coefficient  =0.58  cm/gm«  The  x-ray  energy  that 

gives  this  absorption  coefficient  is  ipO  kev.  The  electron 
energy  determined  from  t^  and  Bo  via.  eq.^5(67^  is  80  kev, 
and  hence,  differs  from  the  absorption  measurement  energy 
by  a  factor  of  two. 

No  meaning  full  error  flags  can  be  placed  on  the 
Ist  half  cycle  data  because  of  the  large  fluctuation  in  the 
x-ray  output  from  shot  tc  shot.  Signals  with  zero  or  very 
small  x-ray  pulses  were  observed  as  often  as  substantial 
ones  so  tnat  the  devl ^vlons  from  the  average  are  as  large 


SCINTILLATOR  PULSE  HEIGHT 


a.s  the  average,  ThereTore  the  neasured  peak  x-ray  energy 

is  conalstant  with  the  electron  energy  in  the  betatron 
field. 

Siunmary 

The  measured  peak  x-ray  energy  is  consistent 
with  the  energy  the  runaway  electrons  would  have  attained 
if  they  had  undergone  continuous  betatron  acceleration 
during  the  time  "t^  ",  'Plierefore  the  runaway  electrons 
were  produced  by  the  betatron  electric  field  rather  than 
by  the  electric  field  associated  with  an  instability  that 
might  have  developed  during  the  discharge. 
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6,6  Electric  Field  Measurements 

Prellmenary  measurements  were  made  of  the  radial 
and  azimuthal  electric  fields  that  develop  during  the 
acceleration  cycle.  The  fields  were  determined  by 
inserting  into  the  acceleration  chamber  two  pairs  of 
double  probes.  The  probes  were  made  from  0,02^''  dia.  copper 
wire  separated  approximately  0.075" (the  Debye  length 

.003"),  D.C,  Isolation  betv/een  the  probes  and  the 
oscilloscope  was  acconpiishod  by  means  of  a  pulse 
transformer  with  a  4*1  turns  ratio.  A  150  ohm  resistor 
was  connected  across  the  pulse  transformer  secondary 
so  that  the  primary  Impedance  was  ^  2400  ohms.  The  Debye 
length  ?!{,  41  probe  separation,  and  the  probe  impedance 
is  high,  therefore  the  secondary  signal  will  be  Just  the 
potential  difference  between  the  probes,  A  typical 
electric  field  signal  together  with  the  total  induced 
current  is  shown  in  Fig,^'*^,  Because  of  excessive  noise 
generated  by  the  capacitor  bank  discharge  all  electric 
flel.d  measurements  were  maae  on  the  2nd  half  cycle  of  the 
betatron  field.  Radial  electric  fields  as  largo  as  1000 
v/cm  have  been  observed.  Peak  values  of  the  dglmuthal 
field  are  very  much  less  than  this,  typically  50-100  v/om. 

It  is  possible  that  the  azimuthal  field  probes  are  in 
reulity  meas^arlng  a  radial  electric  field  because  of 
mechupJcal  mi3-a?vigr,ment ,  therefore  we  will  focus  all  of 


87 


ouj  attention  on  the  radial  field  moasuroments ,  Data  were 
taken  at  several  rcnial  positions  in  the  acceleration 
chamber  and  for  both  values  of  J"  but  for  only  one  value  of 
peak  betatron  field,  8q  ,  and  Ion  mass. 

The  effect  of  the  probe  is  to  decrease  the 
total  Induced  current  from  the  no  probe  condition.  The 
x-ray  emission  time  was  also  measured  with  the  probe  in 
place  and  found  to  be  the  same  as  it  was  without  the  probe. 
By  connecting  the  prlmory  of  the  pulse  transformer 
as  sketched  below  it  is  possible  in  principle  to  make  ion 
density  measurements  during  the  pulsed  discharge. 


C 


SECONDARY 


The  values  of  R  and  C  must  be  chosen  so  that  Rc.  ^  , 
the  duration  of  the  plasma  signal.  The  value  of  R  must 
also  be  chosen  so  that  Ri  44  l/  the  voltage  on  the 
capacitor.  The  maximum  value  of  R  that  could  be  conveniently 
u.3ed  was  ^  1  olim.  Even  with  this  low  value  of  R  no 

clear  Interpretation  of  the  data  could  be  made. 
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VII  RESULTS  km  ANALYSIS  OF  X-RAY  EMISSION  TIME. 


The  results  of  the  measurements  of  the  x-ray 
emission  time  "t^  "  are  plotted  in  Figs,  to  ^  •  The 

electron  energy  in  the  betatron  field  at  "t^  "  is  shown 
in  Figs.  1122  to  21*  The  error  flags  represent  the 
maximujr.  and  minimum  deviation  from  thie  average  at  each 
point.  On  both  the  first  and  second  half  cycles  of  the 
betatron  field  we  find  that  the  x-ray  emission  time 
decreases  as  the  Induction  electric  field,  Increases,  All 
the  dataore  presented  with  Bq  as  a  parameter.  The  value 
of  "t^  "  at  a  given  Bq  and  ,  Increases  with  the 

magnitude  of  Bg  until  ^  l|-000  gauss  and  then  seems 

to  be  independent  of  Bq  .  It  is  usually  the  custom  with 
electron  betatrons  to  extract  the  accelerated  electrons 


when  the  magnetic  field  at  the  orbit  has  reached  its  peak 
value.  This  is  usually  done  by  distorting  the  magnetic 
field  shape.  The  resulting  large  amplitude  betatron 
oscillations  then  cause  the  electrons  to  move  away  from 
their  equilibrium  orbit  and  strike  a  target.  In  this 
experiment  no  external  means  are  used  to  disrupt  the 
electron  acceleration.  The  energetic  electi'ons  strike  the 
vacuum  chamber  because  of  internal  processes  rather  than 
from  some  externally  applied  disturbance.  There  is  no 
x-ray  emission  for  times  between  the  start  of  the  cycle 


and  "t.^  ". 


The  x-ray  emission  begins  at  "t  ^ 


II 


and  is  moat 


INDUCTION  ELECTRIC  FIELD  E^( v/cm,peak) 
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intense  there  (Flp.,^2k^*  Subsequent  x-ray  emisaion  ^(.asta 
for  times  as  long  as  2  /t  sec  after  "t^  ",  indicating  that 
x-vB.ys  corresponding  to  the  peak  energy  of  the  machine  are 
probably  generated. 

As  we  can  see  from  the  experimental  data  there 
is  a  considerable  difference  in  the  x-ray  emission  time, 
and  hence  the  electron  energy,  for  the  l,wo  orientations 
of  Sq  relative  to  ^  . 

In  the  suceeding  paragraphs  we  present  several 
models  to  represent  the  processes  taking  place  during  the 
acceleration  cycle.  The  experimental  data  are  then 
analyzed  j.n  terms  of  these  models  with  the  hop©  of  gaining 
some  insight  into  the  processes  that  take  place  during 
the  acceleration  cycle;  the  processes  that  are  responsible 
for  the  •x’-ray  emission  and  those  that  lead  to  different 
results  for  the  two  directions  of  5q  ^ 

The  simplest  process  that  can  lead  to  premature 
x-ray  emission  comes  from  the  diamagnetic  properties  o^ 
the  beam.  As  we  have  discussed  earlier  (section 
Cl.  Schmidt  has  calculated  the  maximum  bcair  current  that  a 
betatron  field  can  hold  without  collapsing.  Application 
of  the  Schmidt  formula  to  the  vf-ctor  potential  shown  in 
Fig./irB  yields  a  maximiun  bejufi  current  of  approximately 
70  amperes  on  the  let  half  cycle  and  300  amp.  on  the 
second.  Tills  '‘.1  dilation  was  done  for  a  plasma  betatron 
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without  an  azimuthal  magnetic  field.  Just  how  the  formula 
shoulf.  be  modified  to  take  into  account  the  additional 
field  is  not  clear.  However  we  observe  no  change  in 
"t/^  "  as  the  peak  induced  current  la  varied  from  zero 
to  ov^'r  several  hundred  amperes.  The  quantity  that  does  change 
as  the  induced  current  is  varied  is  the  x-ray  emission 
intensity  and  therefore  the  runaway  current.  (  Figs 15  to 
18) 

'i'he  observation  of  the  large  fluctuating  electric 
fields  that  develop  during  the  discharge  leads  one  to 
construct  the  following  model  for  the  interaction  between 
the  runaway  electrons  and  the  background  plasma  through 
which  they  move.  Sudden  application  of  the  betatron 
electric  field  causes  some  of  the  plasma  electrons  to 
runaway.  The  electrons  are  accelerated  to  high  velocity 
very  rapidly  and  hence  move  through  tho  background  plasma 
without  collisions.  When  the  runaway  current  reaches  a 
critical  value  there  is  an  interaction  between  the 
runaways  and  the  plasma.  We  assume  that  some  sort  of  wave 
is  generated.  The  runaway  electrons  then  give  energy  to 
the  wave,  move  across  the  field  and  strike  the  vacuum 
chamber  walls.  This  model  fails  to  account  for  he 
difference  in  "t.  "  that  exists  between  the  two  values  of 

/  (  J  Moreover  since  the  x-ray 

emission  times  are  comparable  with  tho  times  that  are 


START  OF  ACCELERATION  CYCLE  \ 
PEAK  FIELD  OCCURS  Axio'liec.  LATER 


TIME 


SCINTILLATOR  SIGNALS  SHOWING  SUDDEN  ONSET 
AND  DURATION  OF  X-RAY  EMISSION 


FIG. -#34 
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associated  with  ionic  processes  (  oo  -  -x 

=  1?^-^* 

we  would  expect  "tj^  "  to  also  depend  on  the  ion  mass. 

Tiie  x-ray  emission  time  is  essentially  the  same  for  both 
argon  and  krypton  while  there  is  a  factor  of  two  difference 
in  the  ion  masses.  1110  insen'^itiv^  ty  of  the  experimental 
data  to  the  ion  mass  leads  us  to  (  nclude  that  the  x-ray 
emission  time  depends  on  electronic  processes  only. 

As  we  have  just  mentioned  the  times  associated 
with  the  onset  of  x-ray  emission  are  comparable  to  the 
ion  plasma  period  and  ion  cyclotron  period  and  consequently 
are  several  orders  of  magnitude  larger  than  the  values  of 
^ne  same  quantities  foi’  the  electrons,  llius  we  can  also 
conclude  that  collecti ve  electronic  process  are  not 
resp.nslble  for  the  disruption  of  the  betatron  acceleration, 
however  they  may  be  responsible  for  the  number  of  runaways 
that  are  generated.  In  the  following  paragraphs  we  assume 
that  the  runaway  electrons  are  moving  through  the  backround 
plasma  as  single  particles.  We  then  analy.  e  the  single 
particle  motion  more  carefully  than  was  done  in  sections 


Using  a  single-particle  model  in  the  g\ildlng 

(37) 

center  approximation  Meservey  and  Goldberg  have  made 
calculations  of  the  shape  and  position  of  particle  surfaces 
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(locus  of  guldlng"cen te ■■'3 )  for  runaway  electrons  in  a 

figure-eight  ctellarator  with  fixed  rotational  transform 

and  uniform  axial  plasma  current  density.  This  work  was 

(38) 

prompted  by  the  experiments  of  Bernstein  et  al;  who 
reported  on  the  energy  and  the  time  at  which  runaway  electrons 
reach  the  aperature  limiter  in  the  stellarator  for  the  case 
when  the  accelerating  electric  field  E  is  parallel  to  B 
(  2*  =  +  1)  as  well  as  for  the  anti-parallel  case  (  f  =-l). 

Their  experimental  results  agree  qualitatively  with  ours, 
ie,  as  the  magnetic  field  is  Increased  with  E  held  constant, 
the  x-ray  emission  begins  at  later  times.  Similarly  the 
peak  x-ray  energy  and  emission  time  is  greater  for  ^ 

=  +  1  than  it  is  for  J  =  -1.  Meservey  and  Goldberg 
have  also  made  measurements  which  qualitatively  confirm 
the  prediction  of  their  simple  theory,  namely  that,  in  a 
machine  whose  figure  eight  twist  is  negative  "runaway'* 
surfaces  drift  toward  the  center  of  curvatui-e  (inside)  of 
the  stellarator  loop  for  accelerating  field  E  parallel  to 

.Ok 

the  confining  field  3  (  |>  =  +  l )  and  toward  the  outside 
for  E  anti  parallel  to  ^  ,  The  motions  are  reversed 

for  a  macnlne  with  positive  twist  and  with  no 

machine  transform  the  asymjnetry  disappears,  and  the 
particle  surfaces  move  toward  the  outside  for  both 
parallel  and  anti  parallel  cases  .  These  results  have 
led  us  to  consider  a  similar  type  calculation  for  our 
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situation.  In  section  #3*3  we  studied  the  motion  of  a 
single  electron  in  the  combined  betatron  and  azimuthal 
magnetic  fields.  V/e  found  that  in  the  limit 
the  motion  In  the  transverse  plane  (  x,  z  )  was  composed 
of  two  elementary  motions,  a  rapid  motion  with  frequency 
close  tc  the  cyclotron  frequency  superimposed  on  a 
slow  elliptic  motion  with  1"^^, ,  The  slow 

elliptic  motion  is  Just  the  locus  of  the  guiding  center  and 
hence  traces  out  the  particle  surfaces  and  leads  to  a 
rotational  transform.  'Ihls  analysis  was  done  under  the 
assumption  of  zero  plasma  current.  If  we  assume  that  in 
addition  to  the  runaway  electrons  there  is  a  uniform  plasma 
current  flowing,  the  equations  of  motion  ( 3 . 2I4. )  are 


modified  to 


Where  i  ^ 

and  ^  is  t.io  current  density,  (These  equations  with 

(39) 

Uc  =0  agree  with  those  derived  by  Rogers).  Again  if  we 
assume  that  x  and  z  vary  as  wo  find  that  in  the 


limit 


the  motion  in  the  (x,  z)  plane  is  the 


san.a  as  that  discussed  above  but  with  the  fast  and  slow 


frequencies  now  given  by 
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r  _a*  V- 

*•1-  '-* 


T.2 


By  inspection  we  can  see  that  the  electron  motl ''n 
is  independent  of  the  direction  of  the  azimuthal  magnetic 
field  (excluding  the  usual  Larmoi  precession).  Thus  the 
presence  of  the  rotational  transform  and  the  finite  plasma 
current  density  does  not  seem  to  be  responsible  for  the 
difference  between  the  parallel  and  antiparallel  situations 
in  our  experiment. 

So  far  In  our  discussion  no  mention  has  been  made 
of  the  effects  of  the  stray  magnetic  fields  that  are 
generated  by  the  helical  windings  user  to  generate  the 
azimuthal  guide  field.  In  principle  the  currents  used  to 
generate  the  8^  field  should  lie  in  a  plane  that  is 
orthogonal  to  the  "-direction.  The  finite  pitch  of 
the  ^  vindlngs  causes  a  component  of  current  to  flow  in 
the  "  Q  "-direction  thereby  generating  magnetic  fields  in 
the  r  and  z  directions.  In  the  following  section  we 
will  discuss  the  effects  these  pertul>bing  fields  have  on 
the  electron  motion. 
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A  uniform  transverse  field  in  the  "z"  direction 
is  the  simplest  perturbing  field  one  can  have.  Let  " 
be  the  value  of  the  perturbi^tion  field.  Then  the  equations 
of  motion  become 

ff.  +  p  -/rt  ~  \  jj/.  =  -J>,  i~  \ 

Bi  ■'  Qq 

"  *  •  7  1 

/n  ~  u)t^  ^ 


Where  Bq  is  the  value  of  the  betatron  field  at  the  orbit 
radius  (  )  and  terms  containing  (  3' /  f  Lave  also 

been  dropped.  If  we  put 

^  7.4 

-  e'/a^ 

the  equations  take  on  the  familiar  form 

( t-yry  -  \  ^  u/^  ^ 

^  y 

the  effect  of  the  transverse  field  fi'  is  to. displace  the 
position  of  the  equilibrium  orbit  by  an  amount  given  by 

_  -X =.  ( ^  / feg) 

l-/n  -  d'jBa 


Since  Bq  starts  from  zero  we  can  see  that  even  small 
values  of  B  could  shift  the  position  of  the  equilibrium 
orbit  outside  of  the  acceleration  chamber.  We  recall  that 
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the  equations  of  motion  (7«3)  with  o  “0  are  only  valid 
for  times  greater  than  ,  Similarly  equations  (7»3) 
should  only  be  valid  for  i:  >t^.  Urns  is  the  most 

logical  time  to  evaluate  i’/f  ,  Refering  back  to  section 

..8 

#3«4  we  find  that  =;  7  -x  10  sec  with  =5  k  gauss, 

6 

=  1  k  gauss  (peak  field  at  orbit)  and  ji  =0,39  x  10  , 
Therefore  the  field  at  the  orbit  is  =  28  gauss. 

If  we  assume  that  B  =5  gauss,  tno  displacement  of  the 
equilibrium  orbit  will  be  approximately  2  cm,  and  therefore 
it  will  lie  outside  of  the  acceleration  chamber.  Since 

varies  with  tlma  so  will  the  position  of  the  equilibrium 
orbit.  The  velocity  with  which  the  orbit  position  nicves 
can  be  found  by  differentiating  (7.6)  with  respect  to  time. 
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If  we  evaluate  for  the  conditions  given  above  we 

q 

find  that  ^  2  x  lO^cm/sec,  therefore  under  these 

conditions  the  equilibrium  orbit  will  be  located  within 
the  vacuum  chamber  before  the  electrons  have  traveled  very 
far, 

I 

One  may  argue  that  the  value  of  Q  gauss 
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chosen  for  the  calculation  In  the  preceeding  paragraph 
is  entirely  too  small  since  it  is  only  1/10%  of  the  value 
of  the  azimuthal  guide  field.  This  question  can  be  answered 
by  actually  measuring  in  detail  the  magnetic  field  distri¬ 
bution  nroduopd  by  coils.  The  purpose  of  the 

proceeding  paragraphs  was  to  show  the  asymmetric  particle 
behavior  with  respect  to  the  direction  of  Bq  and  not  to  do 
a  detailed  calculation  of  the  single  particle  orbits.  Ihe 
effect  of  reversing  the  direction  of  can  be  seen  by 
examining  (7,6),  If  we  reverse  the  direction  of 
relative  to  ,  the  value  of  the  perturbing  field 

changes  sign,  causing  the  position  of  the  equilibrium 
orbit  to  lie  radially  inward  from  Instead  of  at  a 
radius  Jt  yjlg  (or  vice-versa).  If  the  direction  of 
as  well  as  the  direction  o£  Bff  i  '  reversed, 
thereby  preserving  the  value  of  ^  ,  the  position  of  the 

equilibrium  orolt  does  not  change  as  can  be  seen  from 
(7 *5)  (Both  B  and  change  sign,  preserving  the  sign 
of  the  ratio  b'  .  The  fact  that  the  right  hand  side 

changes  sign  does  not  introduce  as  asymmetry).  This 
effect  has  also  been  observed  experimentally,  le.,  the 
x-ray  emission  time  "t^  ”  does  not  change  if  both 
and  are  reversed. 

The  fact  that  a  small  transverse  field  can  lead 
to  a  lai^ge  displacement  of  the  equilibrium  orbit  by 


equation  7.6  still  causes  concern.  It  Is  possible  that 
there  will  be  no  electrons  available  for  cap’;ure  Into  a 
stable  betatron  orbit  at  the  time  the  position  of  the 
equilibrium  orbit  returns  to  the  interior  of  the  vacuum 
chamber.  This  will  cei’talnly  be  true  if  the  initial 
displacement  is  large  enough.  In  section  ^6.5  we  have 
shown,  via.  absorption  measurements,  that  the  peak  x-ray- 
energy  agrees  with  the  energy  that  the  electrons  should 
have  at  "t^  ",  if  they  have  undergone  continuous 

betatron  acceleration.  To  make  sure  that  a  betatron  field 
is  necessary  to  produce  runaway  electrons  in  this  device 
we  removed  the  flux  concentrators  and  searched  for  x-rays 
originating  from  the  vacuum  chamber.  Experiments  were 
carried  out  under  experimental  conditions  which  lead  to 
runaway  electron  currents  when  was  shaped  to  satisfy 
the  betatron  condition;  the  maximum  values  of  the 
parameters  were  =7.5  k  gauss,  (B^  )=7Agauss,  and 

Eq  =55  v/cm.  No  x-rays  were  detected  under  these 
conditions.  In  an  attempt  to  generate  a  more  uniform  B^ 
fieldjS  coll  support  form  was  constructed.  The  coils 
Instead  of  being  placed  directly  on  the  vacuum  chamber 
were  then  located  approximately  1/4"  away  from  the  flux 
concentrators.  T}-ie  x-ray  emission  under  these  circumstance 
was  greatly  reduced  from  that  with  the  colls  directly 
on  the  vacuum  chamber.  The  presence  of  the  extra  metal  so 


99 


close  to  the  flux  concentrators  distorts  the  betatron 
field  and  hence  to  a  decrease  in  tho  runaway  current. 

Thus  from  this  exnerlment  and  from  the  x-ray  absorption 
measurement  we  can  conclude  that  a  betatron  field  is 
necessary  for  runaway  production  in  this  experiment. 

Because  of  the  method  of  producing  the  field, 
a  transverse  field  In  considerable  excess  of  l/lO^  of  the 
main  field  might  easily  bet  present.  This  would  then,  as 
we  have  discussed  before,  place  the  equilibrium  orbit  far 
outside  the  vacuum  chamber  and  thus  prevent  betatron 
acceleration.  But  we  have  shown  experimentally  that 
betatron  acceleration  has  in  fact  taken  place,  thereby 
casting  suspicion  on  the  model  for  tho  perturbing  field 
that  leads  to  equation  (7.6),  It  is  obvious  that  a  better 
calculation  of  the  perturbing  fields  is  required.  The 
present  c  Iculation  has  value  xn  that  the  asymmetric 
behavior  with  respect  to  the  relative  directions  of 
has  been  revealed, 

Tiie  considerations  of  ohe  last  few  paragraphs 
lead  us  to  reconsider  the  model  for  the  x-ray  emission 
proposed  at  the  beginning  of  this  section.  We  begin  this 
discussion  by  first  calculating  the  longitudinal  electron 
energy  in  the  combined  magnetic  fields,  Tiie  electron 
kinetic  energy  in  the  non-relutlvistlc  limit  is  given  by 
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Vl/p  =  ^  AT.  ATg 


1 

a 


yyyi 


7.9 


In  our  linear  approximation  we  have  replaced  r  and  6  with 

0 

(r^  +  x)  and  (W  +(('  )  reapoctlvely .  “^Ince  ^  * 

to  first  order,  the  electron  energy,  i.s^, 


Wfi  ’•  ytM.  r  t  -  *1 

^  *-  J 
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If  we  put  x=b/2,  the  radius  of  the  acceleration  chamber, 
we  find  that  an  electron  will  hit  the  vacuum  chamber  wall 
if  it  loses  an  amount  of  energy  greater  than 


7.11 


For  the  present  experiment  b=0,625"  and  r  =1.9'‘  so  that 
(7.11)  reduces  to 


=  0.0  s*i 


7.12 


Next  we  calculate  the  energy  stored  In  the  lluctuating 
fields  that  develop  during  the  discharge.  The  stored 
energy  is  given  by 


lOI 
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We  assume  tha^  /£  Is  a  scaler  and  is  Just  6^  .  Then  the 
otored  electric  and  magnetic  energies  are  equal,  thus 


7.14 


where  represents  the  ave’^age  of  El^^over  the  plasma 

volume  and  V  is  the  plasma  volume.  The  magnitude  of  the 
observed  electric  field  is  ^  1000  v/cm  and  the  volume  of 
the  vacuum  chamber  is  60  cm  Hence 

no  jLiOT^fp  =.  ^oules 


The  electric  field  measurements  were  made  with  =4^00 
gauss,  =10  v/cm,  and  ^  =+l,  from  Flg.j/2^  we  find  that 
the  runaway  electrons  had  a  kinetic  energy  of  /b  1^0  kev 
when  they  hit  the  vacuum  chamber  wall.  With  these 
experimental  parameters  the  total  number  of  runaway 
electrons  4  x  10  (Fig.//^),  Therefore  the  total 

energy  given  up  by  the  runaway  electrons  Is 

Njt  =  OSH  1^/70  ^ 

.  c 


Na  =  e.ii  i-  (o 


joules 
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Thus  it  appears  that  the  energy  lost  by  the  runaway 
electrons  can  be  accounted  for  by  the  energy  stored  in  the 
electromagnetic  field.  The  radial  electric  f?eld 
measured  at  four  different  radial  positions  is  shown  in 
Fig.  ^37 .  The  fact  that  at  posit  ions  _1  and  2  the  electric 
field  persists  through  the  current  maximum  rules  out  the 
possibility  that  the  electric  field  is  of  a  (L  ) 

origin.  Very  strong  radial  fields  can  develop  if  there  is 
a  radial  collapse  of  the  electron  stream.  In  ordinary 
betatrons  the  radius  of  the  beam  shrinks  according  to 
f  oC  be  cause  cf  adiabatic  damping.  If  a  similar  type 

of  damping  occurs  in  a  plasma  betatron  large  radial  fields 
can  develop  because  the  beam  will  no  longer  be  space-charge 
neutralized  (the  ions,  since  they  are  considerably  less 
Tioblle  then  the  electrons,  will  be  left  behind,  0,5  ev 
argon  ions  will  move  only  1  mm  in  10  ^sec).  If  we  examine 
we  find  chat  the  radial  electric  field  does  change 
sign  as  the  prube  is  moved  across  the  vacuum  chamber^  This 
is  quite  clear  forJ=+l.  For  =--1  the  field  seems  to  be 
changing  sign  at  x=+,6l4.  cm.  For  x  ><?the  electric  field 
points  radially  Inward  while  for  x  ^  0  the  field  is 
radially  outward  indicating  that  there  is  an  excess 
negative  charge  around  x=0.  Thus  in  the  next  few 
paragraphs  wo  assume  that  in  addition  to  the  betatron 
and  B  fields  there  is  a  radial  electric  field  caused  by 
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an  excess  of  electrons  around  =  O 

We  assume  that  the  excess  electron  density  is 
just  due  to  the  runaways,  Tnen  if  n  c  is  the  charce  density, 


the  electric  field  at  p =(x  +  )  is  riven  by 


r  =  ’■  p 


7.15 


Thus  the  "x"  and  "z"  equations  of  motion  become 


>31  [  I  ~/n 


- 

/m  u)/ 


\  ^  7.16 


7.17 


where  G=  /j.  6^  .  B^y  inspection  we  can  see  that  the 
radial  equation  of  motion  becomes  unstable  when 


7.18 


Hence  for  stable  motion  the  charge  density  within  the  bei 


must  satisfy 


c  /f).  4  ^ 


7.19 


and  is  the  value  of  the  betatron  field  at  the  orbit. 
(This  is  not  correct  at  =B^  for  at  this  time,  the 
magnitude  of  the  betatron  focusing  forces  is  just  equal 
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to  the  magnitude  of  the  terms  In  the  expansion  of 
that  were  dropped.  As  we  shall  see  the  excess  charge 
density  is  also  zero  at  this  time.  Neglecting  both  these 
effects  v;ill  probably  not  cause  too  much  trouble.)  The 
equations  of  motion  can  be  uncoupled  by  differentiating 
(7.16)  with  respect  to  time  twice  .  Tlie  resulting  "x" 
equation  is  then  fourth  order,  i.e., 


it'* 

7.20 

where  r! 

~n  + 

■  ^ ,  ihe  slow 

frequency  Wf  in  the 

/m 

>>  1 

limit  is  obtained 

by  neglecting 

,  ie.j  oJ,  ,  is  the  solution  of 


Ll,  +  ~o 


7.21 


The  fast 
of 


frequency 


is  obtained  from  the  solution 


(  LjL  _  o 

dt^ 


7.22 


Equation  (735)  just  represents  the  electron  cyclotron 
motion  in  the  field  while  equation  (7.21)  rives  the 

slov/  motion  about  the  position  of  the  equilibrium  orbit, 
(Fig. //I),  Therefore  the  temporal  behavior  of  the  beam 
raaius  w ill  be  given  by  the  adiabatic  solution  of  (7,21), 
Equation  (75/)  is  of  the  form 
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4-  Ja^  ^  d  ,  Ja}  >  o 

7.23 


'rhe  approxlm.ite  solution  to  (7.23)  is  given  by 
(W.K.B,  solution). 


(40) 


(y^)  ^  /}  Ji  Ji 


7.24 


whore  A  is  an  arbitrary  constant.  Thus  the  beam  radius 
P  decreases  with  time  according  to 

P  V^o  7.25 


In  an  ordinary  betatron  the  boani  shrinks  leas  rapidly 
since  p  ol  ,  Let  Nj,  be  the  number  r’  runaway 

electrons,  then  N.  is  the  number  of  electrons  captured 
into  stable  betatron  orbits  and  hence  subject  to  eouatlon 
(7.23).  IIius  the  beam  density  is  given  by 


/n. 


NJ! 


7.26 


Wo  assume  that  the  beam  forms  at  t=t^  and  we  also  assum.e 
tl\at  at  t^  the  beam  radius  is  just  the  tnidlus  of  the 
vacuum  chumbor.  Then  the  charge  density  is 
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where  iS  the  value  of  the  betatron  field  at  the  orbit 
at  t^  .  increases  rapidly  with  time,  therefore  for  the 

purposes  of  this  calculation  we  can  neglect  unity  compared 
with  the  ratio.  Then  since  both  the  charge  density  given 
by  (7.27)  and  the  critical  charge  density  given 

by  {7«19)  both  vary  as  B*  we  have  arrived  at  a  condition 
on  N  jj  the  number  of  runaways.  If  we  put  ^ 
we  find  that  the  motion  of  the  runaway  electrons  becomes 
unstable  for  values  of  Nj^  that  satisfy 

/Vi  ^  Nmc  = 

^  7.28 

B|  is  given  by  ( sect ionj^3 .4) 

so  that  (7.4f)  becomes 

///,.'  f -3 ) 

For  B^  =[|..3iigauss,  b=0.625’' 

A/».  =  ^  H.  ^ 

/• 

or  runaways 

Thus  7*29  is  in  approximate  agreement  with  the  experimental 
results  for  the  number  of  runaway  electrons  ( Fig.;i5^35 )  • 
Expression  (7.29)  gives  the  correct  parametric  dependence 
of  with  B^  .  (Quantitatively  the  theoretical 
is  too  large  by  a  factor  of  2-3  for  =+l.  The 
agreement  for  ^  =-1  is  not  as  good,  the  experimental 
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/alue  is  approximately  five  cimes  larper  than  that  givsn 
by  (7,29).  Tlie  dependence  of  Njj  on  E  ^  is  not  explained 
by  this  calculation. 

(i^l) 

Harrison  in  his  analysis  of  the  acceleration 
of  electrons  in  plasmas  predicts  that  electrons  should 
runaway  freely  when  the  applied  electric  field  satisfies 
E  >7  where  Ec  is  now  the  maximum  value  of  the 
fluctuating  electric  field  that  exists  in  the  plasma.  He 
claims  that  may  be  of  the  order  of  10  v/cm,  there I'ore 
much  larger  than  the  experimentally  applied  electric 
fields.  Then  since  IL  only  a  very  small  number  of 
electrons  should  runaway,  i-*ut  we  do  not  observe  iny 
fluctuating  electric  f  ie  Ids  ( v;  1  tt>  frequencies  within  the 
probe  band  width)  at  the  start  of  the  acceleration  cycle. 

The  electric  field  begins  about  l.h  x  10  sec  later,  therefore 

v^e  do  not  believe  tiuit  Hur'rison's  mechanism^Ve spons ible 

for  tlie  small  nuip.ber-  of  runaways  observed  in  this  experiment, 

if  the  efl’ects  of  the  Induced  plasma  current 
and  the  transverse  magnetic  field  are  included  in  the 
analysis  the  maximum  number  of  runawiy  electrons  is 
given  by 


^-/licn  the  beam  forms,  1  e ,  at  t=t^  the  induced  plasma 
current  is  zero,  then  since  capture  into  stable  betatron 


orbits  occurs  at  t^  ,  equation  (7.30)  reduces  to 

N)t  I  (i-yri)  - 

/m  L 

Thus  we  see  that  a  small  transverse  magnetic  field  can  als 
cause  the  number  of  runaways  for^  =+  1  ^ o  be  different 
fromJ=-  1,  Physically,  the  transverse  field  reduces  the 
strength  of  the  focusing  forces  fc'*  one  orientation  and 
for  the  oppos  -a  aase  Increases  It.  Therefore  the  number 
of  runa;;aya  will  be  either  less  than  or  greater  than  the 
number  obtained  without  the  transverse  field. 

At  this  point  in  our  discussion  it  seems  that 
wo  should  review  the  results  of  this  chapter  so  that  we 
can  decide  the  direction  In  which  to  extend  the  analysis. 
V/o  Indicated  at  the  beginning  of  this  section  that  we 
were  attempting  to  explain  a,  the  difference  in  x-ray 
emission  time  for  the  two  orientations  of  Eg  relative  to 
.  and  b,  the  processes  that  cause  the  accelerated 
electrons  to  move  across  the  Bq  field  and  strike  the 
walls  of  the  vacuum  chambei’  without  any  externally 
applied  disturbance.  So  far  we  have  shown  that  a  small 
transverse  magnetic  field  arising  from  the  finite  pitch 
of  the  colls  destroys  the  symmetry  with  respect  to 

the  direction  of  B^  and  thus  could  account  for  the 


] 


7.31 


observed  results.  The  results  of  this  analysis  predicted 
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a  displacement  of  the  position  of  the  equilibrium  orbit 
from  the  dg  =0  <jon'^ition.  Thus  when  the  p?-ectron  motion 
became  unstable  the  time  required  to  reach  the  vacuum 
chamber  wal]  would  depend  on  the  orientation  of  Bq  , 

We  also  found  that  the  self -electric  field  that 
arises  from  the  adiabatic  collapse  of  the  runavvay  electron 
stream  tends  to  reduce  the  focusing  forces  that  are 
generated  by  the  betatron  magnetic  field.  Since  runaway 
electrons  are  captured  into  stable  betatron  orbits  Just 
after  the  start  of  the  acceleration  cycle,  the  self -electric 
field  at  this  t  ime  places  an  upper  limit  on  the  number  of 
runaways  that  can  be  accelerated.  The  number  of  runaways 
determined  from  this  analysis  was  found  to  be  in  fair 
agreement  with  the  experimental  results. 

Thus  we  have  succeeded  in  explaining  a  ,  the 
observed  difference  between  the  two  directions  of  and 
b  ,  the  reason  for  the  relatively  small  number  of  runaways. 
The  processes  that  are  responsible  for  the  disruption  of 
the  accelerated  beam  are  still  unknown.  It  also  appears 
that  we  will  have  to  modify  our  original  aasum{)ticn  that 
the  beam  disruption  was  not  caused  by  electronic  processes. 
As  we  will  discuss  shortly  the  radial  electric  field 
rises  to  its  maximum  value  in  *  10^ sec,  a  time  that  is 
more  consistunt  with  the  beam  plasma  frequency  than  any 
of  the  corresponding  ion  times.  This  could  then  account 
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for  the  observed  Insensitiveness  of  the  x-ray  emission 
time  to  the  Ion  mrss. 

As  we  have  Just  Indicated  the  origin  of  the  very 
large  radial  electric  field  In  the  acceleration  chamber 
is  still  unexplained.  From  Flg.//37  we  see  that  the 
electric  field  rises  to  Its  full  amplitude  in  approximately 

-7 

10  aec»  Thus  a  very  rapid  collapse  of  the  runaway 
stream  could  lead  to  a  beam  density  considerably  in  excess 
of  that  given  by  (7.19).  flie  motion  of  the  runaways  would 
then  be  unstable  and  they  would  move  across  the  field 

to  the  walls  of  the  vacuum  vesser..  Hence  there  are  at 
least  three  possible  ways  for  the  fast  electrons  to 
strike  the  walls.  The  electrons  can  lose  longitudinal 
energy  to  the  electromr  letlc  field  and  then  move  to  the 
walls,  they  can  be  lost  because  of  a  rapid  collapse  of  the 
minor  beam  diameter  or,  they  can  lose  energy  to  the 
electromagnetic  field  and  then  move  to  the  walls  because 
of  the  unstable  motion  in  the  resulting  electric  field. 

In  the  following  paragraphs  we  consider  processes 
whose  effect  on  the  runaway  motion  may  fall  into  one  of 
the  categories  just  described. 

One  mechanism  which  can  lead  to  a  collapse  of 
th*^  runaway  stream  is  radiation.  The  runaway  motion  in 
the  (x,z)  plane  can  be  decomposed  Into  two  oscillating 
electric  v.^,.oles,  90  out  of  phase,  one  along  the  "x"  axis 


X=-0.22"  X= +0.060" 


UPPER  TRACE  =  TOTAL  INDUCED  CURRENT  110  AMR/ CM 
LOWER  TRACE  =  RADIAL  ELECTRIC  FIELD  E^ 
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end  the  oth'>r  along  the  ”2''  axle,  I'he  time  required  for 
the  bsam  radius  to  shrink  to  of  its  initial  value 

is  given  by 


(C.G.S.  7.32 

units ) 


With  B-  =5000  gauss  and  -200  gauss 

4,  ^  (1.1)  f 

■  17,^^ — 


^  «  d.^  %  (6^  sec. 

K 


Hence  ^  Is  much  too  long  for  the  beam 
collapse  to  be  due  to  radiation. 

We  now  ask  if  the  magnitude  of  the  observed 
radial  electric  field  is  large  enougii  to  cause  unstable 
motion.  As  we  have  shown  previously  the  "x"  equation 
of  motion  is 

[ tot\ (t;*)  djL  i-  -V  )  4  = 

it^  ^  ^  dt^  7.20 

and  the  "z"  equation  is  obtained  by  replacing  x  with  2, 
The  slow  motion  about  the  B^  field  is  given  by  the 
solution  of  (in  the  \  limit) 
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^  ■#“ 


a;/ 


7.21 


where  the  effective  field  Index  n  Is  given  by 


/y)  =  •+  ^  e  _ 

We  estimate  the  value  of  following  way. 

We  chose  an  E^and  p  and  then  determine  (  /rf^  )  from 
(7.15).  With  E^=:1000  v/cm,  p  =.11  cra{ evaluated  at  t= 

•  A  /o  /  J 

1  X  10  sec);  we  find  /rj^  =  10  electrons  /cm.  Then  since 

•  (t 

B  =201  ga\„ss  at  t=10  sec. 


,6,  ‘3i 

io  t  r  lifo 

-/A.  ”  -*{ 

17.7  )i  to  ^fH.oS'flO 


/.J7 
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Therefore  since  ,i.  =1,27  >  l-/n  -x.  a.  S' 

The  electron  motion  In  the  (x,z)plane  will  be  unstable. 
For  simplicity  we  take  /r\  (  t- A\'  )=  1,  then  with  =.0 
at  the  time  the  motion  becomes  unstable,  the  solution  to 


{7.21)  Is 


4  * 


7.35 


where  J\^o  is  the  Initial  position  of  the  electron.  When 
Jr  the  electrons  hit  the  vacuum  chamber  wall. 


If  we  take  ftg  =  >4*/14..43  as  before,  we  find  that 


C  <?  5  H 


-I 


[*] 


-  /tn  Be 


C  o 


5  Q 
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4  -  1.  .H3 

4-U/  - 


1(^/6  4.  Of  jc/c'^ 


?  3A 


,  -S 

i^,  =  /^/ 


Thus  when  the  motion  becomes  unstable  the  runaway  electrons 
move  to  the  vacuum  chamber  wall  rapidly.  The  fact  that 
the  x-rays  do  not  appear  at  the  peak  electric  field  is  not 
surprising  in  view  of  the  approximations  made,  ^e,,  the 
assumed  beam  radius,  electric  field  s  trength,  etc.  Much 
more  detailed  electric  field  measurements  are  required  in 
order  to  make  this  last  analysis  more  quantitative. 

In  the  last  paragraphs  we  have  shown  that 
radiation  damping  is  much  too  slow  to  account  for  apparent 
rapid  collapse  of  the  runaway  stream.  We  have  also  shown 
that  the  magnitude  of  the  observed  radial  field  is  indeed 
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lai'ge  enough  to  cause  unstable  notion.  We  have  not  been 
able  to  shed  light  on  the  origin  of  the  radial  field,  i.e., 
does  it  cone  from  a  rapid  collapse  of  the  beam  minor  diameter 
or  is  it  the  result  of  a  longitudinal  energy  loss? 

So  far  in  this  analysis  no  mention  has  been  made 
of  collective  plasma  effects,  with  the  exception  of  the 
sell'  field  of  the  beam.  Since  we  have  not  been  able  to 
understand  "t^  .  ”  on  a  single  particle  picture  we  must  now 
look  at  the  collective  modes.  In  their  analysis  of  the 
stability  of  relativistic  self -focusing  streams  Pinkelstein 

(42) 

and  Sturrock  have  proposed  several  modes  of  instability 
end  studied  each  mode  in  turn  by  means  of  a  simple  model. 

The  first  mode  considered  5s  the  longitudinal  two  stream 
Instability,  If  aT  and  XT  are  the  velocities  of  the  ion 
and  electron  stream  respectively  and  if  the  velocities 


the  dispersion  relation  appropriate  to  infinite  plasmas 
has,  in  the  non-relativlstlc  approximation,  the  form 


(cO  - 


(uJ  ~ 


7.38 
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where  and  coj'  are  the  plasma  frequencies  for  the  ions 
and  electrons,  llie  motion  Is  stable  for  wave  numbers 
which  satisfy 

7.39 

Tho  stream  we  are  interested  in  is  non-re lativistlo 
therefore  u)^  will  be  much  smaller  than  cO,  .  If  the 
electron  stream  were  relativistic  the  opposite  case  could 
hold  since  the  mass  that  enters  into  u).  is  the  longitudinal 
electron  mass  ^  electron 

velocity  will  also  be  much  greater  than  that  of  the  ion, 
therefore,  the  stability  criterion  simplifies  to 


ky  Jk^  <c 

AT 


=:  / _ 


Since  o(  8o  and  o',  Bo  ,  the  critical  wave 

number  is  time  independent.  If  we  put  Jf^  ~  where 

-  y^Jlo  ,  the  critical  beam  density  is  n  =5.6xl0/cm 
(evaluated  at  t=tx).  The  beam  current  that  corresponds  to 
this  density  is  approximately  4.OO  amp.  This  la  more  than 
two  orders  of  magnitude  larger  than  the  observed  runaway 
current.  Therefore  we  can  conclude  that  the  two  stream 
instability  is  not  responsible  for  either  the  beam 
disruption  or  the  small  amount  of  runaway  current, 

Flnkelstein  and  Sturrock  also  propose  several 
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transverse  Instabilities  that  the  electron  and  ion  streams 
will  be  subject  to,  however  their  analysis  does  not  include 
the  strong  azimuthal  guide  field  that  is  present  in  this 
experiment.  The  non-re latlvistlc  stream  without  the 
field  is  subject  to  the  m“0  or  sausage  type  instability. 

TiiC  presence  of  the  strong  field  probably  removes  this 
Instability,  Similarly  the  sinuous  or  kink  instability 
and  the  instability  due  to  a  displacement  of  the  beams 
relative  to  one  another  will  probably  be  supressed. 

The  Negative  Mass  Instability  (N.M.I.),  as 

(43' 

pointed  out  by  Landau,  places  a  much  lower  limit  on  the 
maximum  beam  current  than  the  two  stream  instability.  A 
simple  explanation  of  the  cause  of  the  Instability  is  the 
following.  Consider  an  azlmuthally  uniform  distribution 
of  cold  particles  in  a  beam.  A  small  perturbation  in  the 
beam  density  will  cause  a  perturbation  in  the  electric 
potential  rotating  vlth  the  beam.  Those  particles  ahead 
of  the  potential  bump  will  be  speeded  up  and  those  behind 
it  will  be  slowed  down.  Those  particles  that  were  speeded 
up  move  radially  outward  because  of  the  increased  centrifugal 
force,  and  vice-versa,  1:  ■  amo-.nt  of  radial  motion  depends 
on  the  magnetic  field  shape.  For  weak  focusing  machines 
(  0  l/r\  Li  )»  the  radial  motion  is  large  enough  to 

overcompensate  the  Increase  in  linear  velocity  so  that  the 
angular  velocity  decreases.  Those  particles  that  were 


behind  the  bunp  '^re  slov/ed  down,  move  radially  Inward,  and 
hence  Increaae  their  anf;ular  velocity.  The  n<‘t  effect  is 
that  the  particles  move  toward  the  amrular  position  of  the 
potential  bump  and  the  perturbation  grows.  The  angular 
acceleration  is  opooslte  to  the  force  so  that  the  particles 
behave  as  though  they  have  a  negative  effective  mass. 

Landau  found  that  the  two  stream  N.M.I,  v;ill  b  eupfressed 
if  the  longitudinal  velocity  spread  {  )  in  the  stream 

satisf le  s 


7.1+1 


vhere  9  is  Budker's  parameter,  3,e.,  the  number  of  electrons 
per  length  of  stream  equal  to  the  classical  radius  of  the 
electron,  and  j  is  a  geometrical  factor.  This  criterion 
is  most  serious  at  injection  time,  Iherefore  if  we  put 

A  ^ 

-1  /m  (  i /M.  (  )-  A  ^  1*  and  since 

A  ^ 


•it 


4  the  limiting  lino  density  is 


7.k2 


'3c  -  O.^i  'f./O 


-S’ 


The  measured  \)  is  greater  than  by  only  an  order  of 
magnitude.  Since  there  is  considerable  uncertajnlty  in 

and  In  the  measured  v?  ive  cannot  rule  out  the  N.M.I. 


as  the  cause  for  the  low  runaway  current 
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At  the  beginning  of  this  report  we  discussed  the 

(19) 

runaway  electron  calculations  of  Field  and  Pried,  Their 
calculation  shows  that  the  runaway  current  should  Increase 
linearly  with  time  until  a  threshhold  current 

(  ^JL\  is  reached.  The  runaways  then  begin 

\  /rr?u  J 

to  generate  ion-acoustic  waves  and  the  runaway  current 
drops  to  a  very  low  value.  Therefore  if  this  were  the 
mechanism  for  the  beam  disruption  in  this  experiment,  we 
would  expect  to  see  the  x-ray  emission  start  a  short  time 
after  the  decrease  in  the  runaway  current.  In  terms  of  the 
dimensionless  variables  E=E^  ^  ~ 

runaway  current  Is  found  to  increase  linearly  with  time 
lUltil 

t  J 


'A 


7,43 


ie.,  when  the  driving  term  becomes  cornpf.rable  with  the 
ion-acoustic  wave  term.  The  ratio  of  the  measured  x-ray 
emission  time  to  the  x-ray  emission  time  computed  from 
eq,  (7^5)  is  plotted  jn  Fig,//^  for  two  values  of  ,  The 

plasma  density  used  to  evaluate  ,  was  taken  to  be  the 
measured  runaway  density  since  these  are  the  only  electrons 
that  are  able  to  make  laany  complete  circuits  of  the 
machine  without  colliding  with  the  vacuum  chamber  walls. 
With  =4900  gauss,  the  calculated  x-ray  omission  time  is 
approximately  eight  times  smaller  than  the  measured  time 
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and  for  Bg  =2900  gauss  is  only  about  four  tlr’.es  smaller, 
For^=-l  tne  number  of  runaways  increases  with  increasing 
electric  field.  The  decrease  is  large  enough  to  overpower 
E  '  and  the  net  result  is  that  ,  Increases  with  , 

This  is  not  observed  experimentally  ( Fig, #23-26  )  and 
therefore  the  comparison  for  =-l  is  not  displr  ed.  The 
electric  field  measurements  show  that  the  equilibrium  orbit 
lies  near  the  outside  wall  of  the  acceleration  chamber 
for  J  =-l,  Undoubtly  the  wall  plays  a  role  in  determining 
the  number  of  runaways  so  that  the  data  for  this  case  is 
suspect.  For  J  =+l  the  equilibrium  orbit  lies  near  the 
cancer  of  the  acceleration  chamber,  Hence  the  effects  of 
the  walls  will  be  reduced  considerably. 

It  appears  that  the  mechanism  of  Field  and  Fried 
can  explain  the  beam  disruption  and  subsequent  x-ray 
emission.  However  equation  (7.yi)  was  derived  for  the  case 
when  ^ t  4<  I  ,  In  this  experiment  St  10  so  that 

{7,^3)  does  not  strictly  apply.  Nevertheless  their  calculation 
describes  qualitatively  the  experimental  results. 

We  indicated  in  the  Introduction  that  there  are 
approximately  tv/enty  different  instabilities  that  can 
be  genarated  by  an  electron  beam  that  is  moving  throug}i 
a  plasma,  Hie  i“'3tabi  titles  are  usually  grouped  under 
two  headings,  a,  those  that  arise  because  of  Cerenkov 
excitation,  l.e.,  when  and  those  associated  with 


I?0 


the  Doppler  effect  le,  J€„a/;,  •  An  attractive 

possibility  for  this  experiment  is  the  occurrence  of 
cyclotron  overstability  in  the  electro-static  ion 

(44) 

cyclotron  wave,  'fhis  ionstability  is  generated  by  runaway 
electrons  when  the  electron  energy  satisfies 


/vf^e. } 

(mt,) 

7.44 

The  x-ray  c.nlssion  time 

doe  s 

not  change 

when  the  Ion  mass 

in  the  experiment  is  inci-eased  by  a  factor  of  two.  Thus 
it  seems  that  this  cannot  be  the  mechanism  responsible  for 
the  beam  disruption. 

In  view  of  the  need  for  more  experimental  data 
it  does  not  seem  advisable  at  this  time  to  mention  any 
more  of  the  many  possible  beam-plosma  interactions  that 
could  cone-  -vably  be  responsible  for  the  observed  beam 
disruption. 
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VIII  SUMmRY 

We  jan  now  summarize  the  results  of  this 
experiment,  'Fne  asymmetrical  behavior  with  respect  to  the 
twj  directions  of  Bg,  relative  to  can  be  explained  by 
the  presence  of  a  small  transverse  magnetic  field  arising 
from  the  finite  pitch  of  the  coils. 

Ihe  relatively  small  number  of  runaway  electrons 
is  probably  due  to  either  the  Negativs  Mass  Instability 
or  the  self  electric  field  of  the  runaway  stream. 

The  processes  chat  are  responsible  for  the  beam 
disruption  are  still  unknown,  however  the  mechanism  of 
Field  and  Pried  describes  qualitatively  the  x-ray  omission 


t  ime . 
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IX  CONCLUSIONS 

In  this  experiment  we  have  endeavored  to  make 
the  magnetic  guide  fields  sufficiently  homogeneous  so  that 
the  total  number  of  runaway  electrons  would  not  be  limited 
oecause  of  magnetic  Imperfections.  The  field-coil  geometry 
was  chosen  to  minimize  its  inductance  thus  permitting 
rapid  turn  on  of  the  betatron  guide  field.  Care  was  also 
taken  so  that  neither  too  much  plasma  was  present  at 
injection  nor  v/as  its  azimuthal  distribution  extremely 
non-uniform.  In  spite  of  these  and  other  precautions  taken 
the  runaway  current  was  only  of  the  order  of  one  ampere 
in  agreement  with  that  found  in  other  plasma  betatrons. 
Examination  of  the  experimental  data  shows  that  the 
maximum  beam  current  determined  from  the  Schmidt  criterion 
could  only  be  reached  by  raising  the  value  of  the 
accelerating  electric  field  to  several  thousand  volts/cm. 

In  conclus'  n  it  seems  that  it  is  impossible  to 
generate  large  runav;ay  currents  in  a  plasma  betatron  unless 
some  means  can  be  invented  to  overcome  the  effects  of  the 
self  electric  field  and  or  the  Negatlve-Mass-Instabi llty , 

If  these  effects  can  be  overcome  the  electron  energy  at 
"t^  "  may  be  sufficiently  high  to  be  of  interest. 
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XI  APPIlNDIX  A1 

In  this  section  we  calculate  the  transient 
penetration  of  a  magnetic  field  Into  a  conducting  material. 
For  computational  ease  we  consider  a  semi-inf inito  half 
space  of  material  with  conductivity  <r  and  permeability^, 
Tne  first  problem  we  consider  is  that  of  a  constant  magnetic 
field  applied  at  t=0.  The  diffusion  of  into  the 

conducting  half-space  is  riven  by  a  diffusion  equation. 


o 


11,1 


The  aolutlon  of  this  equation  subject  to  the  condition  that 

=  o  t  lo 

t  >0  11,2 


1b  given  by 

Let  us  now  compute  the  time  (tp  )  at  which  the  magnetic 
field  at  /f-  =  skin  depth,  is  equal 
to  ^  the  value  of  the  field  at  the  surface. 


^  <3.  b 
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/ 


-  O.Qi 


(t 

Thus  with  =  .39  X  10  /K»ft/sec 

Since  the  quarter  period  of  the  betatron  field 
Is  only  )L  sec,  the  step  function  maf^netlc  field 
approximation  Is  not  a  rood  one.  Wo  now  solve  the 
diffusion  equation  with  the  conditions  that 

-  o  ^0 

~  Uft  it"  >  ^ 


We  have  seen  that  the  function 

SjO.t)  -  8,  [i  - 


11.3 


where  K  = 


^  ,  has  the  property  that  B  (  t  )  =0  and 

P 


Thus  the  function 


‘  -  =Y[7feT"‘] 


-J 


11.4 


is  the  function  which  satisfies  the  one  dimensional 
diffusion  equation  and  the  conditions  W  t  o  ,  id/  )=0, 

1  ).  'Ihorcforo  by  applyinn  Duhamel's 

^  (45)  ^ 

Theorem  tlie  solution  to  the  boundry-value  problem 
(4^0  )=0,  P.  ^  (  o^i:  )=  j  {  -t  )  is  given  by 


'X. 

--  jtJ  83  </<•' 


11.5 


With  ?  (  )  =  n  j,  >4^ 


/  B,  (; 


11,6 


If  we  carry  out  the  dlf ferontlation  under  the  integral 
sign  and  let  ^  -  cJt  ,  /f  =  cc>  t  $  ^  ~  ^ 
and  p  = 


the  magnetic  field  is  given  by 


n 

S  J 


I 

d4 


11.7 


liio  integral  can  be  ovaliuitod  numerically.  This  has  been 
done  for  p  =1,  The  result  is  plotted  in  Fig.//J9  together 
witli  the  steady  state  solution 

Ba  ^  K  G  ^  (yf-  P  k7) 


11.8 


BLANK  PAGE 


The  factor  occurs  because  the  most  convenient 

definition  for  the  skin  depth  in  the  transient  analysis 
differs  by  this  factor  from  the  usual  definition  employed 
in  the  steady  state  analysis,  •  Wf’  SQ®  from 

the  plot  that  during  the  first  quarter  period  of  the 
applied  field  the  penetration  is  substantially  less  than 
is  predicted  on  the  basis  of  the  steady  state  solution. 

This  gives  us  some  confidence  that  the  field  shape  during 
this  time  has  a  weak  time  dependence.  However  by  the 
end  of  the  third  quarter  period  of  the  applied  field  the 
penetration  is  well  described  by  the  steady-state 
solution.  Thus  we  would  expect  important  differences 
between  the  field  shapes  at  the  end  of  the  first  and  thj.rd 
quarter  periods,  but  not  between  the  third  and  fifth 
quarter  periods.  This  is  in  agreement  with  our  measurement 
(see  Fig.j(^8).  These  results  show  that  in  a  coil  employing 
flux  concentrators  for  field  shaping,  an  A„C.  steady  state 
method  of  mapoing  the  field  will  not  give  an  i.ccurate 
picture  of  the  transient  field  shape  occurlng  during  the 
first  quarter  period  following  turn  "on”;  a  transient 
measurement  of  the  field  shape  must  be  carried  out. 
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XII  APPENDIX  //2 

In  t!iis  section  v;e  calculate  the  error  in  the 
single  turn-coil  signal  that  arises  from  a  lateral  shift 
In  the  colls.  We  suppose  that  the  center  of  the  pick-up 
coll  Is  displaced  a  distance  "a"  from  the  geometric  axis 
of  the  coil.  We  would  like  to  know  the  error  In  the 
measured  signal  that  Is  produced  by  this  displacement. 

The  voltage  Induced  In  the  single  turn  coll  Is 

given  by 


vjhe  re  the 
elemental 


Ir.tegratlon  Ik  cui'rled  out  along  the 
path  of  integration  Is  given  by 

h  --  Va’  +  i9 


coil. 


12.2 


The 


Therefore 


jar 


v-=- 


r 


Eq 


-lx 


<19 


12.3 


'lx 


With  C(ho4^-\j  -  and  the  equation  for  the 

displacement  circle  a  Wtff'^Txpanded 

to  first  order  In  ’’a',  eq,/^3  becomes 


JTf 

V  =  f  [  /  -  ^-  J$ 

J  ^  12,14. 

0 


Since  4  =  ^  c*/^sT.  we  may  write  for  ^ 


^0  -  I  '*'  i 


12.5 


where  Ao  is  the  valuf  of  the  vector  potential  at  the 
equilibrium  orb.t,  "n"  is  the  field  index  and 

/%Z  <i,0H>e.x  .  I'hus  to  second  order  in  "a"  the 

induced  volt ape  is 


an 


0 


j,  (t-'n)  Jte 


12.6 


12.7 


The  signal  from  the  undisplaced  loop  is  t/,  ~ 

Hence  the  fractional  error  introduced  by  the  shift  is 


i£. 


i/-4_ 

K 


—  -t”  ( I  }  'A — 

i  ^  ^  Jt,*' 


Typical  lateral  displacement  are  it  0,005'%  therefore 
with  Jf,  =1,9",  and  n=0 


i±l 


^  1- 
H 


-<i  -1 


j  r/<'g 
d  il 


=  i.l  %  >0 


-i 


ITius  we  see  that  small  lateral  displacements  of  the  single 
turn  pick-up  coils  will  not  introduce  significent  errors 
into  the  measurement  of  the  radial  distribution  of  the  vect 


potential. 
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XIII  APPENDIX  //3 

Electr'  Quadrupole  Coll  Design 

A  I’our-coil  electric  quadrupole  system  has  been 
designed  to  provide  the  oscillating  R.F,  electric  field 
for  creation  of  a  low  density  plasma  In  the  vacuiun  chamber. 
The  electric  field  resulting  from  four  ring 
charges  (of  evqual  absolute  magnitudes)  in  a  quadrupole 
array  will  vanish  along  a  circular  lino  In  the  median 
plane.  We  want  this  line  to  coincide  v;ith  the  location 
of  the  minimum  in  the  betatron  field  vector  potential. 

The  expression  for  the  median  plane  value  of  the  z  component 
of  the  quadrupole  electric  field  was  computed  as  a  function 
of  the  radial  position  p  measured  from  the  axis  of  the 
field  coll.  The  position  of  the  null  (  )  was  then 

determined  by  interpolation.  Values  were  assigned  to 
and  f  the  distances  of  the  colls  above  and 

below  the  median  plane,  and  to  ,  and  <?a  the  radii  of 


the  coils,  and  Ej  was  evaluated  for  various  combinations 
of  Pi  a,,  '^,j  j  from  the  expression 


13.1 
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where  k  =  _  and  E(k) 

("a-*  p)'  -* 

flrat  kind.  The  reaulta  of 


ia  an  elliptic  Integral  of  the 
the  computation  are  plotted 


in  Fig.  #40. 


RADIUS p  (INCHES) 


Q/  (INCHES) 
FIG.  ^40 
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EARLY-TIME  COLLECTIVE  EFFECTS 
IN  A  PLASMA  BETATRON 


S.J.  Lukas ik 
K.C.  Rogers 


ABSTRACT 


Electrostatic  collective  effects  occurring  during  startup  have  been 
examined  numerically  in  a  simple  model  of  a  plasma  betatron  having  a  uni¬ 
form  guide  field,  vector  potential  A^(r,t)  =  0.5  [r  +  1.3  Tj^/r]  slnout. 
A  region  in  the  neighborhood  of  the  equilibrium  orbit,  populated  by  an 
initially  uniform  zero  temperature  distribution  of  electrons  and  infinitely 
heavy  ions,  was  divided  into  cells.  The  trajectory  of  each  particle  of  a 
cell  was  taken  to  be  that  of  the  central  particle.  The  electrostatic  force 
acting  on  a  particle  during  a  time  step  was  obtained  by  extrapolating  the 
charge  distributions  obtained  during  preceding  time  steps.  Typical  values 
of  parameters  were:  electron  density  0  s  n  s  cm'^,  B  »  50  kilo- 
gauss,  uj  -  10®  sec"  ,  r,  =  5  cm.  For  n  ^  10®  c^^^  the  electrostatic 
forces  on  the  electron  beam  are  negligible  and  pure  betatron  motion  ensues. 
Radial  betatron  oscillations  are  observed  and  the  electron  beam  slowly 
collapses  to  the  equilibrium  betatron  orbit.  For  n  lo’^  cm"^,  the  elec¬ 
trostatic  forces  determine  the  radial  motion  of  the  beam.  The  collapse  of 
the  beam  to  the  equilibrium  betatron  orbit  is  completely  Inhibited.  Radial 
plasma  oscillations  ^  (ne^/e^m) excited  by  the  betatron  oscil¬ 

lations  =  (e/m)B^  sinu)t  are  observed.  Intermediate  values  of  charge 
density  lead  to  beam  behavior  that  shows  a  gradual  transition  from  the 
pure-betatron  case  to  the  pure-plasma  oscillation  case.  In  all  cases  the 
plasma  oscillations  were  found  to  be  stable.  The  presence  of  a  sufficiently 
strong  field  completely  suppressed  the  radial  plasma  oscillations. 
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Early-Time  Collective  Effects  in  a  Plasma  Betatron 

S.J.  Lukasik  and  K.C.  Rogers 
Stevens  Institute  of  Technology,  Hoboken,  New  Jersey 

I  INTRODUCTION 

Previous  calculations  related  to  the  early-time  conditions  In  a 
plasma  betatron  were  concerned  with  the  maximuin  allowable  temperature  of 

the  plasma  that  would  lead  to  the  capture  of  electrons  into  stable  betatron 

(0 

orbits.  '  The  study  consisted  of  calculating  single-particle  trajectories 
for  a  variety  of  initial  kinematic  condit'Ons  and  magnetic  field  shapes  and 
observing  for  what  values  of  the  various  paiameters  the  electrons  escaped 
from  the  region  of  the  vector  potential  well.  For  these  calculations  all 
collective  effects  such  as  the  electric  and  magnetic  seif-fields  of  the 
beam  were  ignored. 

In  the  work  reported  here,  the  effect  of  the  collective  electro¬ 
static  field  that  develops  during  the  start-up  of  a  plasma  betatron  is 
studied.  In  so  doing,  the  initial  temperature  of  the  plasma  has  been 
largely  ignored,  being  assumed  to  be  0°K. 

A  self-consistent  calculation  is  performed  in  the  following  way. 

A  finite  width  beam  Is  divided  into  a  number  of  "macroparticles"  whose  sub¬ 
sequent  motion  is  taken  to  be  identical  to  that  of  a  single  particle  located 
at  their  center.  At  any  time  increment,  the  electrostatic  force  on  a 
particle  is  calculated  using  the  negative  charge  density  distribution  result¬ 
ing  from  the  |X>sitions  of  all  of  the  macroparticles  at  the  end  of  the  previ¬ 
ous  time  increment.  The  positive  charge  is  assumed  to  remain  fixed  in  its 
initial  uniform  distribution.  A  Runge-Kutta  technique  is  used  to  advance 
eacn  of  the  macroparticles  by  a  time  increment.  As  the  new  radial  position 
of  each  macroparticle  is  calculated,  it  is  distributed  into  radial  accumula¬ 
tors,  the  sum  in  each  becoming  the  charge  distribution  to  be  used  during 
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This  approach  has  been  utilized  recently  to  study  several  problems 
involving  collective  effects  In  plasmas.  Bnneman^^^  and  Dawson^^^  con¬ 
sidered  the  large  amplitude  llm’t  of  longitudinal  plasma  oscillations,  the 
thermal izat ion  and  ergodic  behavior  of  an  equilibrium  plasma,  and  the 
development  of  the  two-stream  Instability  by  means  of  a  one-dimensional 
model.  Similarly,  Burger  '  has  studied  the  charge  distribution  in  a  one- 

(5) 

dimensional  plasma  diode  with  arbitrary  boundary  conditions.  Hockney'  ' 
has  examined  the  oscillations  of  a  two-dimensional  plasma.  Dunn  and  Ho^^^ 
have  calcjlated  the  approach  to  equilibriLmi  in  a  one-dimensional  d‘‘ode  as 
the  result  of  an  electron  beam  traversing  and  continually  ionizing  a  back¬ 
ground  gas.  in  all  of  these  numerical  calculations  the  object  has  been  to 
study  the  evolution  in  time  of  a  system  consisting  of  a  number  of  charge 
elements  each  of  which  is  acted  on  by  the  long-range  electrostatic  forces 
of  all  of  the  other  charge  elements. 
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II  EQUATIONS  OF  MOTION 

The  non-relatlvlstlc  Hamiltonian  function  for  a  single  electron  in 
an  axial ly"Symmetr Ic  magnetic  field  Is  given  in  cylindrical  coordinates 
(r,cp,2)  by 

p  ^  P 

This  leads  to  the  equations  of  motion 


(2) 

=  mv  -  erA  =  const 
cp  cp  cp 

0) 

"  “ '  If  4;  [('“z  '''z)*  ♦  (f^  * 

(4) 

where  the  bracketed  quantity  In  eq.  (4)  plays  the  role  of  an  effective 
potent ial 

[(P^  +  eAj®  t  (^  +  eA  )2] 


ef  f  an 


For  this  work  we  have  chosen  the  following  model  for  the  magnetic 


field 


(a)  the  axial  component  of  the  magnetic  field  B 


tp' 

and  the  associated  component  of  the  vector  potential  A  is 

9 

shown  in  Fig.  1.  The  field  has  a  value  Bj  slnuit  out  to  a 
radius  and  a  value  B^  sinout  between  r^  and  r^  . 


For  the  region  r,  r  s:  r. 


,  A^  Is  given  by 


I 


A  (r,t)  = 


1 


^  Pgr  "l^ 

B2(r',t)dr'  «=  ^  1/  (P"0  — J  slnwt 


(5) 


wii— j  »i  ^ 


■-t 


(b)  the  azimuthal  component  of  the  magnetic  field  also 

shown  in  Fig.  1  is  independent  of  time.  By  Ampere's  law 

=  B,r,/r 


9 


^  r 


and  hence 


A^(r)  = 


B^(r')dr'  =  -  q  In  r/fj 


(6) 


where  q  =  B^/B^ 

The  azimuthal  canonical  angular  momentum,  P  ,  can  be  chosen  arbitrarily. 

cp 

The  axial  canonical  momentum,  ,  can  be  expressed  using  '^q.  (2)  as 

'"z  (7) 

where  r(0)  indicates  the  initial  radial  position  of  the  electron  and 
where  has  been  set  equal  to  zero  initially. 

The  equations  of  motion  can  be  conveniently  non-d imens iona 1 ized  by 
introducing  the  following  quantities: 

u)  =  angular  frequency  of  the  betatron  guide  field 

B(t)  =  B  sinujt 
'  o 

=  electron  cyclotron  angular  frequency  in  the  orbit  region 
^o  ”  ° 

X  =  u;t 

p  = 

R  =  r/rj 

3  ”  u-u/ ^ 

o 

2 

a  =  uj  r , 

cp  o  1 

In  temjs  of  these  quantities,  the  radial  and  azimuthal  equations  of  motion 
(3)  and  (4)  for  1  ^  R  ^  p- 1  can  be  written,  using  eqs.  (5)-(7), 
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sin  x]  -  f  R  sin®  X  -  in  (^)  (8) 

*  ”  ^  [“  +  H"*  *  0) 

R 

These  can  be  solved  numerically  using  a  Runge-Kutta  fourth-order  technique. 

Writing  the  radial  equation  (8)  as  two  simultaneous  first-order  equations, 
we  have 


“  S(x,R,S) 


--i 


R  s  i  n  ’  X  -  5.^  )  n  ( 


^oT 


T(x,R,S) 


subject  to  the  initial  conditions 


R  *  R 


S  =  S 


X  =>  X. 


The  evaluation  of  the  azimuth  9(x)  is  simple  when  R(x)  Is  known  since 
all  of  the  quantities  on  the  right-hand  side  of  eq.  (9)  are  then  known. 

According  to  the  Runge-Kutta  procedure  that  has  been  used,  to  fourth  order 
in  Ax, 

AR  =  (a'  +  +  35," '  +  a''')/6  n 


where 


.  '  ' '  '  '  '  u, 

AS  =  (6  +  26  +  26  +  6  )/6 


S(x^  +  t  a  5')^ 


A  '  ’  '  r  /  I  1  ‘  '  I  '  ' 

A  -  S(x^  +  ^,Rq  +  ■*'  2  ^ 


, ,  < « «  III 

S(x  +  Ax,R  +  A  +  6  )Ax 


(1 5a) 


(16'a) 


(1 7a) 


(18a) 


6  =  T(x^  +  +  -^  ,S^  +  ^  6  )A>c  (l6b) 

6  =  T(x^  +  +  -^  ,S^  +  -^  6  )Ax  (17b) 

6*'^  =  T(x^  f  Ax,R^  +  A  ,S^  +  6  )Ax  (l8b) 

Thus,  to  advance  one  time  step  Ax  ,  one  must  calculate  AR  and  AS  using 
eqs.  (13)  and  (14).  This  if  done  by  evaluating  eqs.  (15a)  and  (15b),  (l6a) 
and  (l6b),  (17a)  and  (17b),  and  (l8a)  and  (l8b)  in  that  order. 

Equation  (8)  is  an  expression  for  the  radial  acceleration  of  an 
individual  particle;  I.e.,  it  Is  the  radial  force  per  unit  mass  on  a  particle 
having  as  initial  conditions  eqs.  (12).  Suppose,  however,  that  a  region  of 
space  is  filled  with  a  neutral  plasma  all  of  whose  particles  are  Initially 
accelerated  as  described  by  eq .  (8);  due  to  the  varying  trajectories  of  the 
electrons  and  ions  as  a  result  of  the  differences  in  mass  and  initial  condi¬ 


tions,  there  will  be  a  tendency  for  an  electrostatic  space  charge  to  develop. 
The  resultant  electrostatic  field  will  then  influence  the  subsequent  motion 
of  the  individual  particles  making  up  the  plasma.  The  object  of  the  present 
calculation  is  a  self-consistent  solution  of  eq.  (8)  augmented  by  the  elec¬ 
trostatic  space-charge  force. 

To  derive  this  additional  force  on  the  particle,  note  that  by  Gauss's 
law  the  electrostatic  force  for  a  charge  distribution  having  radial  symmetry 
and  uniform  in  the  z  direction  is  purely  radial  and  is  given  by 


E 


r 


q .  /2r!  e  r 
^  I  n  o 


where  q.^  Is  the  net  positive  charge  inside  the  radius  r  .  Then 


acceleration  =  e  E  /m  =  e  q ,  /2rr  m  e  r 

r  ^  I  n  o 


where,  for  unit  height  in  the  z-direction 


-  7  - 


q .  =  e 

^  I  n 


o 


I  I 

(N  -  N_^)  2  TT  r  dr 


The  N's  are  positive  and  negative  charge  densities.  The  N‘s  are  written 


as 


so  that 


N  *  n  N 
+  +0 


N  “  n  N 
-  o 


F  2  N  0 

e  _  e  o 

m  m  e  r  r 
o 


I  I 

(n  -  n  )  r  dr 


The  quantity  represents  the  initial  uniform  charge  density  per  unit 

volume  and  hence  the  n's  represent  the  charge  density  relative  to  the 
Initial  charge  density.  Taking  the  plasma  f'"equency  u)p  as 
2  2 

uj  =  e  N  /m  e 

p  0  0 


09) 


the  non-dimensional  radial  acceleration  due  to  the  electrostatic  force  is 


d^R 


,2 


I  I 


dx 


2  5r  “  ^  J  (n.  -  n^)  R  dR 


^  lelect  '"l 


(20) 


where 

2  2,2 

Y  =  u)p  /u) 

ThuS;  electrostatic  forces  on  the  motion  of  any  particle  can  be  included  by 
adding  the  right-  nand  side  of  eq.  (20)  to  the  right-nand  side  of  eq.  (8). 
This  amounts  to  replacing  eq.  (11)  by 


^  =  T(x,R,S)  =  [a  +  sin  x 
R 


2  2  2 

R  sm  X  -  a-£.  In  (^)  + 


2  n 


(n_  -  n^)  R  dR 


(21) 


rngumfuwmmB- 


•ilH»  I 


-  8 


The  total  negative  charge  present  is 

00  00 

I  pill 

dr  “eN2n  n(r)r  dr  (22) 

®  ».  ” 
o  o 

The  significant  difference  between  eqs.  (II)  and  (21)  is  that  the 
former  describes  the  trajectory  of  a  single  particle  in  a  given  magnetic 
field  while  the  latter  includes  the  collective  effect  of  the  instantaneous 
positions  of  all  of  the  other  particles  on  the  trajectory  of  a  single 
particle.  The  nr-^erical  solution  of  eq.  (21)  requires  a  fundamentally  dif¬ 
ferent  procedure  since  now  it  >s  necessary  to  determine  simultaneously  the 
trajectories  of  an  ensemble  of  particles. 


’total 


! 

U 


Mi  COMPUTATIONAL  PROCEDURES 
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The  calculation  of  the  sarSy-tlme  behavior  of  a  plasma  betatron 
assumes  an  Infinitely  long  cylindrical  plasma  annulus  filled  uniformly  with 
electrons  In  a  fixed  matrix  cf  positive  charge  as  shown  In  Fig.  2(a).  The 
total  Initial  width  of  the  plasma  region,  hereafter  referred  to  as  the  beam, 

Is  divided  Into  a  number  of  equal  radial  width  sections  called  macropar': ides. 
Each  macroparticle  thus  represents  a  part  of  the  initial  uniform  charge  dis¬ 
tribution.  Throughout  the  subsequent  motion  the  width  of  each  macroparticle 
is  kept  constant.  The  motion  of  each  macroparticle  Is  taken  to  be  that  of 
a  characteristic  electron  located  at  the  mid-radius  of  the  macroparticle. 

Thus  all  the  charged  particles  represented  by  a  macroparticle  are  assumed 
to  follow  the  same  trajectory. 

For  the  calculations  reported  here,  it  has  been  assumed  that 
p  »  Bj/Bg  =  2.3  .  This  means  that  the  vector  potential  well  has  a  (dlm''n- 
slonless)  extent  from  an  inner  radius  R  »  1.00  to  an  outer  radius 
R  =  1.30  .  The  Initial  charge  distribution  has  been  restricted  to  a  region 
near  the  minimum  of  the  vector  potential  well  at  R  =  1.14  whose  width  is 
only  20  percent  of  the  wMth  of  the  well.  This  avoids  complications  arising 
from  the  possible  loss  of  negative  charge  from  the  beam  as  a  result  of  large 
amplitude  transverse  oscillations. 

Since  the  positive  charge  is  assumed  to  remain  fixed  at  its  Initial 
position,  only  the  negative  charge  density  need  be  carried  in  the  calcula¬ 
tions  as  an  explicit  variable.  An  initial  uniform  charge  density  is  specl- 
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fied  by  the  value  of  y  which  Is  proportional  to  ,  the  number  of 
charged  particles  per  unit  volume.  The  radial  dependence  of  the  charge 
distribution  Is  given  by  the  values  of  the  relative  negative  charge  density 
n  in  a  number  of  equal  width  radial  cells.  These  cells  form  a  fixed  radial 


network  that  assists  in  the  charge  "accounting"  that  must  be  performed 


throughout  the  calrulation.  This  network  or  mesh  extends  from  some  lower 
radial  limit  that  is  less  than  the  inner  radius  of  the  beam  to  an  upper 
radial  limit  that  is  greater  than  the  outer  radius  of  the  beam.  It  is 
Important  for  the  successful  operation,  of  the  program  that  the  mesh  always 
be  at  least  as  extensive  as  the  radial  charge  distribution  of  the  beam. 

Since  the  motion  of  a  charged  particie  depends  on  the  radial  depend¬ 
ence  of  the  electrostatic  charge  distribution  (see  eq.  (21)),  It  Is  neces¬ 
sary  to  specify  this  with  as  much  accuracy  as  possible  On  the  other  hand, 
running  time  and  computer  storage  Ccoacity  place  upper  limits  on  the  time 
and  space  resolution  of  the  calculation.  The  time  resolution  Is  determined 
by  the  ;-ize  of  the  time  step  while  the  space  resolution  is  determined  by 
both  the  number  of  macroparticles  into  which  the  beam  is  divided  and  by  the 
radial  mesh  used  ir.  descfibing  the  charge  density.  For  most  of  the  work 
reported  here,  the  beam  was  divided  into  six  macroparticles  and  the  radial 
mesh  used  to  describe  the  charge  density  had  a  spacing  one-tenth  that  of 
the  width  of  a  macropart tc !e .  These  geometric  relationships  are  shown  in 
Fig.  2(b). 

On  the  basis  of  this  description,  the  details  of  the  calculation  that 
are  of  importance  for  evaluating  the  results  can  be  examined.  It  is  especi¬ 
ally  important  to  be  able  to  distinguish  between  computed  results  that  are 
characteristic  of  the  physical  system  under  investigation  and  results  that 
are  simply  characteristic  of  the  numerical  procedures  involved. 

The  first  step  indicated  in  the  flow  chart  in  Fig.  3  is  to  read  in 
the  parameters  needed  to  specify  the  calculation.  Next  the  initial  condi¬ 
tions  for  the  first  time  step  must  be  established.  These  are  of  two  types: 
the  initial  charge  distribution  in  the  beam  and  the  initial  kinematic  con¬ 
ditions  for  the  macroparticles. 


Setting  the  initial  charge  distribution  consists  of  specifying  the 
width  of  each  of  the  radial  intervals  of  the  charge  mesh,  d  ,  the  radial 
position  of  the  center  of  the  innermost  mesh  point  p(i),  zero  values  for 
the  relative  negative  charge  density  in  a  number  of  inner  "buffer"  inter¬ 
vals  n_ ( 1) . . . .n_ (j ,  unit  values  for  the  relative  negative  charge  density 

corresponding  to  the  Initial  position  of  the  beam,  and 
zero  values  for  the  relative  negative  charge  density  In  a  number  of  outer 
"buffer"  Intervals  n  (j  +1)....n  (j  )  .  The  relative  positive  charge 
density  for  ail  time  is  understood  to  have  the  same  distribution,  viz. 


n^(j)  =  0 
n^.(j)  =  i 
n^(j)  =  0 


1  ^  j  j 
J,  +  >  ^  j  ^  j. 


j  +  1  rS'  j  5  j 
•'u  ■'  •'max 


(23) 


The  initial  '•.inematic  conditions  are  the  initial  position  and  velot 
ity  of  each  macroparticle.  Since  z-motion  is  ignored  in  the  model  employed 
here,  these  consist  of  the  radial  and  azimuthal  positions  and  velocities. 

As  shewn  in  Fig.  2(b),  the  initial  radial  positions  are  simply  the  mid-radius 
of  each  macroparticle.  The  initial  radial  velocity  of  each  macroparticle  is 
assumed  to  be  zero.  The  initial  azimuthal  velocity  v  (0)  ,  specified  by  a  , 
has  bean  set  equal  to  zero  in  all  of  the  calculations  that  have  been  per¬ 
formed  . 


After  establishing  the  necessary  in'tial  conditions,  a  series  of 
j  variables  called  accumulators  representing  successive  charge  mesh 
intervals  are  set  equal  to  zero.  As  each  macroparticle  is  advanced  to  a 
new  radial  position  during  a  time  step,  the  values  of  the  appropriate  accum¬ 
ulator  variables  re  suitably  incremented.  Since  the  beam  width  does  not 
in  general  remain  constant,  the  macroparticles  may  overlap  and  so  dui  mg 
the  course  of  a  time  step  an  accumulator  may  receive  contributions  from 
more  than  one  macropart ic le . 


'■t 


The  contents  of  the  accumulators  at  the  end  of  a  time  step  become 
tlie  values  of  n  that  are  used  to  calculate  the  electrostatic  force  on 
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each  fiiacropart  i  c  le  during  the  next  time  step.  Thus  there  must  be  at  least 
two  sets  of  Such  variables  representing  the  radial  charge  distribution. 

One  set  represents  the  charge  distribution  obtained  during  the  previous 
tlrwi  stop  and  being  currently  used  to  determine  the  particle  trajectories; 
.he  second  set  represents  the  accumulators  being  incremented  during  the 
current  time  step  and  whose  values  will  be  used  to  determine  the  particle 
trajectories  during  the  next  time  step. 

In  order  to  advance  a  inacropart  icle  by  one  time  step,  the  Runge- 
Kutta  procedure  requires  that  T(x,R,S)  given  by  eq.  (21)  be  evaluated  at 
the  points  indicated  in  eqs.  {15b),  {l6b),  (17b)  and  (l8b).  The  first  three 
terms  on  the  right  hand  side  of  eq.  (25)  present  no  problem.  However,  the 
last  term  in  eq.  (21)  does  present  a  problem  due  to  the  discreteness  of  the 
charge  distribution  mesh  and  to  the  implicit  dependence  of  n  and  on 

the  time. 

The  integral  in  this  term  is  a  function  of  R  only  through  the  upper 
limit  of  integration,  ft  is  required,  therefore,  to  evaluate  the  quantity 


C(R) 


K 

% 

(n_  -  n^)R'  dR' 


{•^) 


o 

for  each  macroparticle  at  the  four  points  R  =  R.,  Rj  +  0.5A',  R.  -f  , 

and  R.  h  A'*'  at  or  near  the  radial  position  of  tl-e  i^^  macroparticle. 

To  do  this  values  of  the  discrete  variable  ?.{})  corresponding  to  C(R) 
arc  calculated  us'py  tire  djscie.;e  values  of  the  negative  and  positive  rela¬ 
tive  charge  densities  n_(j)  and  n^(j)  .  Based  on  eq.  (24)  these  are 


I 


I 


..alculated  '’rom  'ne  recursion  relations 
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C(j)  =  C(j-l)  +  [n_(j)  »  n^(j)]  p(j)  d 

p(j)  “  P(j-l)  +  d  (25) 

C(0)  =  0 

Proper  handling  of  the  implicit  time  dependence  of  the  integral  term  in 

eq.  (21)  is  also  required  in  order  to  evaluate  eqs.  (15b)-(l8b).  To  do 

this,  V  lues  of  the  integral  C(j)  for  seve'al  past  time  steps  must  be 

saved  and  used  in  an  extrapolation  from  x  to  x  +  Ax/2  and  to  x  +  Ax  . 

^  o  o  o 

In  the  work  reported  here,  values  of  C(j)  for  the  last  three  time  steps 
were  preserved  for  use  in  the  time  extrapolation. 

Advancing  each  macroparticle  by  one  time  step  consists  of  a  loop 
over  all  values  of  the  macroparticle  index  i  ,  for  each  of  which  the 
operations  indicated  by  A  in  Fig.  3  are  performed.  First  the  radial 
position  R|  of  a  macroparticle  is  located  with  respect  to  the  charge 
mesh;  that  is,  the  value  of  j  corresponding  to  the  mesh  interval  con¬ 
taining  R.  is  determined.  Next  the  nine  nearest  values  of  C(J)  at 
locations  symmet  '.aily  distributed  about  K.  are  fitted  by  a  least  squares 
polynomial.  The  degree  of  this  polynomial  is  optional;  the  result  of  using 
polynomials  of  various  orders  will  Le  discussed  later.  This  continuous  func¬ 
tion  fitted  to  the  discrete  values  of  C(j)  is  then  used  in  place  of  C(R) 
when  evaluating  the  right  hand  side  of  eq.  (21). 

In  an  early  version  of  the  program,  difficulty  was  experienced  in 
the  least  squares  fitting  due  to  an  ill-conditioned  matrix  of  coefficients 
in  the  set  of  linear  algebraic  equations  that  must  be  solve'.  This  trouble 
was  avoided  and,  at  the  same  time  computation  time  was  decreased  signifi¬ 
cantly  by  using  the  method  of  oi-thogonal  polynomials  In  the  least  squares 

(8,3) 


fitting 
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An  interpolation  between  the  known  values  of  C(j)  is  obv'Ously 
necessitated  by  the  Runge-Kutta  procedure.  The  use  of  a  least-squares 
polynomial  is  believed  to  be  superior  tc  an  n-point  interpolation  formula 
because  it  provides  a  way  to  compensate  for  the  finite  spatial  resolution 
imposed  by  the  use  of  a  digital  computer .  The  charge  distribution  can  be 
quite  discontinuous,  both  near  the  edge  of  the  beam  and  where  there  is 
only  a  small  overlap  between  macroparticles.  The  latter  situation  is 
shown  schematically  in  Fig,  4.  If  two  macroparticles  overlap  slightly, 
which  is  the  usual  situation,  a  small  step  will  develop  in  the  charge  dis¬ 
tribution.  Such  a  step  does  not  mean  that  there  would  be  a  corresponding 
sharp  discontinuity  in  the  charge  distribution,  in  an  actual  physical  device, 
however.  The  t  sential  continuity  in  nature  is  put  back  into  the  calcula¬ 
tion  by  distributing  this  step  over  the  width  of  the  macroparticle  by  means 
of  the  least-squares  fitting.  At  the  same  time  one  does  not  want  to  spread 
the  step  uniformly  because  it  is  important  to  preserve  the  gradients  of 
the  charge  distribution  since  it  is  these  that  de'jrmine  the  force  on  a 
charged  particle.  By  fitting  a  low  order  polynomial  to  such  a  distribu¬ 
tion  a  representation  of  the  radial  gradients  is  preserved  in  a  way  that 
still  allows  the  use  of  a  Runge-Kutta  integration  technique. 

In  the  light  of  this  discussion,  it  is  also  clear  that  the  number  of 
C(j)  values  used  in  the  least-squares  fit  should  be  such  as  to  approximately 
span  one  macroparticle.  Thus  the  use  of  nine  values  is  related  to  the  choice 
of  the  macroparticle  to  charge  mesh  width  ratio  D/d  ,  of  iO.  Or,  to  put 
the  matter  in  a  different  perspective,  the  choice  of  D/d  is  related  to 
the  number  of  points  and  hence  to  the  amount  of  computing  time  that  one 
wishes  to  employ  in  the  leas t -squares  fit  for  each  macroparticle  at  each 
time  step. 
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Since  the  time  variation  of  the  C{j)  quantities  is  smooth, 
a  simple  polynomial  extrapolation  procedure  could  be  used  to  evaluate 
the  integral  expression  in  eq .  (2^4)  at  the  points  in  time  required  by 
the  Runge-Kutta  procedure.  Thus,  after  locating  each  macroparticle, 
a  least-squares  polynomial  expressing  the  radial  variation  of  C{j) 
was  f'»-ted  to  each  of  the  three  sets  of  C(j)  representing  the  charge 
integral  for  the  three  past  time  steps.  Whenever  C(R)  is  required  at 
a  future  time,  the  fitted  least-squares  polynomials  for  the  three  past 
time  steps  are  evaluated  at  the  desired  radi  points  and  the  resulting 
three  values  are  used  to  extrapolate  C(R)  o  the  desired  time. 

Calling  the  value  of  C(R)  at  a  given  point  Rj^  obtained  using 

the  values  of  C(j)  for  the  most  recent  time  step  x^,CjfRj^)  ,  the  value 

of  C(R)  at  R|^  obtained  using  the  values  of  C(j)  for  the  previous 

t'f'e  step  x^  -  Ax,Cg(Rj^)  ,  and  the  value  of  C(R)  at  Rj^  obtained  using 

aiues  of  C(j)  for  the  next  previous  time  step  x^  -  a:ix,C2(R|^)  , 

then  the  value  of  C(R)  at  R,  for  x  4-  -j- Ax  is 

k  o  2 


W  =r'dV  -  5  W 

The  value  of  C(k)  at  R,  for  x  4-  Ax  is 

K  O 

C3(R,^)  =  3  Cj(R^)  -  3  C2(R^)  +  C^(R^) 

Once  the  integral  C(R)  has  been  evaluated  at  *^he  required  points 
in  space  and  time,  the  calculation  of  AR  and  AS  for  one  time  step  using 
eqs .  (15)-(18)  is  straightforward.  When  the  new  radial  position  has  been 
obtained,  the  appropriate  accumulators  must  be  incremented.  This  requires 
that  the  new  value  of  R.  be  again  located  with  respect  to  the  charge  mesh. 
The  amount  by  which  each  accumulator  is  incremented  is  such  that  the  total 
amount  of  charge  represented  by  the  macroparticle  is  conserved.  Since  each 
macroparticle  represents  a  cylindr(cal  annulus  of  mid-radius  Rj  and  fixed 
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width  D  ,  the  relative  charge  density  per  radial  interval  d  ,  n  (x,R)  , 
must  satisfy  the  relation 

n_(x,R)R| (x)  =  :  I (0)n_(0,R)  (26) 

where  R|(0)  is  the  initial  position  of  the  1^^  macropart i cie  and 
n_(0,R)  “  i  .  Thus,  if  the  width  of  a  macropart  I c ie  is,  for  computational 
convenience,  to  be  kept  fixed,  the  density  of  charge  must  be  adjusted  when 
the  macropart i c ie  departs  from  its  Initial  radial  position.  V/hen 
R|  (><)  >  Rj(0)  ,  due  to  the  cylindrical  geometry  the  macroparticle  represents 
a  greater  volume  and  so  the  charge  density  must  be  proportionally  decreased. 
The  opposite  is  the  case  when  Rj(x)  <  Rj (O)  .  Thus  the  amount  by  which 
each  accumulator  is  incremented  is  not  unity  as  at  the  beginning  of  the 
problem  but  R.(0)/R|(x)  .  Also,  since  for  x  >  0  a  macroparticle  is  not 
in  general  exactly  aligned  with  D/d  charge  mesh  intervals,  the  above 
accumulator  Increment  is  applied  to  only  the  (D/d)-l  intervals  completely 
spanned  by  the  macroparticle.  The  remaining  charge  increment  Rj(0)/R|{x) 
is  divided  between  the  inner  and  outer  accumulators  that  are  only  partially 
spanned  by  the  macropar t i c ie  in  the  ratio  of  the  overlap  of  the  two  mesh 
intervals . 

When  the  last  macroparticle  has  been  advanced  by  one  time  step,  the 
A  loop  in  Fig,  3  is  completed  and  the  calculation  enters  the  B  loop.  This 
consists  of  extracting  the  information  of  physical  interest  from  the  results, 
making  the  necessary  preparations  for  the  next  time  step  and  recycling  the 
program  to  start  the  next  time  step. 

Prior  to  these  steps,  however,  the  total  amount  of  charge  present 
in  the  accumulators  is  calculated  and  compared  with  that  present  initially. 
Due  to  the  approximate  way  in  which  the  accumulators  are  incremented  these 
will  not  in  general  be  identical.  The  contents  in  each  accumulator  are  then 


"corrected"  by  multiplying  by  the  factor  necessary  to  keep  the  amount  of 
negative  charge  in  the  beam  constant  and  thus  preserve  the  overall  neu¬ 
trality  of  the  plasma.  At  this  point  the  relative  negative  charge  dis¬ 
tribution  can  be  printed  if  desired.  This  Information  and  the  radial  and 
azimuthal  position  and  the  radial  velocity  of  each  macroparticle  are  the 
primary  output  of  the  program  for  each  time  step. 

It  is  difficult  to  form  a  clear  picture  of  the  behavior  of  the  beam 
from  this  mass  of  detailed  information.  Therefore;,  the  above  Information 
is  edited  to  yiela  t-he  mean  radial  position  of  the  negative  charge  and  the 
rms  deviation  of  the  negative  charge  from  this  mean.  These  quantities  are 
referred  to  in  the  fi  lures  to  follow  as  the  beam  position  and  the  beam  width. 
Specifically,  they  are  defined  as 


beam  -  r  - 
position 


tbeam  _ 

widthJ 


"*2  2 
n_(R')R'‘=^  dR'  E  n^(j)  pf 
J  "  *  J 

n_(R')R'  dR'  Z  n_(J)  p. 

J  “  "  J 

'n_(R')  (R'-R)^  R*  dR’  Z  n_(j)  (p.-R)^p, 

"  J  J 


(R')R'  dR' 


Z  n  (J)  p 


Finally,  the  contents  of  each  of  the  accumulators  are  transferred  to 
the  n  (j)  array,  the  time  is  incremented,  and  the  program  recycles  to  the 
point  indicated  in  Fig.  3  where  the  accumu lato-'S  are  cleared  and  next  time 
step  is  started. 

The  program  terminates  when  some  maximum  time  value  is  reached,  'ne 
finds  in  many  cases  that  the  beam  position  undergoes  a  sinusoidal  oscilla¬ 
tion  that  is  either  exponentially  growing  or  exponentially  damped.  In  these 
cases  the  really  significant  quantities  are  the  frequency,  phase,  amplitude, 
and  e-folding  rate  of  the  beam  position  oscillations.  Therefore,  the  cal¬ 
culated  beam  position  information  must  be  further  edited  to  produce  these 
quant i t ies . 


This  is  done  by  fitting  the  beam  position  as  a  function  of  time  by 
an  expression  of  the  form 

R(x)  =  Ae  sin(tjux  +  $)  +  S  (29) 

During  the  course  of  the  trajectory  calculation,  some  of  the  values  of  R 

and  X  are  saved  for  input  into  a  fina!  part  of  the  program  where  the 

five  adjustable  constants,  A,  D,  u),  and  S,  are  fitted  by  an  Iterative 
non-linear  least  squares  techn ique ^ This  has  proven  to  be  extremely 
useful  in  studying  the  validity  of  the  numerical  procedures  used  in  the 
calculation.  In  the  work  that  is  reported  here,  about  eighty  set  of 
R,x  values  spanning  about  seven  plasma  periods  have  been  used  to  determine 
the  characteristics  of  the  beam  motion.  Typically  about  thirty  iterations 
are  used  to  determine  the  five  constants.  Comparison  of  the  fitted  R 

values  with  those  used  in  the  fitting  process  shows  that  the  two  agree  to 

at  least  seven  significant  figures.  This  indicates  that  the  time  behavior 
of  R  is  very  closely  approximated  by  eq.  (29).  The  number  of  iterations 
required  to  achieve  such  accuracy  depends,  in  general,  on  the  accuracy  jf 
the  initial  gueis  used  to  start  the  iterative  process.  This  is  a  matter 
of  experience.  It  has  been  found,  for  example,  that  good  initial  values 
for  A,  Q,  and  $  are  not  important  but  that  u)  should  be  specified  to 
within  about  25  percent  of  the  correct  value  if  the  iterative  process  is 
to  converge  properly. 


RESULTS 


On  th^  basis  of  previous  single  particle  trajectory  calculations'  , 
the  effects  of  many  of  the  variables  on  the  early-time  motion  of  a  particle 
is  reasonab 1 y  we  1 1 -understood .  Therefore,  the  following  physically  realis¬ 
tic  values  have  been  chosen: 

p  -  B^/Bg  =  2.3 
u)  =  10^  sec  ' 

3  =  =  10^ 

The  choice  of  p  implies  a  vector  potential  well -depth  of  1.73  percent  for 

0=0.  The  choice  of  3  implies  a  peak  magnetic  field  in  the  0'’bit  region 
2 

of  about  5  w/m  .  As  indicated  previously,  the  Initial  velocity  of  each 
particle  was  assumed  to  be  zero  and  the  initial  position  of  the  beam  was 
chosen  to  be  relatively  narrow  and  located  symmetrically  with  respect  to 
the  minimum  of  the  vector  potential  well. 

The  first  objective  was  to  understand  the  effect  of  the  charge  dens¬ 
ity  on  the  start-up  behavior  of  the  beam.  This  involved  varying  the  value 
of  Y  .  For  all  of  these  calculations,  the  azimuthal  component  of  the  mag¬ 
netic  field  B  was  set  equal  to  zero  (i.e.  q  =  O) .  Figure  5  shows  a  set 

T 
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of  six  macropai  t  icle  trajectories  for  the  case  where  N  =  10  cm"'^  . 

o 

These  trajectories  do  not  differ  significantly  from  the  =  0  case. 

That  is,  they  are  simply  six  single-particle  trajectories;  the  radial  oscil¬ 
lations  are  the  familiar  betatron  oscillations  about  the  equilibrium  orbit 
at  R  =  1 . 14  .  The  abscissa  is  the  dimensionless  time  variable  x  ~  U)t  ; 
for  0)  =  10^  sec  ^  ,  the  dimensionless  unit  of  time  corresponds  to  one  psec. 

The  representation  of  this  Information  in  terms  of  the  beam  parameters 
position  and  width,  defined  by  eqs.  (27)  and  (28),  is  shown  in  Fig.  6. 

Because  the  initial  beam  is  located  symmetrically  about  the  equilibrium  orbit 


the  beam  position  does  not  vary  significantly  with  time.  The  beam  width 

undergoes  a  periodic  oscillation  of  lowjy  decaying  amplitude  corresponding 

to  the  adiabatic  damping  of  the  betatron  oscillations  of  the  individua! 

( i  i ) 

particles  .  The  minimum  beam  width  at  the  cross-over  points  is  not  zero 

as  implied  by  Fig.  S  because  o^  the  finite  width  of  the  macroparticies. 
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The  set  of  trajectories  in  Fig.  ,  or  =  10  cm  cieariy  shows 

the  effect  of  the  electrostatic  force  on  the  beam  behavior.  There  is  a 

tendency  for  the  concentration  of  charge  near  the  equilibrium  orbit  to 

repel  further  macroparticle  penetration.  Despite  the  apparent  disorder  of 

the  trajectories,  the  beam  position  and  beam  width  shown  in  Fig.  8  behave 

in  a  regular  manner.  The  beam  position  undergoes  largei  amplitude  oscii- 
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iations  than  those  for  =  10  c.m  (Fig.  6)  when  a  considerably  more 
symmetric  particle  motion  existed.  The  collapse  of  the  beam  still  occurs 
and  actually  the  collapse  is  more  rapid  at  the  larger  charge  density. 
Nevertheless,  the  essential  features  of  the  s i ng le-part ic ie  betatron  oscil¬ 
lations  are  still  present.  The  frequency  of  the  position  and  width  oscil¬ 
lations  corresponds  to  the  betatron  frequency.  The  frequency  of  these 
oscillations  increases  with  time  and  their  amplitude  decreases  with  time. 

When  the  case  of  a  charge  density  of  -  10*^  cm  '  is  examined, 
the  effect  of  the  electrostatic  forces  on  the  beam  motion  is  even  more 
striking.  The  trajectories  in  Fig.  9  indicate  that  the  beam  collapse  is 
considerably  inhibited.  in  addition,  there  is  now  a  high  frequeii^y  radial 
oscillation  by  each  macroparticle.  From  the  fact  that  these  radial  oscil¬ 
lations  for  e.^ch  macroparticle  are  in  phase,  it  is  clear  that  what  is 
involved  is  a  collective  oscillation  of  the  negative  charge  with  respect 
to  the  fixed  positive  charge.  The  beam  position  and  beam  width  shown  in 
Fig.  10  summarize  these  observations.  The  beam  collapse  Is  smooth  with  no 
evidence  of  betatron  oscillations.  From  the  apparent  approach  of  the 


) 
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beam  width  to  some  .iiinimum  value,  one  infers  that  ar;  equilibrium  between 
betatron  and  electrostatic  forces  has  been  reached.  One  notes,  however 
that  the  beam  position  oscillations,  while  decreasing  in  amplitude,  Still 
show  an  increase  in  frequency  with  time. 

The  period  of  the  initial  radial  beam  position  oscillation  in 
Fig.  10  is  found  to  be  very  nearly  the  plasma  frequency  defined  by  eq.  (19) • 
Thus,  the  effect  of  the  electrostatic  space  charge  force  on  the  early-time 
motion  of  a  finite  width  beam  in  a  plasma  betatron  is  the  following.  The 
rapid  collapse  of  the  beam  to  the  single-particle  equilibrium  orbit  and 
the  subsequent  ad i abat i ca 1 1 y-damped  radial  betatron  oscillations  are 
replaced  by  an  inhibited  spatial  collapse  of  the  beam  and  by  a  radial  oscil¬ 
lation  of  the  position  of  the  beam  with  respect  to  the  fixed  positive  charge 
at  the  plasma  frequency.  The  inhibited  collapse  the  beam  was  expected; 
the  appearance  of  a  radial  charge  oscillation  at  the  plasma  frequency  is 
not  unreasonable  either.  The  rising  magnetic  field  excites  rad  al  betatron 
(single  particle)  oscillations  which  in  turn  drive  the  radial  plasma 
(collective)  oscillations.  The  energy  source  is  the  time-varying  magnetic 
field.  When  the  plasma  frequency  is  much  higher  than  the  radial  betatron 
frequency  one  expects  the  plasma  oscillations  to  occur  essentially  at  the 

'“free"  collective  mode  frequency  which  is  very  nearly  the  plasma  frequency. 
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When  =10  cm  ,  the  charge  density  is  sufficiently  great  that 
the  collapse  of  the  beam  is  completely  inhibited,  as  shown  by  the  macro¬ 
particle  trajectories  in  Fig.  11.  In  this  case  the  motion  is  completely 
characterized  by  the  radial  plasma  oscillations.  The  beam  position  and 
beam  width  shown  in  Fig.  12  confirm  this  observation.  Note  that  for  t-hese 
larger  values  of  the  charge  density  when  the  motion  is  dominated  by  elec¬ 
trostatic  plasma  oscillations,  it  is  necessary  to  change  the  time  scale  of 
the  calculation  with  the  plasma  period.  Thus,  between  Figs.  10  and  12 
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there  is  a  factor  of  100  difference  in  the  charge  density.  This  means 
that  the  plasma  frequency  i higiier  by  a  factor  of  10  in  Fig.  12  and  so 
the  total  amount  of  time  shown  in  Figs.  15  and  12  is  less  by  a  factor  of 
iO  than  that  shown  in  Figs.  $  and  10.  For  the  lower  values  of  charge 
density,  the  plasma  frequency  Is  less  than  the  betatron  frequency  and  so 
the  time  scale  of  the  calculation,  which  is  then  determined  by  the  beta¬ 
tron  periodj,  does  not  depend  on  the  charge  density. 

Several  supplementary  investigations  were  made  to  assess  the 
validity  of  the  computed  results.  One  of  these  consisted  of  solving  a 
series  of  otherwise  identical  problems  where  the  degree  s  of  the  least- 
squares  polynomial  used  to  interpolate  between  the  available  values  of 
C(j)  wa'-  varied  over  the  range  1  s;  s  ^  3  .  It  was  found  that  the  com¬ 
puted  results  were  not  dependent  on  s  .  Aside  from  the  practical  matters 
of  computing  time  and  the  importance  of  truncation  error  in  the  determina¬ 
tion  of  the  coefficients  of  the  least -squares  polynomials,  there  is  a 
logical  necessity  in  keeping  s  small.  As  discussed  previously  the  pur¬ 
pose  of  the  least -squares  interpolation  is  to  smooth  an  otherv<ise  irregular 
charge  distribution;  the  use  of  a  high  degree  polynomial  tends  to  defeat 
this  objective  since  the  interpolation  function  then  beconries  an  increasingly 
"better"  fit  to  the  irregularities  one  desires  to  smooth.  For  all  of  the 
calculations  presented  here,  s  was  chosen  to  be  two.  It  was  felt  that 
this  was  a  reasonable  compromise,  smoothing  C(j)  and  still  preserving 
its  radial  dependence. 

The  effect  of  the  number  of  significant  figures  carried  by  the  com¬ 
puter  was  examined  by  repeating  some  typical  calculations  on  a  computer  that 
worked  to  sixteen  figure  accuracy  instead  of  the  eight  figure  dccuracy  of 
the  machine  used  for  all  of  the  other  calculations  reported  here.  This 
could  be  a  serious  effect  since  as  shown  by  eq .  (25)  the  calculation  of 
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t(j)  requires  taking  the  difference  between  n  and  n^  .  Since  often 
these  are  both  near  unity,  there  can  be  a  loss  of  three  or  four  signifi¬ 
cant  figures  at  this  point.  No  significant  effect  on  such  integrated 
quantities  as  the  beam  parameters  defined  by  eqs.  (27)  and  (28)  was 
observed;  however  such  details  as  the  macroparticle  trajectories  and  the 
charge  distributions  showed  some  small  but  perceptible  differences. 

A  most  important  supplementary  investigation  was  that  of  the  effect 
of  the  time  step  on  the  calculation.  Th.e  earliest  calculations  of  collec¬ 
tive  effects  were  required,  because  of  computer  limitations,  to  ignore  the 
time  dependence  of  the  charge  integral  C(*)  .  it  was  found  that  the  elec¬ 
tron  motion  was  quite  unstable.  Within  a  few  plasma  periods  the  beam 
position  oscillations  became  quite  large;  coupled  unstable  oscillations  in 
the  beam  width  occurred  also.  Eventually,  coherent  beam  motion  ceased  as 
macroparticles  broke  away  from  the  positive  charge  entirely.  However,  it 
was  observed  that  the  magnitude  of  the  growth  rate  of  the  instability 
decreased  as  the  time  step  decreased  indicating  that  the  observed  beam 
instability  had  a  numerical  rather  than  s  physical  origin. 

It  was  at  this  point  that  time  was  made  available  on  a  computer  that 
offered  more  memory  capacity  and  greater  speed.  Not  only  could  the  time 
dependence  of  C(j)  be  treated  but  also  the  time  step  could  be  made  much 
smaller.  The  first  step  was  to  consider  a  linear  extrapolation  of  C(J) 
using  the  results  of  the  previous  two  time  steps.  It  was  found  that  this 
procedure  had  an  enormous  effect  in  decreasing  the  growth  rate  of  the 
instability.  As  before,  it  was  found  that  decreasing  the  time  step  also 
decreased  the  growth  rate. 

However,  the  beam  position  oscillation  was  still  unstable  although 
for  a  sufficiently  small  time  step  the  effect  was  barely  perceptible.  The 
next  step  was  to  perform  a  quadratic  extrapolation  of  C(j)  using  the 


results  of  the  previous  three  time  steps.  When  this  was  done  It  was 
found  that  the  beam  position  oscillation  v^as  no  longer  unstable.  The 


coefficient  of  the  exponential  factor  sn  eq.  (29)  was  found  to  be  nega¬ 
tive  indicating  that  the  plasma  oscillations  were  damped.  Again  it  was 
found  that  the  rate  of  this  process  decreased  as  the  time  step  decreased. 
Therefore,  thr  interpretation  has  been  made  on  the  basis  of  both  the 
linear  and  quadratic  extrapolation  of  C(j)  thai  the  observed  plasma 
oscillations  are  stable  in  the  limit  as  the  time  step  goes  to  zero. 

The  residual  growth  or  decay  rate  for  a  small  finite  time  step  is  believed 
to  be  a  purely  numerical  effect.  These  results  are  summarized  in  Table  I. 
For  all  of  the  calculations  reported  here,  quadratic  extrapolation  of  C(j) 
was  employed  and  the  time  stap  was  either  5  x  10  ^  or  1/220  of  the  plasma 
period  if  that  was  smaller. 

Finally,  the  effect  of  a  B  field  on  the  results  was  studied. 

1 2  -3 

A  charge  density  =  10  cm  was  selected  since  this  represents  a 

case  v;here  the  electrostatic  forces  completely  dominate  the  bean  oehavior, 

-2 

When  q  =  B,/B_  =  10  ,  corresponding  to  B_(r,)  =  5^0  gauss,  there  was 

virtually  no  effect  on  the  radial  position  oscillation.  The  beam  width 
remained  constant  at  its  initial  value  as  shown  in  Fig,  12  for  q  =  0  . 
However,  when  q  =  1  i,e,  when  the  B^  field  is  essentially  equal  to 
the  betatron  field  in  the  orbit  region,  the  radial  plasma  oscillations 
were  completely  suppressed.  The  macroparticle  radial  positions  were 
constant  in  time  and  the  beam  width  was  consta.it,  therefore. 


Table  I 


Dependence  of  time 

step  and  time-dependence  of 

12  -3 

Tor  N  =  10  cm 

0 

C(J)  on  results 

time-dependence 

plasma  period 

e-foldinq  time 

of  C(j) 

time  step 

plasma  period 

none 

11 

1.2 

22 

2.3 

44 

4.5 

i i near 

22 

41 

44 

340 

no 

5,600 

220 

38,000 

440 

150,000 

quadratic 

22 

-4o 

no 

-2,500 

220 

-27,000 

2( 


V  CONCLUSIONS 

On  the  basis  or  the  results  reported  here^  the  following  conclusions 
can  be  drawn: 

9-3 

1.  For  charge  densities  less  than  about  10  cm  ,  the  electrostatic  field 
that  develops  during  the  start-up  of  a  plasma  betatron  is  ineffective  in 
preventing  the  collapse  of  the  beam  to  the  single-particle  equilibrium 
orbit 

2.  For  charge  densities  in  the  range  10^  -  10^^  cm  the  beam  behavior 
Indicates  a  transition  from  that  of  a  single-particle  pure  betatron 
character  to  a  collective  character  dominated  by  long-range  electrostatic 
effects 

1 0  -3 

3.  For  charge  densities  greater  than  about  10  cm  ,  a  radial  plasma 
oscillation  appears.  This  consists  of  an  oscillation  of  the  center  of 
charge  of  the  electrons  with  respect  to  the  relatively  fixed  positive 
charge.  The  frequency  of  this  charge  oscillation  i':  the  plasma  frequency 
k.  This  radial  plasma  oscillation  appears  to  be  stable  if  the  space  and 
time  dependence  of  the  charge  integral  term  is  handled  sufficiently  care¬ 
fully  aiid  if  the  time  step  is  sufficiently  small 

5.  The  addition  of  an  azimuthal  component  of  the  magnetic  field  compar¬ 
able  to  the  betatron  field  at  the  orbit  suppresses  the  radial  plasma 
oscillation  completely 
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13  ABSTRACT 


Runaway  electrons  frenerated  by  induction  acceleration  in  a  plasma  betatron  with 
axially  symmetric  guide  field  B^(r,z,t)=B^  (r^/r)"x  f(g)  sin  2nt/T  (n=0.6,T=l6usec , 
r=4,8cm)  together  w^tn  auxilliary  azimuthal  guide  field  R, (r )®(k/r) ,{2  kG<B  (r)< 

5kG)  have  been  studied  for  a  wide  range  of  induction  eleccric  fields  Eg(27$  v/cra< 
(E)g^26v/cm) .  Experiments  have  been  carried  out  in  argon  and  krypton  at  ambient 
pressures  of  0.5-r  Torr.  The  ratio  of  the  applied  electric  field  to  the  critical 
electrical  field  for  runaway  was  typically  '^'500  so  that  most  of  the  electrons  shoul 
run  away.  Runaway  currents  of  ''-lA  were  observed  immersed  in  a  background  conductio 
current  of  25-100  A.  Thi“  runaway  electrons  strike  the  vacuum  chsmiber  walls  at  a  ti 
t^  after  betatron  acceleration  has  begun.  The  resulting  x-rays  are  detected  with  a 
scintillation  counter.  The  x-ray  time  t  hae  been  measured  as  a  function  of  E  and 
Bgfor  several  half-cycles  of  the  guide  f?eld.  It  is  found  that  t  decreases  with  , 
E  but  increases  with  B^.  It  is  also  found  that  t  depends  on  the  relative  oriental 
tion  of  Eg  with  respect  to  B.,  This  is  shown  to  be  due  tc  a  small  transverse  mag¬ 
netic  field  associated  with  the  Bg  guide  field.  The  small  runaway  current  is  con¬ 
sistent  with  nat  determined  from  the  self  electric  field  generated  by  the  adiabatic 
constriction  of  the  runaway  stream  or  the  Negative-Mass-Instability,  The  reason  for 
the  beam  disruption  and  subsequent  x-ray  emission  are  unexplained  but  the  mechanism 
of  Field  and  Fr..ed  appears  to  describe  the  process. 

The  betatron  guide  field  was  generated  in  a  low  inductance  (50  x  10'"'^h)  single¬ 
turn  coil  using  flux-concentrators .  Techniques  used  to  construct  the  flux-concentrh 
tors  and  measure  the  resulting  pulsed  magnetic  field  distribution  are  also  described 
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